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Abstract 

We review how nuclear forces emerge from low-energy QCD via chiral effective field theory. The presentation 
is accessible to the non-specialist. At the same time, we also provide considerable detailed information 
y—{ (mostly in appendices) for the benefit of researchers who wish to start working in this field. 
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1. Introduction and historical perspective 

The theory of nuclear forces has a long history (cf. Table [T]). Based upon the seminal idea by Yukawa [T], 
first field-theoretic attempts to derive the nucleon-nucleon (NN) interaction focused on pion exchange. 
While the one-pion exchange turned out to be very useful in explaining NN scattering data and, in particular, 
the properties of the deuteron [2 , multi-pion exchange was beset with serious ambiguities ^ |4j that could 
not be resolved in a satisfactory way. Thus, the "pion theories" of the 1950s are generally judged as failures — 
for reasons we understand today: pion dynamics is constrained by chiral symmetry, a crucial point that was 
unknown in the 1950s. 

Historically, the experimental discovery of heavy mesons [S] in the early 1960s saved the situation. The 
one-boson-exchange (OBE) model [BHTTj emerged, which still today is the most economical and quantitative 
phenomenology for describing the NN interaction [T2HI1]- The weak point of this model, however, is the 
scalar-isoscalar "sigma" or "epsilon" boson, for which empirical evidence remains controversial. Since this 
boson is associated with the correlated (or resonant) exchange of two pions, a vast theoretical effort was 
launched to derive the 27r-exchange contribution of the nuclear force, which creates the intermediate-range 
attraction. During this effort, which occupied more than a decade, dispersion theory (Stony Brook [l5 l fTB ] 
and Paris [ITj [18] potentials) as well as field theory (Partovi-Lomon model [19], Bonn potential [TTl I20j) 
were invoked. 

The nuclear force problem appeared to be solved; however, with the discovery of quantum chromo- 
dynamics (QCD), all "meson theories" were relegated to the status of models and the attempts to derive 
the nuclear force had to start all over again. 

The problem with a derivation of nuclear forces from QCD is that this theory is non-perturbative in the 
low-energy regime characteristic of nuclear physics, which makes direct solutions very difficult. Therefore, 
during the first round of new attempts, QCD-inspired quark models pTH5T] became popular. The positive 
aspect of these models is that they try to explain hadron structure and hadron-hadron interactions on 
an equal footing and, indeed, some of the gross features of the NN interaction are explained successfully. 
However, on a critical note, it must be pointed out that these quark-based approaches are nothing but 
another set of models and, thus, do not represent fundamental progress. For the purpose of describing 
hadron-hadron interactions, one may equally well stay with the simpler and much more quantitative meson 
models. 

A major breakthrough occurred when the concept of an effective field theory (EFT) was introduced and 
applied to low-energy QCD. As outlined by Weinberg in a seminal paper [32], one has to write down the 
most general Lagrangian consistent with the assumed symmetry principles, particularly the (broken) chiral 
symmetry of QCD. At low energy, the effective degrees of freedom are pions (the Goldstone bosons of the 
broken symmetry) and nucleons rather than quarks and gluons; heavy mesons and nucleon resonances are 
"integrated out". So, the circle of history is closing and we are back to Yukawa's meson (pion) theory, 
except that we have finally learned how to deal with it: broken chiral symmetry is a crucial constraint that 
generates and controls the dynamics and establishes a clear connection with the underlying theory, QCD. 

The idea of chiral symmetry has an interesting history of its own. The modern understanding is that 
this symmetry arises because the up and down quarks happen to have relatively small masses. However, 
chiral symmetry and its significance for low-energy hadron (pion) physics was discovered long before QCD. 
In 1960, based upon concepts proposed by SchwingerJ53] , Gell-Mann and Levy [33] developed the sigma 
model, which is a linear realization of chiral symmetryjj One major problem researchers had been struggling 
with in the 1950 's was that the pion- nucleon scattering length came out two orders of magnitude too large 
when the (renormalizable) pseudo-scalar (75) TriV interaction was used. This unrealistic prediction was 
due to very large contributions from virtual anti-nucleon states (the so-called "pair terms" or "Z-graphs"). 
Similar problems occurred in the 27r-exchange contribution to the NN interaction. In the sigma model, the 
large pair terms are canceled by processes involving the (fictitious) a boson. In this way, the linear sigma 
model demonstrates how imposing chiral invariance fixes the problem with low-energy n-N scattering. 



^For a pedagogical introduction into chiral symmetry and the sigma model, see |35l . 
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Table 1: The Theory of Nuclear Forces: Seven Decades of Struggle 



1935 


Yukawa: Meson Theory 


1950's 


The "Pion Theories" . 




One-pion exchange: good; multi-pion exchange: disaster. 


1960's 


Many pions = muhi-pion resonances: cr, p, a;, ... 




The One-Boson-Exchange Modeh success. 


1970's 


Divers two-pion-exchange models: 




Partovi-Lomon, Stony Brook, Paris, Bonn. 


1980's 


Nuclear physicists discover QCD: 




Quark Models. 


1990's 


Nuclear physicists discover EFT; Weinberg, van Kolck, . . . 


and beyond 


Back to Meson (Pion) Theory! But, constrained by Chiral Symmetry. 



However, the fictitious character of the a particle as well as the reliance on the perfect cancelation of huge 
terms are uncomfortable features. In 1967, motivated by the current algebra approach to soft pion physics, 
Weinberg [36] worked out what has become known as the non-linear sigma model, which does not include 
a a anymore and has pions and nucleons interact via pseudo-vector (derivative, 757^9^) coupling besides 
a new (non- linear) irirNN term also involving a derivative ( "Weinberg- Tomozawa term" [37l |38]). The 
derivative (equivalent to momentum) guarantees that the interaction vanishes when the momentum goes to 
zero providing a natural explanation for the weakness of the interaction by soft pions which does not rely 
on the cancelation of large terms. Following suggestions by Schwinger, Weinberg [33] developed, soon after, 
a general theory of non-linear realizations of chiral symmetry, which was further generalized in an elegant 
way by Callan, Coleman, Wess, and Zumino [40] . 

Even though the original work on chiral symmetry was obviously all performed by particle physicists, it 
must be stated — to the honor of nuclear physics — that there have been some far-sighted nuclear physicists 
who early on understood and appreciated the significance of chiral symmetry for low-energy hadron inter- 
actions. One of them was Gerry Brown, who as early as 1968 published with two co-workers a paper [IT] 
on three-nucleon forces, where the consequences of chiral symmetry are fully exploited. In 1970, Brown 
wrote a remarkable Comment [42] and, in 1979, he published a book chapter entitled "Chiral symmetry and 
the nucleon-nucleon interaction" [33] . Moreover, in the more sophisticated relativistic meson models of the 
past [ini HiHlS] the pseudo-vector coupling was applied in ttNN vertices (instead of the simpler pseudo- 
scalar one that was commonly in use) in recognition of chiral symmetry. However, this chiral patch work, 
even though it points into the right direction, cannot be perceived as a serious chirally invariant theory. 
Moreover, one has to face the problem that the derivative coupling is not renormalizable in the conventional 
sense. 

Therefore, ideas were still needed for how to implement chiral symmetry consistently in the theory of 
pionic and nuclear interactions and how to deal with the renormalization issue. In his contribution to the 
'Festschrift' in honor of Schwinger of 1979 [35] H?]) Weinberg proposed to consider the most general possible 
Lagrangian including all higher-derivative terms that are consistent with chiral symmetry (besides the other 
commonly accepted symmetry principles) . For this theory to be manageable, one needs to assume some sort 
of perturbative expansion such that only a finite number of terms contribute at a given order. This expansion 
is provided by powers of small external momenta over the chiral symmetry breaking scale, A-^ 1 GeV. 
The higher-derivative terms supply the counterterms that make possible an order-by-order renormalization, 
which is the appropriate renormalization procedure for an effective field theory. Weinberg's suggestions were 
soon picked up by Gasser, Leutwyler, and associates who worked out, to one loop, the cases of wtt [48] and 
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■kN scattering with great success. 

But there was stiU the problem of the nuclear force which is more difficult, since nuclear interactions 
do not vanish in the chiral limit {q — > 0; rriu/di'mTT — > 0) and require a non-perturbative treatment because 
of the existence of nuclear bound states. In a series of papers published around 1990 [SHHS^ . Weinberg 
picked up the nuclear force issue and suggested to calculate the NN potential perturbatively in the chiral 
expansion and then iterate it to all orders in a Schroedinger or Lippmann-Schwinger equation to obtain the 
nuclear amplitude. Here, the introduction of four-nucleon contact terms is crucial for renormalization. 

Following the Weinberg proposal, pioneering work was performed by Ordonez, Ray, and van Kolck |53[I54) 
who applied time-ordered perturbation theory to construct a NN potential up to next-to-next-to-leading 
order (NNLO). The results were encouraging and nuclear EFT quickly developed into one of the most popular 
branches of modern nuclear physics. The Munich group used covariant perturbation theory and dimensional 
regularization to calculate the perturbative NN amplitude without [5S] and with A(1232)-isobar degrees of 
freedom [56] at NNLO. Besides this, the Munich group worked out important loop contributions of higher 
order [571 - I6T] . A relativistic approach was also taken by the Brazil group [62l |63]. The Bochum-Jiilich 
group devised a method of unitarity transformations to eliminate the energy-dependence of time-ordered 
perturbation theory amplitudes and calculated the NN potentials up to NNLO [Ml [SS] . The Idaho group 
managed to construct a chiral NN potential at next-to-next-to-next-to-leading order (N'^LO) and showed 
that only at this order can one achieve the precision necessary for reliable few-nucleon and nuclear structure 
calculations [661 - 169] . Progress extended beyond the NN interaction, as nuclear many-body forces based 
upon chiral perturbation theory were also developed [EH I7D1175] . 

During the past decade or so, chiral two-nucleon forces have been used in many microscopic calculations 
of nuclear reactions and structure [7iH55] and the combination of chiral two- and three-nucleon forces has 
been applied in few-nucleon reactions [TT] I55H^ , structure of light- and medium-mass nuclei 'M UlOlj , and 
nuclear and neutron matter |1021 1103] — with a great deal of success. The majority of nuclear structure 
calculations is nowadays based upon chiral forces. 

Consequently, it may be of interest to the community to have a good understanding of these forces and 
their background. It is therefore the purpose of this report to provide an accessible review on how nuclear 
forces emerge from low-energy QCD via chiral effective field theory. A pedagogical introduction into the 
phenomenology and the traditional view of nuclear forces can be found in Refs. [TT lll04j . Alternative reviews 
on various aspects of the modern perspective are published in Refs. [ 105H108] . 

This article is organized as follows. In Section [2] we sketch the foundations of an EFT for low-energy 
QCD including the effective Lagrangians. Section [3] provides an overview on nuclear forces derived from 
chiral EFT. The two-nucleon force is then discussed in detail in Section[4]and many-body forces in SectionjSj 
The extension of the theory through the introduction of A(1232)-isobar degrees of freedom and higher order 
contributions to nuclear forces are considered in Section |6] Finally, Section [7] contains our conclusions. The 
appendices provide many mathematical details which may be useful to researchers who wish to start working 
in the field. 

2. Effective field theory for low-energy QCD 

Quantum chromodynamics (QCD) is the theory of strong interactions. It deals with quarks, gluons and 
their interactions and is part of the Standard Model of Particle Physics. QCD is a non-Abelian gauge field 
theory with color SU{3) the underlying gauge group. The non-Abelian nature of the theory has dramatic 
consequences. While the interaction between colored objects is weak at short distances or high momentum 
transfer ("asymptotic freedom"); it is strong at long distances ( ^ 1 fm) or low energies, leading to the 
confinement of quarks into colorless objects, the hadrons. Consequently, QCD allows for a perturbative 
analysis at large energies, whereas it is highly non-perturbative in the low-energy regime. Nuclear physics 
resides at low energies and the force between nucleons is a residual color interaction similar to the van der 
Waals force between neutral molecules. Therefore, in terms of quarks and gluons, the nuclear force is a 
very complicated problem that, nevertheless, can be attacked with brute computing power on a discretized, 
Euclidean space-time lattice (known as lattice QCD). In a recent study [109] . the neutron-proton scattering 
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lengths in the singlet and triplet S"- waves have been determined in fully dynamical lattice QCD, with a 
smallest pion mass of 354 MeV. This result is then extrapolated to the physical pion mass with the help of 
chiral perturbation theory. The pion mass of 354 MeV is still too large to allow for reliable extrapolations, 
but the feasibility has been demonstrated and more progress can be expected for the near future. In a lattice 
calculation of a very different kind, the nucleon-nucleon {NN) potential was studied |110l llll"j . The NN 
potential is extracted from the equal-time Bethe-Salpeter amplitude with local interpolating operators for 
the nucleons. The central part of the potential shows a repulsive core plus attraction of intermediate range. 
This is a very promising result, but it must be noted that also in this investigation still rather large pion 
masses are being used. In any case, advanced lattice QCD calculations are under way and continuously 
improved. However, since these calculations are very time-consuming and expensive, they can only be used 
to check a few representative key-issues. For everyday nuclear structure physics, a more efficient approach 
is needed. 

The efhcient approach is an effective field theory. For the development of an EFT, it is crucial to identify 
a separation of scales. In the hadron spectrum, a large gap between the masses of the pions and the masses 
of the vector mesons, like p(770) and a;(782), can clearly be identified. Thus, it is natural to assume that 
the pion mass sets the soft scale, Q ^ mjr, and the rho mass the hard scale, ^ irip, also known as the 
chiral-symmetry breaking scale. This is suggestive of considering an expansion in terms of the soft scale 
over the hard scale, Q/A^. Concerning the relevant degrees of freedom, we noticed already that, for the 
ground state and the low-energy excitation spectrum of an atomic nucleus as well as for conventional nuclear 
reactions, quarks and gluons are ineffective degrees of freedom, while nucleons and pions are the appropriate 
ones. To make sure that this EFT is not just another phenomenology, it must have a firm link with QCD. 
The link is established by having the EFT observe all relevant symmetries of the underlying theory. This 
requirement is based upon a 'folk theorem' by Weinberg •32]: 

If one writes down the most general possible Lagrangian, including all terms consistent with 
assumed symmetry principles, and then calculates matrix elements with this Lagrangian to any 
given order of perturbation theory, the result will simply be the most general possible S-matrix 
consistent with analyticity, perturbative unitarity, cluster decomposition, and the assumed sym- 
metry principles. 

In summary, the EFT program consists of the following steps: 

1. Identify the soft and hard scales, and the degrees of freedom appropriate for (low-energy) nuclear 
physics. 

2. Identify the relevant symmetries of low-energy QCD and investigate if and how they are broken. 

3. Construct the most general Lagrangian consistent with those symmetries and symmetry breakings. 

4. Design an organizational scheme that can distinguish between more and less important contributions; 
a low-momentum expansion. 

5. Guided by the expansion, calculate Feynman diagrams for the problem under consideration to the 
desired accuracy. 

In the following (sub)sections, we will elaborate on these steps, one by one. Since we discussed the first step 
already, we will address now step two. 

2.1. Symmetries of low- energy QCD 

In this section, we will give a brief introduction into (low-energy) QCD, its symmetries and symmetry 
breakings. More detailed presentations of this topic are provided in Refs. |112H114] . 

2.1.1. Chiral symmetry 

The QCD Lagrangian reads 

>Cqcd - q{irV, - M)q - \g^.,aGr (2.1) 
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with the gauge-covariant derivative 

- 5p - *5y^p,a (2.2) 

and the gluon field strength tensoi[^ 

In the above, q denotes the quark fields and Ai the quark mass matrix. Further, g is the strong coupling 
constant and Afj_,a are the gluon fields. The arc the Gell-Mann matrices and the fate the structure 
constants of the SU (3)coior Lie algebra (a, 6, c = 1, . . . , 8); summation over repeated indices is always implied. 
The gluon-gluon term in the last equation arises from the non-Abelian nature of the gauge theory and is 
the reason for the peculiar features of the color force. 

The masses of the up (u), down (d), and strange (s) quarks are |115j : 

mu = 2.5±0.8MeV, (2.4) 
md = 5 ± 0.9 MeV, (2.5) 
TTis = 101 ± 25 MeV. (2.6) 

These masses are small as compared to a typical hadronic scale, i.e., a scale of low- mass hadrons which are 
not Goldstone bosons, e.g., rup = 0.78 GeV ~ 1 GeV. 

It is therefore of interest to discuss the QCD Lagrangian in the limit of vanishing quark masses: 

^°QCD = QH^'V^q - J^^..at/r • (2.7) 
Defining right- and left-handed quark fields, 

qR = PRq, QL^PLq, (2.8) 

with 

^^fl. = ^(1+75), Pl = ^(1-75), (2.9) 
we can rewrite the Lagrangian as follows: 

-Cqcd = QRil'^T^t.qR + QLil'^V^qL - \gt.uMgr . (2.10) 

Restricting ourselves now to up and down quarks, we see that Cqqy) invariant under the global unitary 
transformations 

,«^(-)^»««^cxp(-,ef5)(-) ,2.11, 



and 

Ql 



(-)-».,. = =p(-.ef|)(-), (2.12, 

where {i = 1,2,3) are the generators of S'C/(2)flavor, the usual Pauli spin matrices with commutation 
relations 

.2 2 J 2 ^ ' 

and gji and gL are elements of SU{2)fi and SU{2)l, respectively. In conclusion: The right- and left-handed 
components of massless quarks do not mix. This is SU{2)ii x SU(2)l symmetry, also known as chiral 
symmetry. Noether's Theorem implies the existence of six conserved currents; three right-handed currents 

R-^QrY^Qr with a^i?f=0 (2.14) 



^For SU (N) group indices, we use Latin letters, . . . , a,b, c, . . . ,i, j, k, . . ., and, in general, do not distinguish between sub- 
scripts and superscripts. 
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and three left-handed currents 

Lt^qLl^'^QL with d^Ll = 0. (2.15) 
It is useful to consider the following linear combinations; namely, three vector currents 

= + Lf = g7^|(? with d^Vf^O (2.16) 

and three axial-vector currents 

= if =97^75 with d^Al^O, (2.17) 

which got their names from the fact that they transform under parity as vector and axial-vector current 
densities, respectively. The vector transformations are given by 

«)^cxp(-.6r|)(;;), ,2.18, 

which are isospin rotations and, therefore, invariance under vector transformations can be identified with 
isospin symmetry. 

There are six conserved charges, 

(fx B!1 = j (fx q^j^{t,x)'^qR{t,x) with = (2.19) 

and ^ 

^ j (fxLl= j (fxq{{t,x)'^qL{t,x) with "^=0, (2.20) 

QY = / (fx = / (fx q\t,x)^q{t,x) with ^ = (2.21) 

(fxA'l^j(fxq\t,x)-i5'^q{t,x) with ^=0. (2.22) 
The Qf and Qf satisfy the commutation relations of the Lie algebra of SU{2)jj x SU{2)]^ (chiral algebra), 

[Qf , Qf] = je»jfcQ« [gf , gj^] = *6»^"^g^, [gf , Qf] = o . (2.23) 

For QY and Qf, the commutation relations read, 

[gr , gj] = ze^^-^-gr , [Qf, Qf] = ^^'''qL [qY, Qf] - ^^'''Qt ■ (2.24) 

For reasons of completeness, we mention that massless u and d quarks satisfy an even larger symmetry 
group, namely, S'C/(2)/j x SU{2)l x U{1)v x U{\)a- While the U{l)v symmetry corresponds to quark 
number conservation, the U{\)a is broken on the quantum level ("J7(1)a anomaly") and is not a symmetry 
of the system. 

2.1.2. Explicit symmetry breaking 

The mass term —qAiq in the QCD Lagrangian Eq. (2.1 1 breaks chiral symmetry explicitly. To better 
see this, let's rewrite Ai for the two-flavor case. 



or, alternatively. 



and 



M 



niu 
rrid 

= ^[mu + md){Q ^ j + :^[mu - md) { Q _^ 
= ^{mu + rrid) I + ^{rriu - nid) T3 . (2.25) 



The first term in the last equation in invariant under SU{2)v (isospin symmetry) and the second term 



vanishes for = irid- Thus, isospin is an exact symmetry if m„ = nid- However, both terms in Eq. (2.251 



break chiral symmetry. Since the up and down quark masses [Eqs. ( |2.4[ ) and (2.5)] are small as compared to 
the typical hadronic mass scale of ~ 1 GeV, the explicit chiral symmetry breaking due to non- vanishing quark 
masses is very small. Note also that, as a consequence of the non- vanishing quark masses, the axial- vector 



current, Eq. (2.17), is not conserved anymore. 



2.1.3. Spontaneous symmetry breaking 

A (continuous) symmetry is said to be spontaneously broken if a symmetry of the Lagrangian is not 
realized in the ground state of the system. There is evidence that the (approximate) chiral symmetry of the 
QCD Lagrangian is spontaneously broken — for dynamical reasons of nonperturbative origin which are not 
fully understood at this time. The most plausible evidence comes from the hadron spectrum. 



Since the conserved quantity Qf, Eq. (2.22), commutes with the Hamiltonian and has negative parity, 
one naively expects the existence of degenerate hadron multiplets of opposite parity, i.e., for any hadron 
of positive parity one would expect a degenerate hadron state of negative parity and vice versa. However, 
these "parity doublets" are not observed in nature. For example, take the p-meson which is a vector meson 
of negative parity (J^ — 1^) and mass 776 MeV. There does exist a l"*" meson, the ai, but it has a mass of 
1230 MeV and, therefore, cannot be perceived as degenerate with the p. On the other hand, the p meson 
comes in three charge states (equivalent to three isospin states), the p^ and the p", with masses that differ 
by at most a few MeV. Thus, in the hadron spectrum, SU{2)v (isospin) symmetry is well observed, while 
axial symmetry is broken: SU{2)^ x SU{2)]^ is broken down to SU{2)y. As a consequence of this, the 
vacuum (QCD ground state) is invariant under vector transformations, i. e., QY\0) — 0, while this is not 
the case for axial transformations, Qf\0) ^ 0, where |0) denotes the vacuum. 

A spontaneously broken global symmetry implies the existence of (massless) Goldstone bosons wit h the 
quantum numbers of the broken generators |1161 1117j . The broken generators are the Qf of Eq. (2.22 1 



which are pseudoscalar. The Goldstone bosons are identified with the isospin triplet of the (pseudoscalar) 
pions, which explains why pions are so light. The pion masses are not exactly zero because the up and down 
quark masses are not exactly zero either (explicit symmetry breaking). Thus, pions are a truly remarkable 
species: they reflect spontaneous as well as explicit symmetry breaking. Goldstone bosons interact weakly 
at low energy. They are degenerate with the vacuum and, therefore, interactions between them must vanish 
at zero momentum and in the chiral limit (m-,^ — > 0). 



2.2. Chiral effective Lagrangians 
2.2.1. Relativistic formulation 

The next step in our EFT program is to build the most general Lagrangian consistent with the (broken) 
symmetries discussed above. An elegant formalism for the construction of such Lagrangians was developed 
by Callan, Coleman, Wess, and Zumino (CCWZ) [40 who worked out the group-theoretical foundations 
of non-linear realizations of chiral symmetryj^ It is characteristic for these non-linear realizations that, 
whenever functions of the Goldstone bosons appear in the Lagrangian, they are always accompanied with 
at least one space-time derivative. The Lagrangians given below are built upon the CCWZ formalism. 

As discussed, the relevant degrees of freedom are pions (Goldstone bosons) and nucleons. Since the 
interactions of Goldstone bosons must vanish at zero momentum transfer and in the chiral limit (m-^ — >■ 0), 
the low-energy expansion of the Lagrangian is arranged in powers of derivatives and pion masses. The hard 
scale is the chiral-symmetry breaking scale. Ay, 1 GeV. Thus, the expansion is in terms of powers of Q/A^ 
where Q is a (small) momentum or pion mass. This is chiral perturbation theory (ChPT). 

The effective Lagrangian can formally be written as, 

^off = -Ctttt + ^TfAf + ■ • ■ , (2.26) 



'An accessible introduction into the rather involved CCWZ formalism can be found in Ref. I113| . 
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where ^Ctttt deals with the dynamics among pions, Ct^n describes the interaction between pions and a nucleon, 
and the ehipsis stands for terms that involve pions and two or more nucleons. The individual Lagrangians 
are organized as follows: 

/:.. = £?2+£W + ... (2.27) 

and 

Cj,N = C^-^N + ^'fli + -^Sv + ■ ■ • : (2.28) 

where the superscript refers to the number of derivatives or pion mass insertions (chiral dimension) and the 
ellipsis stands for terms of higher dimensions. 

To construct chiral Lagrangians, we introduce the following SU{2) matrix U in flavor space which collects 
the Goldstone pion fields, tt: 

U = l+ r-TV- — - (r . 7r)3 + —^tv" + ..., (2.29) 

where /^r denotes the pion decay constant. In this expansion, the coefficient of the term linear in the pion 
field TT is fixed such as to produce the correct kinetic term in the tttt Lagrangian, below, and the coefficient 
of the quadratic term is chosen to satisfy the unitary condition U'^U = 1, at second order in the pion 
field. However, the coefficient a of the third order is arbitrary because of our freedom of choice for the 
interpolating pion fields (constrained only by U unitary and det U = 1). The coefficient of the fourth order 
is then dictated by the unitarity condition, U^U = 1, at fourth order. Note that (on-shell) observables must 
not depend on the choice for the pion fields or, in other words, they must not depend on the (unphysical) 
parameter a (and there are more such parameters as you continue to higher orders in the above expansion 
of U). Therefore, diagrams with vertices that involve three or four pions must always be grouped together 
such that the a dependence drops out. For more on this issue, see the calculation of the two-loop 27r and 
Stt contributions, below, where those vertices enter. Popular choices for the pion fields are the exponential 
parametrization U — exp(iT • tt//^) which corresponds to a = 1/6 and the so-called sigma representation 
U — {a + IT ■ Tv) //tt with a = -^//^ — tt-^ which is equivalent to a = 0. 
The leading order (LO) tttt Lagrangian is now given by [48] 

C'^^l ^ ^-fir [d^Ud^U^ +ml{U + U^)] , (2.30) 

where tr denotes the trace in flavor space and is the pion mass. Since Goldstone bosons can interact 
only when they carry momentum, the interaction between pions comes in powers of dfJJ . Only even powers 
are allowed because of Lorentz invariance. Note that the U field transforms under global chiral rotations via 

U^gLUg^ (2.31) 



with gn and g^ elements of SU(2)ji and SU{2)i^, respectively, cf. Eqs. (2.11) and (2.121. Since we consider 
global chiral rotations, gu and g^ do not depend on space-time and, therefore, dfj_U transforms in the same 



way as U. Thus, the first term in Eq. (2.30) is clearly chiral invariant. The second term breaks chiral 
symmetry explicitly with the coefficient chosen such as to reproduce the correct mass term, which can be 
seen by inserting U and expanding in numbers of pion fields: 

= ^-d,7T-d^n-^-mln' (2.32) 
+ ^^(^ ■ ^^^^)(^ ■ ^''^) - ^^'^M^ • + + ^(^') ' (2-33) 

^/tt Jit ^Jtt 

where we dropped the constant term f^rn^, since it does not contribute to the dynamics. 

Baryon fields can also be incorporated into the effective field theory in a chirally consistent manner. 
Gasser et al. [49] have derived the LO relativistic irN Lagrangian to be 

= * {n^'D, - Mn + Y^^15U^) * (2.34) 
11 



with 

D,, = d^ + r^ (2.35) 

the chirally covariant derivative which introduces the so-called chiral connection (an analogy to a gauge 
term) 

- l[C^,d,C]^\ied,^ + Cd,e) (2.36) 

= ^^-(^x^M^) + 0{7v')^ (2-37) 

representing a vector current that leads to a coupling of even numbers of pions with the nucleon. Besides 
this, the Lagrangian includes a coupling term which involves the axial vector 

= ^{i\^^C}=^{i^^,^-C^,,^^) (2.38) 

= -^T-d^7V+ ^J^iT ■ n){7V ■ d^TT) + ^7T^ (t • O^Tt) + ©(tT^) , (2.39) 

which couples an odd number of pions to the nucleon. The definition of ^ used in the above is 

+ 4'"-Ti^'^->-' + '^^^ (2.40) 



Note that = 1; therefore, a^(^^^) = {d^O^^ + ^d^C^ = 0, which is applied in Eqs. (2.36) and (2.38). 



Thus, more explicitly, the LO relativistic TriV Lagrangian, Eq. (2.34), reads 



Ci'l, = (^t^^d, - Mn - ^ rr ■ {n X d,7r) - |^7^5T ■ d,7v + . . }j ^ . (2.41) 

In this equation, is the relativistic four-component Dirac spinor field representing the nucleon, denotes 
the nucleon mass and 17^1 the axial-vector coupling constant. Numerical values will be given later. The term 
proportional to 5A/2/7r is the familiar axial- vector coupling of one pion to the nucleon, while the nonlinear 
term proportional to 1/4/^ is known as the Weinberg- Tomozawa coupling (a 2tt contact term) [571 155] . 
which is crucial for tt-N s-wave scattering at threshold. 

2.2.2. Heavy baryon formalism 

The relativistic treatment of baryons in chiral perturbation theory leads to problems. The reason for 
the problems is the fact that the time-derivative of a relativistic baryon field generates a factor E M 
(where M denotes the baryon mass) which is not small as compared to the chiral-symmetry breaking scale 

« 1 GeV; in fact, Mjv/A^ « 1. Note also that the nucleon mass does not vanish in the chiral limit. 
The consequence of all this is that the one-to-one correspondence between the expansion in pion loops, on 
the one hand, and the expansion in terms of small external momenta and pion masses, on the other side, is 
destroyed pi lTTSllTTO] . 

A solution to these problems has been proposed by Jenkins and Manohar [120) using effective field theory 
techniques originally developed by Georgi |121j for the study of heavy quark systems. The basic idea is to 
treat the baryons as heavy static sources ( "extreme non- relativistic limit" ) such that the momentum transfer 
between baryons by pion exchange is small as compared to the baryon mass. The expansion is performed 
in terms of these small momenta over the baryon mass and has become known as heavy baryon (HB) chiral 
perturbation theory (HBChPT). We will now briefiy sketch this approach. 

The four-momentum of the heavy baryon is parametrized as 

pt" = Mw^ -I- (2.42) 

where is the four-velocity satisfying = \ and is a small residual momentum, v ■ I <^ M. Defining 
projection operators 

p^^l±2^^ p+ + p- = l^ (2.43) 
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we introduce the so-called velocity-dependent fields 

N = e**^'"^P+^' , h = e^'^'^'T-^ , (2.44) 

such that the relativistic four-component Dirac spinor field representing the baryon can be written as 



4' 



-lAIv-x 



(2.45) 



The exponential factor in Eq. (2.441 eliminates the kinematical dependence on the baryon mass. To make 
the meaning of the fields N and h more transparent, consider a positive-energy plane wave solution of the 
Dirac equation 

ipp{x) ^ u{p,s)e~'P''' , u{p,s)^ E + M I 



2M 



(T-p 

E+M 



Xs 1 



(2.46) 



where I is the two-dimensional identity matrix, ^ E = \Jp'^ + M^, and Xs a Pauli spinor describing the 
spin state of the baryon. Assuming for the four-velocity = (1, 0, 0, 0), which is what we have in the rest 
frame of the baryon, the explicit expressions for the wave function components and are 







E + M f X. 
2M 

E + M 
2M 



-i{E-M)t+ip-x 



Xs 





-i{E-M)t+ip-x 



a-p 
E+mXs 



-i{E-M)t+ip-x 



(2.47) 

(2.48) 



where the represents a column vector that consists of two zeros. The arrow indicates the properly normal- 
ized spinor at leading order. So, for = (1,0,0,0), Np represents the large/upper component and hp the 
small/lower component of the Dirac wave function ij^p. Note also that the energy of Np and hp is different 
from 'il}p] namely, they carry the small residual energy [cf. Eq. (2.42)] 



Zn = £^ - M w 



P 

2M 



(2.49) 



Thus, Iq is zero in leading order and p^/2M at next-to-leading order (NLO). 

To derive the leading-order heavy-baryon Lagrangian, we start from the leading-order relativistic one. 



(2.50) 



and insert Eq. (2.45), assuming v'^ = (1,0,0,0) for simplicitjQ One immediate result is that the e 
factor of Eq. (2.45) generates a term 70M which kills the troublesome baryon mass term (— M) in the upper 
component projection yielding 



C 



(2.51) 



where the ellipsis stands for additional expressions involving the lower-component field h. Since N contains 
only upper components, Eq. (2.51) collapses to 



= ^ A 



a -u] N 



(2.52) 



plus 1 /M corrections which are generated by expressing h in terms of N via the equations of motion for N 
and h. The first such 1/AI corrections are shown in Eq. (2.56), below. 



*For a full-fledged derivation, see Ref. |113) 



13 



Equation (2.521 is the leading order TriV Lagrangian in the HB formahsm (as indicated by the hat). 
Expanding in numbers of pion fields, this Lagrangian reads 



(1) _ 



TXN 



1 

471 
,9.4(4a- 1) 



T • (tt X 5o7r) - T ■ {a ■ V)tv 



(t • tt) tt • ((7 • V)7r 



9Aa 2 



T • ((7 • V)7r 



N- 



where the ellipsis stands for terms involving four or more pions. 
At dimension two, the relativistic t:N Lagrangian reads 



(2.53) 



(2.54) 



The various operators O^"^^ are given in Ref. [122] . The fundamental rule by which this Lagrangian — as well 
as all the other ones — are assembled is that they must contain all terms consistent with chiral symmetry 
and Lorentz invariance (apart from other trivial symmetries) at a given chiral dimension (here: order two). 
The parameters Ci are known as low-energy constants (LECs) and are determined empirically from fits to 
ttN data. 

The HB projected ttN Lagrangian at order two is most conveniently broken up into two pieces, 



(2) 
■N 



(2) 

■N, fixed 



(2) 

■N, ct I 



(2.55) 



with 



and 



C 



(2) 

ttJV, ct 



^(2) 

■-ttN. fixed 



N 



N ■ 



D ■ D + ■ D,uo} 



AM 



2Mn 
2Mn AMnU 



N 



N 



WA 



• ( V doTv ~ doTv V 



V -(tt X Vtt) — {tt X Vtt)- V 



N 
N 



2 ci m 



C3 u^u'' 



2 2ci 22 f 9^ \ ^ 

Acim^ - y^m^TT + I C2 - I j^idoTT ■ don 



C4 



C3 



4M 



N 



a ■ {u X u) 



Jtt 



C4 ■ 



AM 



N 



N 



(2.56) 



N 



(2.57) 



where we neglected the isospin-breaking cs-term proportional to (m,j — Wd); the ellipsis represents terms 
involving more pions. 

Note that C'i^lf j^^^^ is created entirely from the HB expansion of the relativistic l'^^I^ and thus has no 

free parameters ("fixed"), while is made up by the new ttN contact terms proportional to the Ci 

parameters (plus those l/Mj^ corrections which happen to have the same mathematical structure as 
terms). 

At dimension three, the relativistic ttN Lagrangian can be formally written as 



(2.58) 
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with the operators, of'\ Usted in Refs. |122l 1123] : not ah 31 terms wih be of interest here. The new LECs 
that occur at this order are the di. Similar to the order two case, the HB projected Lagrangian at order 
three can be broken into two pieces. 



C 



(3) 
-N 



c 



(3) 



t:N , fixed 



c 



(3) 

-N, ct ' 



(2.59) 



with 41 fixed and il^^^^^ given in Refs. [THHlS]. 



^(3) 



2.3. Nucleoli contact Lagrangians 

Two-nucleon contact interactions consist of four nucleon fields (four nucleon legs) and no meson fields. 
Such terms are needed to renormalize loop integrals, to make results reasonably independent of regulators, 
and to parametrize the unresolved short-distance dynamics of the nuclear force. For more on the role of 
contact terms, see Section [431 

Because of parity, nucleon contact interactions come only in even powers of derivatives, thus. 



Cnn — C 



(0) 



C 



(2) 



c 



(4) 



-NN — '-NN ^ "-NN ^ '-NN 

The lowest order (or leading order) NN Lagrangian has no derivatives and reads [SOI [ST] 



r(0) 

'-NN 



-CsNNNN - -CriNaN) ■ (NaN) , 



(2.60) 



(2.61) 



where N is the heavy baryon nucleon field (N — A^^). Cs and Ct are unknown constants which are 
determined by a fit to the NN data. The second order NN Lagrangian can be stated as follows [51j 



■(2) 
NN 



C[ {N^Nf + {'VNNf -~C2{N^N)-{'VNN)-C3NN 



NV^N + W^NN 



NWN ■ (ViV X aN) + {WN)N ■ [Na x VA^) 



-iC'^NN{^N ■ G X V7V) - iC'^{NGN) ■ (VA^ x V7V) 
- iC!,6,kSji + C'^SuSkj + C^%40 [NakmNaidjN + dJVakNdjNaiN] 
' (C'l^SikSji + C'liSiiSkj + C'i2Sijdki ) Nukd.NdjNuiN 
1 



C[^[5,k5ji + SaSk,) + C[^5,,5ki [d.NcikdjN + d^Nakd.N] NaiN . 



(2.62) 



For a thorough discussion of second order contact Lagrangians, see Refs. |1241I125] . 

Similar to Cs and Ct of Eq. (2.61 ), the C[ of Eq. (2.62) are unknown constants which are fixed in a fit 
to the NN data. Obviously, these contact Lagrangians blow up quite a bit with increasing order, which is 
why we do not give >C^)v explicitly here. The NN contact potentials that emerge from these Lagrangians 
are given in Section [4. 3[ 

Besides the above contact Lagrangians involving two nucleons, there exist contact interactions among 
three or more nucleons representing nuclear many-body forces (cf. last term of Eq. (2.66), below, and 
Section [5]). 

2.4. Summary: Effective Lagrangians organized by interaction index A 

To summarize, the effective Lagrangian needed to derive nuclear forces includes the following parts: 



N 



c 



NN ■ 



(2.63) 



where the ellipsis stands for terms that involve two nucleons plus pions and three or more nucleons with or 
without pions, relevant for nuclear many-body forces (cf. last two terms of Eq. (2.66), below, and Section[5]). 

In previous sections, we organized the Lagrangians by the number of derivatives or pion-mass insertions. 
This is the standard way, appropriate particularly for considerations of tt-tt and n-N scattering. As it turns 
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out (cf. Section 3.1), for interactions among nucleons, it is sometimes useful to also consider the so-called 

n 



index of the interaction, 



A = d+ - - 2. 
2 



(2.64) 



where d is the number of derivatives or pion-mass insertions and n the number of nucleon field operators 
(nucleon legs). We will now re- write the HB Lagrangian in terms of increasing values of the parameter A. 
The leading-order Lagrangian reads. 



1 



1 



£^=0 = -a.Tr-a^TT- "-2_2 
2 ^ 



1 - 4a 



ido- r-{a- V)7r - t • (tt x Sott) 



8a- 1 
N 



9 4 



+N 
1 



ffA(4a- 1) 



(t • tt) tt • ((T • V)7r 



5Aa 2 



T • ((7 • V)7r 



^C'sNNNN - -CriNaN) ■ (NaN) + 



and subleading Lagrangians are. 



A=l 



= N' 



2Mn 4M^/, 



a ■ (v doTT - doTv v) 



4cim^ 



V -(tt X Vtt) - (tt x Vtt)- V 

2c- / "2 



N 



+ -H(a„ir-a"«-) 



1 



C4 



4A/Ary 2/2 



1 



C3 

'^{NN)N It - {a- V)-k\ N - ^E{NN){NtN) ■ (NtN) + 
4/7r L J 2 



^A=2 
^A=4 



(3) 
ttN 



c 



(2) 
WW 



(4) 



(2.65) 



(2.66) 

(2.67) 
(2.68) 



where the ellipses represent terms that are irrelevant for the derivation of nuclear forces up to fourth order. 



3. Nuclear forces from EFT: Overview 

In the beginning of Section [2] we listed the steps we have to take for carrying out the EFT program 
of a derivation of nuclear forces. So far, we discussed steps one to three. What is left are steps four (low- 
momentum expansion) and five (Feynman diagrams). In this section, we will say more about the expansion 
we are using and give an overview of the Feynman diagrams that arise order by order. 

3.1. Chiral perturbation theory and power counting 

Effective Lagrangians have infinitely many terms, and an unlimited number of Feynman graphs can be 
calculated from them. Therefore, we need a scheme that makes the theory manageable and calculable. This 
scheme which tells us how to distinguish between large (important) and small (unimportant) contributions 
is chiral perturbation theory (ChPT). 

In ChPT, graphs are analyzed in terms of powers of small external momenta over the large scale: 
{Q / K^Y , where Q is generic for a momentum (nucleon three-momentum or pion four-momentum) or a 
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pion mass and ^ 1 GeV is the chiral symmetry breaking scale (hadronic scale, hard scale). Determining 
the power v has become known as power counting. 

For the moment, we will consider only so-called irreducible graphs; the problem of reducible or iterative 
diagrams and their relevance for the NN system will be discussed later (cf. Section 4.5). By definition, an 
irreducible graph is a diagram that cannot be separated into two by cutting only nucleon lines. Following 
the Feynman rules of covariant perturbation theory, a nucleon propagator is Q^^, a pion propagator Q^^, 
each derivative in any interaction is Q, and each four-momentum integration Q^. This is also known as 
naive dimensional analysis. Applying then some topological identities, Weinberg obtained for the power of 
an irreducible diagram involving A nucleons iSC 



vw A-2C + 2L 



A, 



(3.1) 



with 



A, 



n, 

— - 2. 
2 



(3.2) 



where C denotes the number of separately connected pieces and L the number of loops in the diagram; di is 
the number of derivatives or pion-mass insertions and rii the number of nucleon fields (nucleon legs) involved 
in vertex i; the sum runs over all vertices i contained in the diagram under consideration. Note that A^ > 
for all interactions allowed by chiral symmetry. Purely pionic interactions have at least two derivatives 
{di > 2,ni = 0); interactions of pions with a nucleon have at least one derivative (di > l,ni = 2); and 
nucleon-nucleon contact terms {ui = 4) have di > 0. This demonstrates how chiral symmetry guarantees a 
low-energy expansion. 



The Weinberg formula Eq. (3.1 ) works well for connected diagrams with A < 2 nucleons and any number 



of pions, but there are problems when applied in systems with A > 3. 

To illustrate the problem, consider one-pion exchange (L = 0, Aj = 0) between two nucleons in an A = 2 
system (C = 1) which yields vw = 0. Now, when the same interaction occurs in an A = 3 environment, 
then C — 2, since one nucleon is not interacting, and Eq. (3.1) produces — —3. The reason for this 



result is that, by widespread convention |115| . particle states are normalized to Dirac 5-functions, 



(3.3) 



which carry dimension Q~^. A nucleon line passing through a diagram without interaction is represented 
by such a three-momentum-conserving 5-function. 

The above-illustrated A-dependence of the power v^/ is undesirable. As indicated, the reason for the 
problem is the dimension (—3) introduced for each nucleon through the normalization Eq. (3.3). Therefore, 
to fix the problem, we add (+3A) to Eq. (3.1) and subtract 6 to have the A = 2 case unaltered. Thus, we 
introduce the new power v — vw + 3A — 6, which reads explicitly. 



-2 + 2A - 2C + 2i + 2^ A, 



(3.4) 



This definition of the power v will be applied throughout this report. It is also what is used in other 



works |107j . Another way to get the formula Eq. (3.4 1 is to define v as the difference between vy/ and the 
smallest possible power, Urnm- v = vy/ — i-',nin- The minimal power z^min is obtained from Eq. (3.1) for the 
case of no loops (L = 0), A^ = for all vertices, and the maximal number of separately connected pieces, 
which is C = A — 1 if at least two nucleons interact; thus, t'min = 6 — 3^. 

Notice that, even though A still appears in Eq. (3.4), the formula is essentially A independent: when 
a non-interacting nucleon is added, A and C go up by one, which cancels. From the above derivation, it 



should also be clear that Eq. (3.4) is suitable only for systems with A> 2 nucleons while Eq. (3.1 ) may be 



used for connected graphs with A < 2. 

An alternative method for developing a reasonable power formula has been presented by Friar |126j . 
who considers the expectation value of an m-nucleon operator (to < A) in A-nucleon space. Working in 
configuration space generates additional phase-space factors of power 3(^ — 1) to be added to Eq. (3.1). 
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Thus, Friar's power is vp = i^w + — 3 = v + 3. The additional power of three by which Friar differs from 



Eq. (3.4) is due to an additional momentum-space integration that converts the units of a momentum-space 
potential into units of energy, which is reasonable. Finally, we note that normalizing the nucleon states in 
a box to a Kronccker-(5 (instead of using the Dirac (5-function normalization of the continuum) should lead 



to Eq. (3.4) in a straightforward fashion. 

In any case, the most important observation from power counting is that the powers are bounded from 
below and, specifically, v > 0. This fact is crucial for the convergence of the low-momentum expansion. 



Moreover, the power formula Eq. (3.4) allows to predict the leading orders of connected multi-nucleon 



forces. Consider a m-nucleon irreducibly connected diagram (m-nucleon force) in an A-nucleon system 



(to < A). The number of separately connected pieces is C = A — to + 1. Inserting this into Eq. (3.4| 
together with L = and = yields ly = 2m — 4. Thus, two-nucleon forces (to = 2) start a.t ly — 0, 

three-nucleon forces (to = 3) at = 2 (but they happen to cancel at that order), and four- nucleon forces at 
v — 4 (they don't cancel). More about this in the next sub-section and Section [5] 

For later purposes, we note that for an irreducible NN diagram (A ^ 2, C ~ 1), the power formula 
collapses to the very simple expression 

ly = 2L + J2^^■ (3-5) 

i 

In summary, the chief point of the ChPT expansion is that, at a given order v, there exists only a finite 
number of graphs. This is what makes the theory calculable. The expression (Q / A^Y'^^ provides a rough 
estimate of the relative size of the contributions left out and, thus, of the accuracy at order ly. In this sense, 
the theory can be calculated to any desired accuracy and has predictive power. 

3.2. The hierarchy of nuclear forces 

Chiral perturbation theory and power counting imply that nuclear forces emerge as a hierarchy controlled 
by the power Fig. [T] 

In lowest order, better known as leading order (LO, ly — 0), the NN amplitude is made up by two 
momentum-independent contact terms (^ Q"), represented by the four-nucleon-leg graph with a small-dot 
vertex shown in the first row of Fig. [ij and static one-pion exchange (IPE), second diagram in the first 
row of the figure. This is, of course, a rather crude approximation to the two-nucleon force (2NF), but 
accounts already for some important features. The IPE provides the tensor force, necessary to describe the 
deuteron, and it explains NN scattering in peripheral partial waves of very high orbital angular momentum. 
At this order, the two contacts which contribute only in ^-waves provide the short- and intermediate-range 
interaction which is somewhat crude. 

In the next order, ^ = 1, all contributions vanish due to parity and time-reversal invariance. 

Therefore, the next-to-leading order (NLO) is = 2. Two-pion exchange (2PE) occurs for the first 
time ("leading 2PE") and, thus, the creation of a more sophisticated description of the intermediate-range 



interaction is starting here. Since the loop involved in each pion-diagram implies already ly — 2 [cf. Eq. (3.5 1], 
the vertices must have = 0. Therefore, at this order, only the lowest order ttNN and ttttNN vertices 
are allowed which is why the leading 2PE is rather weak. Furthermore, there are seven contact terms of 
0{Q^), shown by the four-nucleon-leg graph with a solid square, which contribute in S and P waves. The 
operator structure of these contacts include a spin-orbit term besides central, spin-spin, and tensor terms. 
Thus, essentially all spin-isospin structures necessary to describe the two-nucleon force phenomenologically 
have been generated at this order. The main deficiency at this stage of development is an insufficient 
intermediate-range attraction. 

This problem is finally fixed at order three (ly = 3), next-to-next-to-leading order (NNLO). The 2PE 
involves now the two-derivative ttttNN seagull vertices (proportional to the LECs) denoted by a large solid 
dot in Fig. [l] These vertices represent correlated 2PE as well as intermediate A(1232)-isobar contributions. 
It is well known from the meson phenomenology of nuclear forces [181 120! that these two contributions are 
crucial for a realistic and quantitative 2PE model. Consequently, the 2PE now assumes a realistic size 
and describes the intermediate-range attraction of the nuclear force about right. Moreover, first relativistic 
corrections come into play at this order. There are no new contacts. 
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2N Force 



3N Force 4N Force 




NLO 



itx m\ 

+... +... 

Figure 1: Hierarchy of nuclear forces in ChPT. Solid lines represent nucleons and dashed lines pions. Small dots, large solid 
dots, solid squares, and solid diamonds denote vertices of index A = 0, 1, 2, and 4, respectively. Further explanations are 
given in the text. 

The reason why we talk of a hierarchy of nuclear forces is that two- and many-nucleon forces are created 
on an equal footing and emerge in increasing number as we go to higher and higher orders. At NNLO, the 
first set of nonvanishing three-nucleon forces (3NF) occur HOIITI], cf. column '3N Force' of Fig.[lj In fact, at 
the previous order, NLO, irreducible 3N graphs appear already, however, it has been shown by Weinberg I52j 
and others [TDl 11271 1128] that these diagrams all cancel. Since nonvanishing 3NF contributions happen first 
at order (Q/A^)^, they are very weak as compared to 2NF which start at {Q/ K^)^ . 

More 2PE is produced aX v ^ A, next-to-next-to-ncxt-to- leading order (N'^LO), of which we show only 
a few symbolic diagrams in Fig. [T] Two-loop 2PE graphs show up for the first time and so does three-pion 
exchange (3PE) which necessarily involves two loops. 3PE was found to be negligible at this order [57l [58] . 
Most importantly, 15 new contact terms ^ arise and are represented by the four-nucleon-leg graph with 
a solid diamond. They include a quadratic spin-orbit term and contribute up to Z3-waves. Mainly due to 
the increased number of contact terms, a quantitative description of the two-nucleon interaction up to about 
300 MeV lab. energy is possible, at N'^LO (for details, see below). Besides further 3NF, four-nucleon forces 
(4NF) start at this order. Since the leading 4NF come into existence one order higher than the leading 3NF, 
4NF are weaker than 3NF. Thus, ChPT provides a straightforward explanation for the empirically known 
fact that 2NF > 3NF > 4NF .... 

4. Two-nucleon interactions 

The last section was just an overview. In this section, we will fill in all the details involved in the ChPT 
development of the NN interaction; and 3NF and 4NF will be discussed in Section [5j We start by talking 





NNLO 



N^LO 
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(a) 




(b) (c) (d) (e) 




(f) (g) (h) (i) 

Figure 2: One-pion exchange contributions. Diagram (a) represents the leading-order, while graphs (b)-(i) contribute to 
renormalization. Notation as in Fig. [l] Dia grams that result from interchange of nucleon lines and/or time reversal are not 
shown. 



about the various pion-exchange contributions. 

Pion-exchange contributions in ChPT 

Based upon the effective pion Lagrangians of Section 2.2 we will now derive the pion-exchange contri- 
butions to the NN interaction order by order. 

As noted before, there are infinitely many pion-exchange contributions to the NN interaction and, thus, 
we need to get organized. First, we arrange the various pion-exchange contributions according to the number 
of pions being exchanged between the two nucleons: 

K - V,^ + V2.. + V3^ + ... , (4.1) 

where the meaning of the subscripts is obvious and the ellipsis represents An and higher pion exchanges. 
Second, for each of the above terms, we assume a low-momentum expansion: 

V,. - Vl'j+V!'j+Vi!^+V}t' + ... (4.2) 
= V^^^ + V^!^+V,^^'> + ... (4.3) 
Vs. = V^::^ + ..., (4.4) 

where the superscript denotes the order v and the ellipses stand for contributions of fifth and higher orders. 
Due to parity and time-reversal, there are no first order contributions. Moreover, since n pions create 



L = n — I loops, the leading order for n-pion exchange occurs at = 2n — 2 [cf. Eq. (3.5)] 



In the following subsections, we will discuss Vitt, V27r, and Vstt, one by one and order by order. 
4.I.I. One-pion exchange 

One-pion exchange (IPE) diagrams are shown in Fig. [2j At leading order (LO, = 0), we have the 
well-known static IPE, diagram (a) of Fig. [2] which is given by [for notation, see Eq. (4.8) below]: 



T. ^ 9a ai-q a2-q , , 
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On-shell and in the center-of-mass system (CMS), there are no relativistic corrections at any order. Off- 
sheh corrections come into play when the relativistic IPE is iterated, i.e., in four-dimensional planar box 
diagrams. These corrections will be taken into account in the evaluation of those box diagrams up to the 
given order (see below). 

At second order (NLO), the IPE gets renormalized due to one-loop graphs and counter term insertions, 
shown in the second and third row of Fig. [I Graphs (b) and (c) renormalize the nucleon and graphs (d) 
and (e) the pion lines. Diagrams (f)-(i) renormalize the pion- nucleo n coupling. In the one-loop graphs, 
all vertices are from the leading order Lagrangian C^^^ , Eq. (2.65), while counter term insertions stem 



from C^~^ . Note that graph (f) vanishes because it involves an odd power of the loop momentum that is 
integrated over. In graph (i), the solid square includes the dis-vertex from [which is part of L^^^ , 



Eq. (2.67)]. This correction, which is known as the Goldberger-Treiman discrepancy, can be taken care of 
by replacing 

9A — > 9A~ 2di8m^ . (4.6) 

For more details on the NLO corrections to the IPE see Ref. [129] . At NNLO, there are further one-loop 
corrections, but with one subleading vertex, which renormalize qa- 

Finally, at N^LO, there are two-loop corrections with leading vertices only, as well as one-loop and tree 
diagrams including (sub-)subleading vertices. It has been shown [55, 57 that these contributions renormalize 
various LECs and the pion mass rriTr, but do not generate any ttN form-factor-like functions. Furthermore, 
a correction to the Goldberger-Treiman discrepancy arises at N^LO. 

We use gA — 1.290 (instead of gA ~ 1.276 [1301 ) to account for the Goldberger-Treiman discrepancy. 
Via the Goldberger-Treiman relation, g^^NN — 9aMn / Jt^, our value for gA together with /^r = 92.4 MeV 
and Mn — 938.918 MeV implies g'^j^j^/Air — 13.67 which is consistent with the empirical value g'^j^jq/Air — 
13.65 ± 0.08 obtained from nN and NN data analysis [1311 fl32 . The renormalizations of m^r, and Mm 
are taken care of by working with their physical vales. 



In summary, the familiar expression for IPE, Eq. (4.5), is appropriate to at least fourth order 



4-1-2. Two-pion exchange 

The exchange of two or more pions always involves loop diagrams, which implies that we are faced with 
a non-trivial problem. To conduct such calculations, a scheme of field-theoretic perturbation theory as well 
as regularization and renormalization methods for dealing with divergent loop integrals must be adopted. 
Ordonez, Ray, and van Kolck [S31[53], who were the first to calculate chiral 27r-exchange to NNLO, used 
time-ordered perturbation theory and a Gaussian cutoff function for regularization. The Juelich group [107j 
developed a method involving unitary transformations (to get rid of the energy-dependence of time-ordered 
perturbation theory) and applied so-called spectral function regularization. Both, the Brazil and the 

Munich [551 1601 161| groups use covariant perturbation theory and dimensional regularization, but neverthe- 
less, their works differ substantially in detail. We will follow here the method chosen by the Munich group 
since we believe it to be the most efficient and elegant one. In this approach, one starts with the relativistic 



versions of the ttN Lagrangians (cf. Section 2.2) and sets up four-dimensional (covariant) loop integrals. 
Relativistic vertices and nucleon propagators are then expanded in powers of l/M^- The divergences that 
occur in conjunction with the four-dimensional loop integrals are treated by means of dimensional regular- 
ization, a prescription which is consistent with chiral symmetry and power counting. The results derived 
in this way are the same obtained when starting right away with the HB versions of the -kN Lagrangians. 
However, as it turns out, the method used by the Munich group is more efficient in dealing with the rather 
tedious calculations and particularly useful in conjunction with the planar box diagram. To give the reader 
a taste of the rather involved calculations, we present in [Appendix Bj the explicit evaluation of the NLO 
diagrams shown in Fig. [3] 

The results will be stated in terms of contributions to the momentum-space NN amplitude in the CMS, 
which takes the general form 

V{p',p) = Vc +T1-T2WC 

+ [Vs +'ri- T2WS ] Si - ^2 
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(Q/A,)^ 



Figure 3: Leading two-pion exchange contributions to the A'^A'^ interaction. Notation as in Fig.^ 



+ [Vls + Ti- T2 Wls] • [q X k) ) 

+ [Vt + Ti ■ T2 Wt ] Si ■ q $2 ■ q 

+ [VaL+Ti- T2 W„L ] l?l • (g X fc ) (72 • (9 X fc ) , 

where p ' and p denote the final and initial nucleon momenta in the CMS, respectively; moreover, 



k 
S 



P -P 



- 0-2 ) 



is the momentum transfer, 
the average momentum, 
the total spin. 



(4.7) 



(4.8) 



and 2 and Ti_2 are the spin and isospin operators, respectively, of nucleon f and 2. For on-energy-shell 
scattering, Va and Wa {a = LS,T,aL) can be expressed as functions of q and k (with g = |(7| and 
k=\k\), only 

Our notation and conventions are similar to the ones used by the Munich group (SSHHOlini] except for two 
differences: our spin-orbit potentials, Vls and Wls, differ by a factor of (+2) and all other potentials differ 
by a factor of (—1) from the Munich amplitudes. Our definitions arc more in tune with what is commonly 
used in nuclear physics. 

In all expressions given below, we will state only the nonpolynomial contributions to the NN amplitude. 
Note, however, that dimensional regularization typically generates also polynomial terms which are, in part, 



infinite or scale dependent (cf. Appendix B|. These polynomials are absorbed by the contact interactions 
to be discussed in a later section. 

Next-to-leading order (NLO). The leading two-pion exchange appears at second order {v — 2, next-to- 
leading order, NLO) and is shown in Fig. [s] Since a loop creates already v = 2 [cf. Eq. (3.5)], the vertices 
involved at this order have index = 0, i.e., they are from the leading order Lagrangian C^^'^ , Eq. (2.65), 
where the ttN vertices carry only one derivative. These vertices are denoted by small dots in the figures. 
The rather complicated evaluation of these diagrams is presented in [Appendix B[ 

Concerning the planar box diagram in Fig. [Sj we should note that we include only the non-iterative part 
of this diagram which is obtained by subtracting the iterated IPE contribution Eq. (4.25) or (4.26), below, 
but using Ml/Ep « Mf^/Epn « Mn at this order (NLO). Summarizing all contributions from irreducible 
two-pion exchange at second order, one obtains: 



Wn = - 



Vt = 



384712 

f 



where 



and 



Aml{hg\ - Ag\ - 1) + q\22,g\ - \^g\ - 1) 
64^2/^ ' 



q ImT^ 



w . 



\/4to2 



(4.9) 
(4.10) 

(4.11) 
(4.12) 



As will be demonstrated in Sec. |4.2[ below; this part of the 2PE is rather weak and insufficient to properly 
describe the NN interaction at intermediate range. 
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Figure 4: Two-pion exchange contributions to the A'^A'^ interaction at order t hree in small momenta (NNLO). Basic notation as 
in Fig.[l] Large solid dots denote vertices from the Lagrangian C^~^ , Eq. (2.66i, proportional to the LECs c^. Symbols with 
an open circles are relativistic 1/Mjv corrections which are also part of C^~^ . Only a few representative examples of 1/Mm 
corrections are shown. Note that all football diagrams shown in this figure vanish. 



Next-to-next-to-leading order (NNLO). The two-pion exchange diagrams of order three {v = 3, next-to- 
next-to-leading order, NN LO) a re very similar to the ones of order two, except that they contain one 
insertion from C^^^ ^ Eq. (2.66). The resulting contributions are typically either proportional to one of 
the low-energy constants Ci or they contain a factor 1/Mm- Notice that relativistic 1/Mm corrections 
derive from vertices and nucleon propagators. In Fig. |4j we show in row one the diagrams with one vertex 
proportional to (large solid dot), and in row two and three a few representative graphs with a l/Mjy 
correction (symbols with an open circle). The number of l/Af^r correction graphs is large and not all 
are shown. Note that all football diagrams vanish at this order, because the loop integrals involve odd 
powers of the time-component of the loop momentum. Again, the planar box diagram is corrected for a 
contribution from the iterated IPE. If the iterative 2PE of Eq. (4.25), below, is used, the expansion of the 
factor M]q/Ep = Mm — /2Mm -f . . . is applied and the term proportional to {—p'^/2Mm) is subtracted 
from the third order box diagram contribution. Then, one obtains for the full third order |55| : 



Vc 
Wc 

Vt 
Wt 
Vls 



^9a [ dlml 
167r/4 \ WMnw^ 



2m^(2ci-C3)-gMc3-f 



35i 



W^A{q) 
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1287rM^/^ 
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[Ami + 



5^(4m2 +3g2)] w^A{q)] 
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w'A{q) 



with 



and 



Ml) = ^ arctan 

Zq 2m7r 



(4.13) 
(4.14) 
(4.15) 
(4.16) 
(4.17) 
(4.18) 

(4.19) 
(4.20) 



This contribution to the 2PE is the crucial one, because it provides an intermediate-range attraction of 



proper strength (Sec. 4.2). The iso-scalar central potential, Vc, is strong and attractive due to the LEG C3, 
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Figure 5: One-loop 27r-exchange contributions to the A'^ interaction at order four. Notation as in Fig. |4] Moreover, symbols 
with a large solid dot and an open circle denote 1/Mjv corrections of vertices proportional to C j. Sym bols v^ith two open circles 
mark relativistic corrections. Both corrections are part of the Lagrangian C^~^, Eq. (2.67 1. Representative examples 

for all types of one-loop graphs that occur at this order are shown. 



which is negative and of large magnitude (cf. Table [2] below). Via resonance saturation, C3 is associated with 
TT-TT correlations ('cr meson') and virtual A-isobar excitations, which create the most crucial contributions 
to 2PE in the frame work of conventional meson theory |11[|20| . The configuration-space expressions, which 
correspond to the above momentum-space potentials, are given in Ref. [55j. where also a detailed comparison 
with meson-exchange potentials is conducted. 

If the iterative 2PE defined in Eq. (4.261, below, is applied, the l/M^ terms are slightly different. As 



derived in Appendix C the changes are taken care of by adding to Eqs. (4.13 )-(4.16 1 the following terms 



Vc 
Wc 

Vt 
Wt 



;Vs 



i9\ 



5127r/4M, 



N 



-(m^ -I- w'^A{q)) 



-,Ws = - 



2567r/4Mjv 



(4.21) 
(4.22) 
(4.23) 
(4.24) 



Next-to-next-to-next-to-leading order (N^LO). This order {ly — 4) is very involved. The contributions can 
be subdivided into two groups, one-loop graphs. Fig. [sj and two-loop diagrams. Fig. [6j Applying Eq. (3.5), 
it is easy to verify that all contributions result in ^' = 4. We have relegated the comprehensive mathematical 
expressions of this order to [Appendix D The net effect of this contribution is small indicating a trend 
towards convergence (cf. Sec. 4.2). 
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Figure 6: Two-loop 27r-exchange contributions at order four. Basic notation as in Fig.[T] The oval stands for all one-loop ixN 
graphs some of which are shown in the lower part of the figure. The solid square represe nts vertices pr oportional to the LECs 
di introduced in /^'^i which is part of 



Eq. I 2.67 1. More explanations are given in 



Appendix D 



Iterated one-pion- exchange. Besides all the irreducible 2PE contributions presented above, there is also the 
reducible 2PE which is generated from iterated IPE. This "iterative 2PE" is the only 2PE contribution 
which produces an imaginary part. Thus, one wishes to formulate this contribution such that relativistic 
elastic unitarity is satisfied. There are several ways to achieve this. 
Kaiser et al. [55] define the iterative 2PE contribution as follows, 



(KBW) 



2ir,it 



Ml 

En 



(4.25) 



with ViTT given in Eq. (4.5) 



As it will turn out (cf. Section 4.4), in the development of a chiral NN potential, it is useful to adopt 
the relativistic scheme by Blankenbecler and Sugar [T33] (BbS). In this approach, the iterated IPE has the 
following form: 

,(EM),^, ^_ f d^p" Ml Vi^{p',p")Vi^{p",p) 



(2^)3 Ep 



p^ — p" 



(4.26) 



Equations (4.25) and (4.26) state the iterative 2PE to all orders. On the other hand, the covariant box 



diagram is calculated order by order and, therefore, requires a subtraction of the iterative 2PE order by 



order. For this, the expansi on M ^/Ep = Mjv — p /2Mm + ... is applied in Eq. (4.251 and M^/E^ 



Mn —p" I^Mn + ... in Eq. (4.261. At NLO, both choices for the iterative 2PE collapse to the same, while 



at NNLO there are differences which give rise to the correction terms Eqs. (4.21 )-( 4.24) when Eq. (4.261 is 



used (see Appendix C). 



4- 1.3. Three-pion exchange 

Since the exchange of three pions involves two loops, three-pion exchange (3PE) starts at order four 
(N^LO). Two loops generat e v = 4 and, therefore, at this order, all vertices have to be from the leading 
order Lagrangian £.^^^ , Eq. (2.65 1. One can distinguish between three groups of diagrams, namely, diagrams 
proportional to g\ (first and second row of Fig. [?]) , g\ (third row of the figure) , and g\ (fourth row) . The 
graphs in the first row involve the 3ttNN contact vertex from the leading ttN Lagrangian and the An vertex 
from the leading tttt Lagrangian (both belong to C^^^). Note that these vertices contain the unphysical 
parameter a which was introduced to parametrize the interpolating pion fields, Eq. (2.29). Therefore, this 
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Figure 7: Leading three-pion exchange contributions to the NN interaction. Notation as in Fig. ^ Diagrams that result from 
interchange of nucleon lines and/or time reversal are not shown. 



group of graphs has to be calculated together such that the parameter a drops out, as it should, since 
measurable quantities must not depend on a. 

The 3PE contributions at N'^LO have been calculated by the Munich group and found to be negligible [571 
I58j . This is not surprising since only the leading vertices are involved which are known to be weak. As 
discussed, for similar reasons, the leading 2PE also turned out to be rather small (even though not negligible). 
Note, however, that the subleading 3PE contributions which involve one vertex from the Lagrangian 
Eq. ( |2.66 1, proportional the LEG c;, have been found to be sizable [59]. They contribute at fifth order 



(N^'^LO). Our current analysis is restricted to N^LO where only the leading 3PE contributes which wc will 
ignore because of its negligible strength. 

4-2. Perturbative NN scattering in peripheral partial waves 

Nucleon-nucleon scattering in peripheral partial waves is of special interest — for several reasons. First, 
these partial waves probe the long- and intermediate-range of the nuclear force. Due to the centrifugal 
barrier, there is only small sensitivity to short-range contributions and, in fact, the N'^LO contact terms 



make no contributions for orbital angular momenta L > 3 (cf. Section 4.3 ). Thus, for F and higher waves and 
energies below the pion-production threshold, we have a window in which the NN interaction is governed by 
chiral symmetry alone (chiral one- and two-pion exchanges) , and we can conduct a relatively clean test of how 
well the theory works. Using values for the LECs from ttN analysis, the NN predictions are even parameter 
free. Moreover, the smallness of the phase shifts in peripheral partial waves suggests that the calculation 
can be done perturbatively. This avoids the complications and the possible model-dependence that the 



non-perturbative treatment of the Schroedinger equation, necessary for low partial waves (Section 4.5), is 
beset with. Because of the importance of peripheral NN scattering, many calculations based upon chiral 
dynamics can be found in the literature [S51 UHl El IHIl II31HI3S] ■ We will follow Ref. [SJ. 

Since we will compare the predictions with neutron-proton {np) phase shifts, we will specifically calculate 
np scattering in this sub-section. Defining, 

yi^(m^) = - — ^, (4.27) 
4/^ + rn^ 

the correct IPE for np scattering is given by 

V[7\v\v) = -14.(m^o) + (-1)^+1 2yi^(m^±) , (4.28) 
where / denotes the isospin of the two- nucleon system. We use the pion masses given in Table [2] and 

Mat = ^j^^p^" ^ 938.9182 MeV . (4.29) 
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Table 2; Hadron masses and low-energy constants (LECs). m^± and m^o denote the charged- and neutral-pion masses, 
and Mp and M„ are the proton and neutron masses, respectively, in units of MeV. qa is the axial-vector coupling constant 
(dimensionless) a nd the pion decay constant (in units of MeV). The Ci belong to the dimension-two ttN Lagrangian, 
Eqs. (|2.57[ l and l|2.66| l, and are in units of GeV ^, whereas the di are associated with the dimension-three Lagrangian, 
Eqs. | |238f and ^2m^ and are in units of GeV . Column "NN Potential" lists the values used for a NN potential at N'^LO 
presented in Section |4.6.3| while column "Peripheral perturbative NN" shows the parameters applied in the peripheral NN 
scattering calculations of Section |4.2| Finally, the last column displays values from empirical determinations (see text for 
comments). 



NN Potential 



(Sec. 4.6.3) 



Peripheral 
perturba tive NN 
(Sec. Ill 



Empirical 



Mp 
M„ 
QA 

Cl 
C2 
C3 

_ C4 _ 

di + d2 

ds 

_ d5 _ 

di4 — di5 



139.5702 
134.9766 
938.2720 
939.5653 
1.29" 
92.4 
-0.81 
2.80 
-3.20 
5.40 
3.06 
-3.27 
0.45 
-5.65 



139.5702 
134.9766 
938.2720 
939.5653 
1.29" 
92.4 
-0.81 
3.28 
-3.40 
3.40 
3.06 
-3.27 
0.45 
-5.65 



139.57018(35) [HS] 
134.9766(6) ^TE\ 
938.272013(23) [TTS] 
939.565346(23) [TTi] 
1.2759(45) [120] 
92.2 ±0.2 TH] 
-0.81 ± 0.15'' 
3.28 ±0.23^= 
-4.69 ± 1.34'' 
3.40 ±0.04'' 
3.06 ±0.21^= 
-3.27 ±0.73^= 
0.45 ±0.42^= 
-5.65 ±0.41= 



"Confer discussion of Goldberger-Treiman discrepancy in Section [4.1.11 
''Table 1, Fit 1 of Ref. [T3S]. 
^Table 2, Fit 1 of Ref. 



Also in the iterative 2PE, we apply the correct np IPE, i.e., in Eq. (4.251 we replace Vi^ with Vi^^\ Thus, 
the perturbative relativistic T-matrix for np scattering, taking the exchange of up to two pious into account, 
is calculated in the following way, 



r(p',p) = v,^7\p\p) + vt\p\p} 



(KBW.np) 



{p',p) + vUp',p). 



(4.30) 



As discussed, the expression for Vi^^ is good to any order we will consider in this article (cf. Section 4.1.1 1 
and V2^f^'"'^'' includes all orders. There is no need to break the latter term up into orders, because the 



admixture of (very small) higher order contributions from 



KBW,np) 



will not affect the accuracy we are 



working at. The most important term in the above equation is V2t^, the irreducible 2PE contributions, for 
which we have the low-momentum expansion. 



'(2) 



'(3) 



'(4) 



(4.31) 



with the various V2^''^ given in Section 



4.1.2 



In the calculation of Vj^^., we use the average pion mass 



rriTT = 138.039 MeV and, thus, neglect the charge-dependence due to pion-mass splitting in irreducible 
diagrams. The charge-dependence that emerges from irreducible 27r exchange was investigated in Ref. |137) 
and found to be negligible for partial waves with L > 3. 



For the T-matrix given in Eq. ( 4.30 1 , we calculate phase shifts for partial waves with i > 3 and Tiab < 300 
MeV (see Ref. |:67! for the details of this calculation). The LECs used in this calculation are shown in 
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Figure 8: _F-wave phase shifts of neutron-proton scattering for laboratory kinetic energies below 300 MeV. We show the predic- 
tions from chiral pion exchange to leading order (LO), next-to-leading order (NLO), next-to-next-to-leading order (NNLO), and 
next-to-next-to-next-to-leading order (N3LO). Note that in the NNLO and N3LO curves cannot be distinguished on the 
scale of the figure. The solid dots and open circles are the results from the Nijmegen multi-energy np phase shift analysis [139] 
and the VPI/GWU single-energy np analysis SM99 I140| . respectively. 



Table [2] column "Peripheral perturbative NN" . Note that many determinations of the LECs, Ci and di, can 
be found in the literature. The most reliable way to determine the LECs from empirical ttN information 
is to extract them from the ttN amplitude inside the Mandelstam triangle (unphysical region) which can 
be constructed with the help of dispersion relations from empirical t:N data. This method was used by 
Biittiker and Meifiner |138) . Unfortunately, the values for C2 and all di parameters obtained in Ref. |138| 
carry uncertainties, so large that the values cannot provide any guidance. Therefore, in Table [2] only ci, 
C3, and C4 are from Ref. |138j . while the other LECs are taken from Ref. |123j where the TriV amplitude in 
the physical region was considered. To establish a link between ttN and NN, we apply the values from the 
above determinations in our calculations of the NN peripheral phase shifts. In general, we use the central 
values; the only exception is C3, where we choose a value that is, in terms of magnitude, about one standard 
deviation below the one from Ref. |138j . With the exception of C3, phase shift predictions do not depend 
sensitively on variations of the LECs within the quoted uncertainties. 

In Figs. [8] and [9j we show the phase-shift predictions for neutron-proton scattering in F and G waves, 
respectively, for laboratory kinetic energies below 300 MeV. The orders displayed are defined as follows: 



Leading order (LO) is just IPE, first term on the r.h.s. of Eq. (4.301. 

Next-to-leading order ( NLO) i nclud es th e first two terms on the r.h.s. of Eq. (4.301 (IPE & iterative 

Eqs. (4.91 and (4.10)]. 



2PE) plus ^ [Section 



4.1.2 



Next-to-next-to-leading order (NNLO) consists of NLO plus 1^2^ [Section 
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Figure 9: Same as Fig. [8] but for G waves. 



It is clearly seen in Figs. |8] and [9] that the leading 2n exchange (NLO) is, in general, rather small, insufficient 
to explain the empirical facts in most partial waves. In contrast, the next order (NNLO) is very large; in 
some cases, several times NLO. This is due to the ttttNN contac t inte ractions proportional to the LECs 
Ci that are introduced in the sub-leading Lagrangian £^=^, Eq. (2.661. These contacts are supposed to 
simulate the contributions from intermediate A-isobars and correlated 27r exchange which are known to be 
large and crucial for a realistic model for the NN interaction at intermediate ranges (see, e. g., Ref. jlO]). 
At N^LO a clearly identifiable trend towards convergence emerges. In G waves (except for see 



Appendix D.3 for a discussion of this issue), N'^LO differs very little from NNLO implying that we have 



reached convergence. Also ^F^ and ^i^4 appear fully converged. However, in ^^^2 and ^F^, N^LO differs 
noticeably from NNLO, but the difference is much smaller than the one between NNLO and NLO. This is 
what we perceive as a trend towards convergence. Individual N'^LO contributions to peripheral phase shifts 
are shown in Appendix D.3| 



In Figs. [10] and 11 we conduct a comparison between the predictions from chiral one- and two-pion 
exchange at N'^LO and the corresponding predictions from conventional meson theory (curve 'Bonn'). As 
representative for conventional meson theory, we choose the Bonn meson-exchange model for the NN inter- 
action |20| . since it contains a comprehensive and thoughtfully constructed model for 27r exchange. This 27r 
model includes box and crossed box diagrams with NN, NA, and AA intermediate states as well as direct 
TTTT interaction in S- and P-waves (of the tttt system) consistent with empirical information from ttN and tttt 
scattering. Besides this the Bonn model also includes (repulsive) w-meson exchange and irreducible diagrams 
of TT and p exchange (which are also repulsive). However, note that in the phase shift predictions displayed 
in Figs. [TO|and[TT] the "Bonn" curve includes only the In and 2tt contributions from the Bonn model; the 
short-range contributions are left out since the purpose of the figure is to compare different models/theories 
for 7r + 27r. In all waves shown (with the usual exception of^G^) we see, in general, good agreement between 
N^LO and Bonn. In and above 150 MeV and in above 250 MeV the chiral model to N^LO 
is more attractive than the Bonn 27r model. Note, however, that the Bonn model is relativistic and, thus. 



29 





100 200 300 100 200 300 

Lab. Energy (MeV) Lab. Energy (MeV) 

Figure 10: F-wave phase shifts of neutron-proton scattering for laboratory kinetic energies below 300 MeV. We show the results 
from one-pion-exchange (IPE), and one- plus two-pion exchange as predicted by ChPT at next-to-next-to-next-to-leading order 
(N3LO) and by the Bonn Full Model 20j (Bonn). Note that the "Bonn" curve does not include the repulsive uj and irp exchanges 
of the full model, since this figure serves the purpose to compare just predictions by different models/theories for the tt + 2it 
contribution to the A' interaction. Empirical phase shifts (solid dots and open circles) as in Fig. Is] 



includes relativistic corrections up to infinite orders. Thus, one may speculate that higher orders in ChPT 
may create some repulsion, moving the Bonn and the chiral predictions even closer together |141j . 

The 27r exchange contribution to the NN interaction can also be derived from empirical ttN and tttt 
input using dispersion theory, which is based upon unitarity, causality (analyticity), and crossing symmetry. 
The amplitude NN — > tttt is constructed from ttN — ^ ttN and ttN ttttN data using crossing properties 
and analytic continuation; this amplitude is then 'squared' to yield the NN amplitude which is related to 
NN by crossing symmetry |T6j . The Paris group [TTl [18] pursued this path and calculated NN phase shifts 
in peripheral partial waves. Naively, the dispersion-theoretic approach is the ideal one, since it is based 
exclusively on empirical information. Unfortunately, in practice, quite a few uncertainties enter. First, 
there are ambiguities in the analytic continuation and, second, the dispersion integrals have to be cut off 
at a certain momentum to ensure reasonable results. In Ref. |20j . a thorough comparison was conducted 
between the predictions by the Bonn model and the Paris approach and it was demonstrated that the 
Bonn predictions always lie comfortably within the range of uncertainty of the dispersion-theoretic results. 
Therefore, there is no need to perform a separate comparison of our chiral N'^LO predictions with dispersion 
theory, since it would not add anything that we cannot conclude from Figs. [TO|and[TT| 

Finally, we need to compare the predictions with the empirical phase shifts. In F waves the N'^LO 
predictions above 200 MeV are, in general, too attractive. Note, however, that this is also true for the 
predictions by the Bonn tt + 2tt model. In the full Bonn model, besides tt -|- 2tt, (repulsive) lu and 
exchanges are included which move the predictions right on top of the data. The exchange of a a; meson or 
combined irp exchange are Stt exchanges. Three-pion exchange occurs first at chiral order four. It has been 
investigated by Kaiser |S71|SH] and found to be negligible, at this order. However, 3n exchange at order five 
appears to be sizable |59j and may have impact on F waves. Besides this, there is the usual short-range 
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Figure 11: Same as Fig. |10[ but for G waves. 



phenomenology. In ChPT, this short-range interaction is parametrized in terms of four-nucleon contact 
terms (since heavy mesons do not have a place in that theory). Contact terms of order four (N'^LO) do 
not contribute to F-waves, but order six does. In summary, the remaining small discrepancies between the 
N'^LO predictions and the empirical phase shifts may be straightened out in fifth or sixth order of ChPT. 

4-3. NN contact terms 

The successful test of the chiral pion-exchange contributions in peripheral partial waves in the previous 
section has shown that we got the long- and intermediate-range parts of the nuclear force right. However, for 
a "complete" nuclear force, we have to describe correctly all partial waves, including the lower ones. In fact, 
in calculations of NN observables at low energies (cross sections, analyzing powers, etc.), the lower partial 
waves with L < 2 are the most important ones, generating the largest contributions. The same is true for 
microscopic nuclear structure calculations. These lower partial waves are dominated by the dynamics at 
short distances. Therefore, we need to look now into the short-range part of the NN potential. 

In conventional meson theory, the short-range nuclear force is described by the exchange of heavy mesons, 
notably the aj(782). Qualitatively, the short-distance behavior of the NN potential is obtained by Fourier 
transform of the propagator of a heavy meson, 

piq-r p-m^r 

d'l^-^ • (4-32) 

ChPT is an expansion in small momenta Q, too small to resolve structures like a p(770) or w(782) meson, 
because Q « Tnp.u- But the latter relation allows us to expand the propagator of a heavy meson into 

a power series, 

-J—^\(l-% + 9^^ + ,\^ (4.33) 
mi + mi \ mi m% J 
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where the w is representative for any heavy meson of interest. The above expansion suggests that it should 
be possible to describe the short distance part of the nuclear force simply in terms of powers of Q /m^, which 
fits in well with our over-all power expansion since Q/rriu, ~ Q/A^. 

A second purpose of contact terms is renormalization. Dimensional regularization of the loop integrals 
that occur in multi-pion exchange diagrams typically generates polynomial terms with coefRcients that are, 
in part, infinite or scale dependent (cf. Appendix B). Contact terms pick up infinities and remove scale 
dependences, which is why they are also known as counter terms. 

The partial-wave decomposition of a power has an interesting property. First note that Q can only 
be either the momentum transfer between the two interacting nucleons q or the average momentum k [cf. 
Eq. (4.81 for their definitions]. In any case, for even i^, 

Q-^/llcos^), (4.34) 

where fm stands for a polynomial of degree m and 9 is the CMS scattering angle. The partial-wave 
decomposition of Q'^ for a state of orbital-angular momentum L involves the integral 



Q''PL{cos9)dcoa9 



/id (cos 9)P]^ (cos 9)d cos 9 , 



where is a Legendre polynomial. Due to the orthogonality of the P^, 



for 



(4.35) 



(4.36) 



Consequently, contact terms of order zero contribute only in S'-waves, while order-two terms contribute up 
to P- waves, order- four terms up to £)-waves, etc.. 

Due to parity, only even powers of Q are allowed. Thus, the expansion of the contact potential is formally 
given by 



(2) 



(4) 



(4.37) 

where the superscript denotes the power or order. 

We will now present, one by one, the various orders of NN contact terms which result from the contact 
Lagrangians presented in Section [2.3| 



.3.1. Zeroth order (LO) 



The contact Lagrangian -C^jy, Eq. (|2.61|), which is part of C 
contact potential, 



Eq. ( |2.65[ ), leads to the following NN 

(4.38) 



and, in terms of partial waves, we have 



(0),1 



^o) 



Ciso =47r(C5-3CT) 
= 47r {Cs + Ct) ■ 



(4.39) 



.3.2. Second order (NLO) 

The contact Lagrangian -CjvAf (|2-62[), which is part of C'^~^, Eq. (|2.67|, generates the following NN 



contact potential 



+ (Cs + Ci k^) (71 • (72 
+ C^i^-lS-iqxk) 

+ Ce ((71 ■ q){S2- q) 
+ C'r ((71 • k) ((72 • k) . 
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(4.40) 



The coefficients Ci used here i n the contact potential are, of course, related to the coefficients C[ that occur 
in the Lagrangian C)^^, Eq. (2.62 1. The relation, which is unimportant for our purposes, can be found in 
Refs. [illMj- 

There are many ways to perform the partial-wave decomposition of the above potential. We perceive 
the method presented by Erkelenz, Alzetta, and Holinde |142j as the most elegant one. Thus, one obtains 
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.3.3. Fourth order (N^LO) 
The contact potential of order four reads 



3 



+ (Ds + £'6 fc"* + D7 + ZJg (g X ) CTi • CT2 

+ (Z?ii g2 + i?,2 e) (a, . • q) 

+ (Z?i3 + I?14 k^) ■ k) (g2 ■ k) 

+ Di5 ( (7i ■ {qxk) $2- [q-x k) 



(4.41) 



(4.42) 



(4.43) 



The rather lengthy partial- wave expressions of this order are relegated to [Appendix E 
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4.4- Definition of NN potential 

We have now rounded up everything needed for a reahstic nuclear force — long, intermediate, and short 
ranged components — and so we can finally proceed to the lower partial waves. However, here we encounter 
another problem. The two-nucleon system at low angular momentum, particularly, in S waves, is character- 
ized by the presence of a shallow bound state (the deuteron) and large scattering lengths. Thus, perturbation 
theory does not apply. In contrast to tt-tt and n-N, the interaction between nucleons is not suppressed in 
the chiral limit (Q — > 0). Weinberg [SHI EI] showed that the strong enhancement of the scattering amplitude 
arises from purely nucleonic intermediate states ("infrared enhancement"). He therefore suggested to use 
perturbation theory to calculate the NN potential (i.e., the irreducible graphs) and to apply this potential 
in a scattering equation to obtain the NN amplitude. We will follow this prescription and discuss potential 
problems in the next subsection. 

Since the irreducible diagrams that make up the potential are calculated using covariant perturbation 



theory (cf. Section 4.1), it is consistent to start from the covariant Bethe-Salpeter (BS) equation |143) 



describing two-nucleon scattering. In operator notation, the BS equation reads 

T = V + VgT (4.44) 

where T denotes the invariant T-matrix (T — iM with Ai the invariant amplitude) for the two-nucleon 
scattering process, V the sum of all connected two-particle irreducible diagrams, and Q is (— i) times the 
relativistic two-nucleon propagator. The BS equation is equivalent to a set of two equations 

T ^ V + VgT (4.45) 
V = V + V{g-g)V (4.46) 
= V + Vi, (^ - g) Vi, -t- . . . , (4.47) 

where g is a, covariant three-dimensional propagator which preserves relativistic elastic unitarity. We choose 
the propagator g proposed by Blankenbecler and Sugar (BbS) |133) (for more details on relativistic three- 



dimensional reductions of the BS equation, see Ref. [H]). The ellipsis in Eq. (4.471 stands for terms of 
irreducible Stt and higher pion exchanges which we neglect. 

Note that when we speak of covariance in conjunction with (heavy baryon) ChPT, we are not referring 
to manifest covariance. Relativity and relativistic off-shell effects are accounted for in terms of a Q/Mm 
expansion up to the given order and up to the number of pions we take into consideration. Thus, Eq. ( |4.47[ ) 
is evaluated in the following way, 

V « V(on-shell) + Vi, G Vi, - V^^ g V^^ (4.48) 
~ Vi. + + Vi^gVi, -yi,.9^i. (4.49) 
= Vi. + VL , (4.50) 



where Vi^r denotes the relativistic (off-shell) IPE, while Vi^ is the on-shell IPE given in Eq. (4.5 ). V2,r stands 



for the irreducible 27r exchanges calculated (on-shell) in Section 4.1.2 but without the covariant planar box 



diagram, which in the above equations is given by ViT^QVi-n- Furthermore, Vi^gVi^ is the iterated IPE 



discussed in Section 4.1.2 Thus, the term (Vi^ Q Vi^r — Vi^ g Vi^) represents the irreducible part of the box 
diagram contribution. Finally, Vj-n- subsumes all 27r exchanges without the iterated Vi^r and is, therefore, 
also known as the irreducible 2PE. Since all contributions are calculated on shell, the potential has no energy 
dependence. 

Adding the short-range contact terms Vet to the above, yields the full TV TV potential, 

V^Vi^ + V^^ + . (4.51) 



Notice that the pion-exchange part of this potential differs from the perturbative amplitude, Eq. (4.30), 
by the absence of the iterative 2PE. The latter is generated automatically when the potential is inserted 
into a Schroedinger or Lippmann-Schwinger equation (see below). We also note that adding the contact 
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interactions to the pion interactions will generate loop corrections of the contact terms. This leads to a 
m,r-dependent renormalization of the contact parameters and can be ignored if mjr-dependence is not an 
issue. 

As discussed, the irreducible 2PE, V^^^, is organized according to increasing orders, 



'(2) 



and was calculated in Section 



4.1.2 



/(2) 



V, 



'(3) 



V, 



'(4) 



(4.52) 

^2^^'' is given by the contributions of Eqs. (4.9) and (4.10), ¥2^"^ is 



made up from Eqs. (|4.13|)-(|4.24|), and Vj^^'' is contained in Appendix D 



The contact potentials come in even orders, 



Vr t - K 
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(4) 



and were presented in Section |4.3[ 

In summary, the NN potential V, calculated to certain orders, is given by: 
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The potential V satisfies the relativistic BbS equation, Eq. (4.45), which reads explicitly, 



V{p " 



1 



np",p) 



(4.53) 



(4.54) 
(4.55) 
(4.56) 
(4.57) 



(4.58) 



with Epii = 



p' 



(27r)3 ■ '-f ' Epfi p^ - p" 
. The advantage of using a relativistic scattering equation is that it automatically 



includes relativistic corrections to all orders. Thus, in the scattering equation, no propagator modifications 
are necessary when raising the order to which the calculation is conducted. 
Defining 



and 



V(j>',p) = 



T{p\p)^ 




V(j>',p) 



T{p',p) 



N 



En 



(4.59) 



(4.60) 



where the factor l/(27r)'^ is added for convenience, the BbS equation collapses into the usual, nonrelativistic 
Lippmann-Schwinger (LS) equation, 



T{p',p} - V{p',p) + / dV V{p',p") 



M. 



N 



p2 _ p//- _|_ jg 



T{p",p) 



(4.61) 



Since V satisfies Eq. (4.61 ), it can be used like a usual nonrelativistic potential, and T may be perceived as 
the conventional nonrelativistic T-matrix. In applications, it is more convenient to use the K-matrix instead 
of the T-matrix and to have the LS equation decomposed into partial waves: all these technical issues are 
explained in detail in Appendix A of Ref. (13j where also the formulas for the calculation of np and pp (the 
latter with Coulomb) phase shifts are provided. The partial wave decomposition of the operators by which 
the potential is represented can be found in section 4 of Ref. |142) , and numerical methods for solving the 
LS equation are explained in Ref. |144j . 
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.5. Renormalization 

.5.1. Regularization and non-perturbativ e reno rmalization 
Iteration of V in the LS equation, Eq. (4.61 1, requires cutting V off for high momenta to avoid infinities. 



This is consistent with the fact that ChPT is a low-momentum expansion which is vahd only for momenta 
Q ^ « 1 GeV. Therefore, the potential V is multiplied with the regulator function f{p',p), 



with 

such that 



V{p',p)^V{p',p)f{p',p) 
/(p',p)=exph(p7A)2«-(p/A)2" 

V{p',p)fip',p)^Vip',p)\l- 



2n 



2n 



Typical choices for the cutoff parameter A that appears in the regulator are A 



(4.62) 
(4.63) 

(4.64) 



0.5 GeV < A^ « 1 GeV. 



Equation (4.64) provides an indication of the fact that the exponential cutoff docs not necessarily affect 
the given order at which the calculation is conducted. For sufficiently large n, the regulator introduces 
contributions that are beyond the given order. Assuming a good rate of convergence of the chiral expansion, 
such orders are small as compared to the given order and, thus, do not affect the accuracy at the given order. 



In calculations, one uses, of course, the exponential form, Eq. (4.63), and not the expansion Eq. (4.64). On a 
similar note, we also do not expand the square-root factors in Eqs. ( 4.59|[4^60 ) because they are kinematical 
factors which guarantee relativistic elastic unitarity. 

It is pretty obvious that results for the T-matrix may depend sensitively on the regulator and its cutoff 
parameter. This is acceptable if one wishes to build models. For example, the meson models of the 
past [TTl [2U] always depended sensitively on the choices for the cutoff parameters which, in fact, were 
important for the fit of the NN data. However, the EFT approach wishes to be fundamental in nature and 
not just another model. 

In field theories, divergent integrals are not uncommon and methods have been developed for how to deal 
with them. One regulates the integrals and then removes the dependence on the regularization parameters 
(scales, cutoffs) by renormalization. In the end, the theory and its predictions do not depend on cutoffs or 
renormalization scales. 

So-called renormalizable quantum field theories, like QED, have essentially one set of prescriptions that 
takes care of renormalization through all orders. In contrast, EFTs are renormalized order by order. 

The renormalization of perturbative EFT calculations is not a problem. The problem is nonperturbative 
renormalization. This problem typically occurs in nuclear EFT because nuclear physics is characterized by 
bound states which are nonperturbative in nature. EFT power counting may be different for nonperturbative 
processes as compared to perturbative ones. Such difference may be caused by the infrared enhancement of 
the reducible diagrams generated in the LS equation. 

Weinberg's implicit assumption |47l I50| was that the counterterms introduced to renormalize the per- 
turbatively calculated potential, based upon naive dimensional analysis ("Weinberg counting"), are also 
sufficient to renormalize the nonperturbative resummation of the potential in the LS equation. In 1996, 
Kaplan, Savage, and Wise (KSW) |145j pointed out that there are problems with the Weinberg scheme if 
the LS equation is renormalized by minimally-subtracted dimensional regularization. This criticism resulted 
in a fiurry of publications on the renormalization of the nonperturbative NN problem |146H162] . The lit- 
erature is too comprehensive to discuss all contributions. Let us just mention some of the work that has 
particular relevance for our present discussion. 

If the potential V consists of contact terms only (a.k.a. pion-less theory), then the nonperturbative 
summation (4.61) can be performed analytically and the power counting is explicit. However, when pion 



exchange is included, then (4.61 ) can be solved only numerically and the power counting is less transparent. 



Perturbative ladder diagrams of arbitrarily high order, where the rungs of the ladder represent a potential 
made up from irreducible pion exchange, suggest that an infinite number of counterterms is needed to 
achieve cutoff independence for all the terms of increasing order generated by the iterations. For that 
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reason, Kaplan, Savage, and Wise (KSW) |145) proposed to sum the leading-order contact interaction to 
all orders (analytically) and to add higher-order contacts and pion exchange perturbatively up to the given 
order. Unfortunately, it turned out that the order by order convergence of IPE is poor in the ^Si-^Di 
state |146j . The failure was triggered by the singularity of the IPE tensor force when iterated to 

second order. Therefore, KSW counting is no longer taken into consideration (see, however, [160| ). A 
balanced discussion of possible solutions is provided in 150 . 

Some researchers decided to take a second look at Weinberg's original proposal. A systematic inves- 
tigation of Weinberg counting in leading order has been conducted by Nogga, Timmermans, and van 
Kolck ^152) in momentum space, and by Valderrama and Arriola at LO and higher orders in configura- 
tion space |151l 11531 1154) . A comprehensive discussion of both approaches and their equivalence can be 
found in Refs. [rFflfTOT] . 



The LO NN potential is given in (4.54[) and consists of IPE plus two nonderivative contact terms that 



contribute only in S waves. Nogga et al. jl52] find that the given counterterms renormalize the S waves 
(i.e., stable results are obtained for A — >■ oo) and the naively expected infinite number of counterterms is 
not needed. This means that Weinberg power counting does actually work in S waves at LO (ignoring the 
dependence of the contact interaction discussed in Refs. |1451I150] ). However, there are problems with a 
particular class of higher partial waves, namely those in which the tensor force from IPE is attractive. The 
first few cases of this kind of low angular momentum are ^Pq, ^P2, and ■^I?2, which need a counterterm for 
cutoff independence. The leading (nonderivative) counterterms do not contribute in P and higher waves, 
which is why Weinberg counting fails in these cases. But the second order contact potential provides 
counterterms for P waves. Therefore, the promotion of, particularly, the '^Pq and contacts from NLO to 
LO would fix the problem in P waves. To take care of the problem, a N^^LO contact, i.e. a term from 
V^^\ needs to be promoted to LO. Partial waves with orbital angular momentum L > 3 may be calculated 
in Born approximation with sufficient accuracy and, therefore, do not pose renormalization problems. In 
this way, one arrives at a scheme of 'modified Weinberg counting' |152j for the leading order two-nucleon 
interaction. 

4-5.2. Renormalization beyond leading order 

As shown below, for a quantitative chiral NN potential one needs to advance all the way to N'^LO. 
Thus, the renormalization issue needs to be discussed beyond LO. Naively, the most perfect renormalization 
procedure is the one where the cutoff parameter A is carried to infinity while stable results are maintained. 
This was done successfully at LO in the work by Nogga et al. |152j described above. At NNLO, the 
infinite-cutoff renormalization procedure has been investigated in [158) for partial waves with total angular 
momentum J < 1 and in |154) for all partial waves with J < 5. At N^LO, an investigation of the ^5*0 state 
exists |157j . From all of these works, it is evident that no counter term is effective in partial- waves with 
short-range repulsion and only a single counter term can effectively be used in partial- waves with short-range 
attraction. Thus, for the A — > c» renormalization prescription, even at N'^LO, we have either one or no 
counter term per partial-wave state. This is inconsistent with any reasonable power-counting scheme and, 
therefore, defies the principals of an EFT. 

A possible way out of this dilemma was proposed already in jl52j and reiterated in a recent paper by 
Long and van Kolck |159j . In the latter reference, the authors examine the renormalization of an attractive 
1/r^ potential perturbed by a l/r^ correction. Generalizing their findings, they come to the conclusion 
that, for any attractive 1/r" potential (with n > 2), partial waves with low angular momentum L must 
be summed to all orders and one contact term is needed for each L to renormalize the LO contribution. 
However, there exists an angular momentum Lp {Lp w 3 for the nuclear case, cf. Ref. [152| ). above which 
the leading order can be calculated perturbatively. In short, naive dimensional analysis (NDA) does not 
apply at LO below Lp. However, once this failure of NDA is corrected at LO, higher order corrections can 
be added in perturbation theory using counterterm that follow NDA j l59| . 

Reference |159j used just a toy model and, therefore, a full investigation using the chiral expansion is 
needed to answer the question if this renormalization approach will work for the realistic nuclear force. 
A first calculation of this kind for the S waves was recently performed by Valderrama |163| . The author 
renormalizes the LO interaction nonperturbatively and then uses the LO distorted wave to calculate the 2PE 
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contributions at NLO and NNLO perturbatively. It turns out that perturbative renormalizability requires 
the introduction of three counterterms in ^Sq and six in the coupled ^Si Di channels. Thus, the number 
of counterterms required in this scheme is larger than in the Weinberg scheme, which reduces the predictive 
power. For a final evaluation of this approach, also the results for P and D waves are needed, which are not 
yet published. 

However, even if such a project turns out to be successful for NN scattering, there is doubt if the 
interaction generated in this approach is of any use for applications in nuclear few- and many-body problems. 
In applications, one would first have to solve the many-body problem with the renormalized LO interaction, 
and then add higher order corrections in perturbation theory. However, it was shown in a recent paper |164) 
that the renormalized LO interaction is characterized by a very large tensor force from IPE. This is no 
surprise since LO is renormalized with A — >■ (X) implying that the IPE, particularly its tensor force, is totally 
uncut. As a consequence of this, the wound integral in nuclear matter, k, comes out to be about 40%. The 
hole-line and coupled cluster expansions are know to converge oc with n the number of hole-lines or 

particles per cluster [165H167] . For conventional nuclear forces, the wound integral is typically between 5 
and 10% and the inclusion of three-body clusters (or three hole-lines) are needed to obtain converged results 
in the many-body system [551 [571 [T55] . Thus, if the wound integral is 40%, probably, up to six hole-hnes 
need to be included for approximate convergence. Such calculations are not feasible even with the most 
powerful computers of today and will not be feasible any time soon. Therefore, even if the renormalization 
procedure proposed in |159j will work for NN scattering, the interaction produced will be highly impractical 
(to say the least) in applications in few- and many-body problems because of convergence problems with 
the many-body energy and wave functions. 

4-5.3. Back to the beginnings 

The various problems with the renormalization procedures discussed above may have a simple common 
reason: An EFT that has validity only for momenta Q < A^ is applied such that momenta Q ^ A^^ are 
heavily involved (because the regulator cutoff A — ?> oo). A recent paper by Epelbaum and Gegelia |169) 
illustrates the point: The authors construct an exactly solvable toy-model that simulates a pionful EFT 
and yields finite results for A — )• oo. However, as it turns out, these finite results are incompatible with the 
underlying EFT, while for cutoffs in the order of the hard scale consistency is maintained. In simple terms, 
the point to realize is this: // an EFT calculation produces (accidentally) a finite result for A — >■ oo, then 
that does not automatically imply that this result is also meaningful. 

This matter is further elucidated in the lectures by Lepage of 1997 [170]. Lepage points out that it makes 
little sense to take the momentum cutoff beyond the range of validity of the effective theory. By assumption, 
our data involves energies that are too low — wave lengths that are too long — to probe the true structure 
of the theory at very short distances. When one goes beyond the hard-scale of the theory, structures are 
seen that are almost certainly wrong. Thus, results cannot improve and, in fact, they may degrade or, in 
more extreme cases, the theory may become unstable or untunable. In fact, in the NN case, this is what 
is happening in several partial waves (as reported above). Therefore, Lepage suggests to take the following 
three steps when building an effective theory: 

1. Incorporate the correct long-range behavior: The long-range behavior of the underlying theory must 
be known, and it must be built into the effective theory. In the case of nuclear forces, the long-range 
theory is, of course, well known and given by one- and multi-pion exchanges. 

2. Introduce an ultraviolet cutoff to exclude high-momentum states, or, equivalent, to soften the short- 
distance behavior: The cutoff has two effects: First it excludes high-momentum states, which are 
sensitive to the unknown short-distance dynamics; only states that we understand are retained. Second 
it makes all interactions regular at r = 0, thereby avoiding the infinities that beset the naive approach. 

3. Add local correction terms (also known as contact or counter terms) to the effective Hamiltonian. These 
mimic the effects of the high-momentum states excluded by the cutoff introduced in the previous step. 
In the meson-exchange picture, the short-range nuclear force is described by heavy meson exchange, 
like the p(770) and w(782). However, at low energy, such structures are not resolved. Since we must 
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Figure 12: Phase shifts of np scattering as calculated from NN potentials at different orders of ChPT. The black dotted 
line is LO(500), the blue dashed is NLO(550/700) [T72], the green dash-dotted NNLO(600/700) [172], and the red solid 
N^LO(500) '68], where the numbers in parentheses denote the cutoffs in MeV. Partial waves with total angular momentum 
J < 1 arc displayed. Empirical phase shifts (solid dots and open circles) as in Fig. [s] 



include contact terms anyhow, it is most efficient to use them to account for any heavy-meson exchange 
as well. The correction terms systematically remove dependence on the cutoff. 

A first investigation in the above spirit has been conducted by Epelbaum and Meifiner [155| in 2006. 
The authors stress that there is no point in taking the cutoff A beyond the breakdown scale of the EFT, 
A-^ « Trip « 1 GeV, since the error of the calculation is not expected to decrease in that regime. Any value 
for the cutoff parameter A is acceptable if the error associated with its finite value is within the theoretical 
uncertainty at the given order. The authors conduct an investigation at LO (including only the counter 
terms implied by Weinberg counting) and find that, starting from A w 3 fm~^, the error in the NN phase 
shifts due to keeping A finite stays within the theoretical uncertainty at LO. 

4-5.4- Concluding the renormalization issue 

Crucial for an EFT are regulator independence (within the range of validity of the EFT) and a power 
counting scheme that allows for order-by-order improvement with decreasing truncation error. The purpose 
of renormalization is to achieve this regulator independence while maintaining a functional power counting 
scheme. After the comprehensive tries and errors of the past, it appears that there are two renormaliza- 
tion schemes which have the potential to achieve the above goals and, therefore, should be investigated 
systematically in the near future. 

In scheme one, the LO calculation is conducted nonperturbatively (with A — > oo as in |152| ) and sublead- 
ing orders are added perturbatively in distorted wave Born approximation. As mentioned above, Valderrama 
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Figure 13: Same as Fig. |12[ but J = 2 phase shifts and J < 2 mixing parameters are shown. 



has started this in 5* waves |163j , but results in higher partial waves are needed to fully assess this approach. 
Even though at this early stage any judgment is speculative, we take the liberty to predict that this approach 
will be only of limited success and utility — for the following reasons. First, it will probably require about 
twice as many counterterms as Weinberg counting and, therefore, will have less predictive power. Second, 
this scheme may converge badly, because the largest portion of the nuclear force, namely, the intermediate- 
range attraction appears at NNLO. Third, as discussed in Ref. |164| . this force may be problematic (and, 
therefore, impractical) in applications in nuclear few- and many-body systems, because of a pathologically 
strong tensor force that will cause bad convergence of energy and wave functions. Finally, in the work that 
has been conducted so far within this scheme by Valderrama, it is found that only rather soft cutoffs can be 
used. 

The latter point (namely, soft cutoffs) suggests that one may then as well conduct the calculation 
nonperturbatively at all orders (up to N'^LO) using Weinberg counting, which is no problem with soft 
cutoffs. This is scheme two that we propose to investigate systematically. In the spirit of Lepage, the 
cutoff independence should be examined for cutoffs below the hard scale and not beyond. Ranges of cutoff 
independence within the theoretical error are to be identified using 'Lepage plots' |170j . A very systematic 
investigation of this kind does not exist at this time and is therefore needed. However, there is comprehensive 
circumstantial evidence from the numerous chiral A^A'' potentials constructed over the past decade [53 1^51 - 



I171j (see Figs. 14 and 15 below) indicating that this investigation will most likely be a success. The 



potentials discussed in the following section are all based upon Weinberg counting. 
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Table 3: Columns three and four show the x^/datum for the reproduction of the 1999 np database [173] (subdivided into 
energy intervals) by families of np potentials at NLO and NNLO constructed by the Juelich group 11721 . The x'^ /datum is 



stated in terms of ranges which result from a 


variation of the cutoff parameters used in the regulator functions. The values of 


these cutoff parameters 


in units of MeV are 


given in parentheses. Ti^b denotes the kinetic ener 


5y of the incident neutron in 


the laboratory system. 








Tiab bin 


# of np 


— Juelich np potentials — 


(MeV) 


data 


NLO 


NNLO 






(550/700-400/500) 


(600/700-450/500) 


0-100 


1058 


4-5 


1.4-1.9 


100-190 


501 


77-121 


12-32 


190-290 


843 


140-220 


25-69 


0-290 


2402 


67-105 


12-27 



4-6. Constructing quantitative chiralNN potentials 
4.6.1. What order? 

As discussed, the NN potential can be calculated up to various orders, cf. Eqs. ( 4.54 )-( 4.57), and the 
accuracy increases as the order increases. That triggers the obvious question: To what order do we have 
to go for an accuracy that we would perceive as necessary and sufficient for, e. g., reliable microscopic 
nuclear structure calculations? We had a first look at the order by order evolvement in Section |4.2| where 
we compared phase shift predictions with empirical values in peripheral partial waves. But that was only 
a small part of the story. The most important partial waves are the lower ones, since they generate the 
largest contributions in most applications. Now we have potentials at hand that are defined through all 
partial waves, and so we can investigate the issue of how well ChPT converges in the important lower partial 
waves. This is demonstrated in Figs. [12] and [T3j where we show the J < 2 phase parameters for potentials 
constructed at LO, NLO, NNLO, and N'^LO. These figures clearly reveal substantial improvements in the 
reproduction of the empirical phase shifts as the orders go up. 

There is an even better way to confront theory with experiment. We can calculate observables of NN 
scattering and compare to the experimental data. It is customary to state the result of such a comparison 
in terms of the /datum where a value around unity would signify a perfect fit. 

Let's start with potentials developed to NLO and NNLO. In Table [s] we show the x^/datum for the fit 
of the world np data below 290 MeV for families of np potentials at NLO and NNLO constructed by the 
Juelich group |172j . The NLO potentials produce the very large x^/datum between 67 and 105, and the 
NNLO are between 12 and 27. The rate of improvement from one order to the other is very encouraging, 
but the quality of the reproduction of the np data at NLO and NNLO is obviously insufficient for reliable 
predictions. 

Based upon these facts, it has been pointed out in 2002 by Entem and Machleidt [HZ] that one has 
to proceed to N'^LO. Consequently, the first N'^LO potential was published in 2003 |68) . 

At N^^LO, there are a total of 24 contact terms (24 parameters) which contribute to the partial waves with 
L < 2 (cf. Section 4.3 and Appendix E). These 24 LECs are essentially free constants which parametrize the 
short-ranged phenomenological part of the interaction. In Table [4] column 'Q^/N'^LO', we show how these 
terms are distributed over the partial waves. Most important for the improved reproduction of the NN 
phase shifts (and NN observables) at N^^LO is the fact that contacts appear for the first time in D-waves. 
D-waves are not truly peripheral and, therefore, IPE plus 2PE alone do not describe them well (Fig. 13). 
The D-wave contacts provide the necessary short-range corrections to get the Z?-phases right. Besides this, 
at N^LO, another contact is added to each P-wave, which leads to substantial improvements, particularly, 
in and ^Pi above 100 MeV (cf. Fig. [l2|. 

In Table|4j we also show the number of parameters used in the Nijmegen partial wave analysis (PWA93) |139) 
and in the high-precision CD-Bonn potential [13]. The table reveals that, for S and P waves, the number 
of parameters used in high-precision phenomenology and in EFT at N'^LO are about the same. Thus, the 
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Table 4: Number of parameters needed for fitting tiie np data in phase-shift analysis and by a high-precision A^A^ potential 
versus the total number of A''A'' contact terms of EFT based potentials to different orders. 
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EFT approach provides retroactively a justification for the phenomenology used in the 1990's 
to obtain high-precision fits. 

At NLO and NNLO, the number of contact parameters is substantially smaller than for PWA93 and 
CD-Bonn, which explains why these orders are insufficient for a quantitative potential. The 24 parameters 
of N'^LO are close to the 30-1- used in PWA93 and high precision potentials. Consequently (see following 
sections for details), at N'^LO, a fit of the NN data is possible that is of about the same quality as the one 
by the high-precision NN potentials p ^ lT3| .fT74|ll75] . Thus, one may perceive N^LO as the order of ChPT 
that is necessary and sufficient for a reliable TV TV potential. 

4-6.2. Charge- dependence 

So far we considered only neutron-proton scattering. However, as stressed repeatedly by the Nijmegen 
group [1391 1176) . for a precise reproduction of the TV TV data, i.e., the np and the pp data, it is crucial to 
take charge-dependence into account and to construct separate np, pp, and nn potentials which differ by 
the subtleties of charge-dependence. 

By definition, charge independence (or isospin symmetry) is invariance under any rotation in isospin 
space. A violation of this symmetry is referred to as charge dependence or charge independence breaking 
(CIB). Charge symmetry is invariance under a rotation by 180*^ about the y-axis in isospin space if the 
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positive 2-direction is associated with the positive charge. The violation of this symmetry is known as 
charge symmetry breaking (CSB). Obviously, CSB is a special case of charge dependence. 

GIB of the strong NN interaction means that, in the isospin 1=1 state, the proton-proton {1^ = 
+1), neutron-proton (/^ = 0), or neutron-neutron (/^ = —1) interactions are (slightly) different, after 
electromagnetic effects have been removed. CSB of the NN interaction refers to a difference between 
proton-proton (pp) and neutron-neutron (nn) interactions, only. For reviews, see Refs. PH, 11041 11371 11771 - 



CIB is seen most clearly in the ^Sq NN scattering lengths. The latest empirical values for the singlet 
scattering length a and effective range r are: 

= -17.3 ± 0.4 fm [nS], = 2.85 ± 0.04 fm [178]; (4.65) 

= -18.95 ± 0.40 fm [mi [IM], r^„ = 2.75 ± 0.11 fm [HH]; (4.66) 
a„p = -23.740 ± 0.020 fm [13], r„p = 2.77 ± 0.05 fm [13]. (4.67) 

The values given for pp and nn scattering refer to the nuclear part of the interaction as indicated by the 
superscript N; i. e., electromagnetic effects have been removed from the experimental values. 
The above values imply that charge-symmetry is broken by 

AacsB = a^p -a^„ = 1.65 ± 0.60 fm, (4.68) 
ArcsB = r^p -r^„= 0.10 ± 0.12 fm; (4.69) 

and the following CIB is observed: 

AaciB = \ (<p + a^n)- anp = 5.6 ± 0.6 fm, (4.70) 
Arc/B= ^(rfp + r,^„)-r„p = 0.03 ± 0.13 fm. (4.71) 

In summary, the TV TV singlet scattering lengths show a small amount of CSB and a clear signature of CIB. 

The current understanding is that — on a fundamental level — charge-dependence ('isospin violation') is 
due to a difference between the up and down quark masses and electromagnetic interactions. As first 
discussed in Refs. [183H185] . EFT is a suitable tool to also deal with isospin violations. In the two-flavor 



case, the mass term in the QCD Lagrangian, Eq. (2.1|, can be written as [cf. Eq. (2.25l] 



with 



qMq = -^q {niu + rud) I q - {mu - md) T-i q 

= -^q(TO„ +TOrf) (/- eT3)gf (4.72) 

^^ m.-m„ ^l^ (4.73) 
niu + md 3 



cf. Eqs. (2.4) and (2.5). The first term of Eq. (4.72) conserves isospin but breaks chiral symmetry. It 
is responsible for the nonvanishing pion mass and leads to terms oc (to^)" (with n > 1) in the effective 
Lagrangian. The second term breaks isospin symmetry and generates terms proportional to (em^)". Note 
that the isospin-breaking effects are much smaller than the value for e suggests, because the relevant scale 
is rather than (tou -|- m^). 

For isospin violating effects caused by the electromagnetic interaction, the small parameter = Ana 
1/10 can be used (where a = 1/137.036 denotes the fine-structure constant). It is then suggestive to consider 

e^e^^ (4.74) 
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Table 5: Isospin breaking contributions to the A'^A'^ interaction 



Order 


o n f 1 Vi 1 1 1 o n d 

\ '\Jl.l.\il. ± 1^ UL Li IVJllo 




NL0 (v — 2) 


_Pion-mass splitting m IPb, 
Static Coulomb potential. 




NNL0 {v = -i) 


CSB contacts without derivatives, 






Charge-dependence of the pion-nucleon coupling constant in IPE 


^ em2/A2). 




CIB contacts without derivatives. 






Charge-dependence of the pion-nucleon coupling constant in IPE [~ 


eV(4^)2], 




Pion-mass splitting in NLO 2PE, 






Nucleon-mass splitting in NLO 2PE and LS-equation, 






7r7-exchange, 






Relativistic corrections to the Coulomb potential Q^/Af^), 






Further electromagnetic corrections. 





as the general expansion parameter. For photon loops, one may further assume 



(47r)2 



(4.75) 



We will now briefly discuss various isospin violating contributions to the NN interaction listed in Table [5j 
For a more comprehensive and systematic study, the interested reader is referred to Ref. [186) . 

It is well known that the pion-mass difference is mainly due to electro-magnetic interactions among the 
quarks and, thus, according to Eq. (4.751, of order = 4. Therefore, the effect of pion-mass difference on 
IPE can be estimated to be 



Am: 



A2 

X 



A2 

X 



A2 

X 



(47r)2 



A2 

X 



(4.76) 



So, it is of order v — 2 or NL0 (cf. Table [5]) with the slash signifying that the extended power counting 
scheme, Eqs. (4.74) and (4.751, is applied. This CIB effect is taken care of by replacing the 7=1 isospin- 
symmetric IPE potential, Eq. (4.5), by 



T/(PP) _ T/(n") 



Vi^(m^o) 



for and nn scattering and 



T/(np,-f=l) 



(4.77) 



(4.78) 



for I = 1 np scattering, with Vi^(m^) defined in Eq. (4.27). If the pion masses were all the same, the above 
potentials would be identical. However, due to their mass splitting (Table [2]), the I — 1 np potential is more 
attractive than the pp one, and this effect is known to explain about one half of the CIB scattering-length 
difference, Eq. (4.70). For completeness, we also give the / = np IPE potential which is 



.0) 



(4.79) 



Due to the smallness of the pion mass, IPE is also a sizable contribution in all partial waves with L > 0; 
and due to the pion's relatively large mass splitting (3.4%), IPE creates relatively large charge-dependent 
effects in all partial waves (cf. Tables V and VI and Fig. 4 of Ref. [E]). Therefore, all modern phase shift 
analyses [1391 ED] and all modern NN potentials |12l [13l [174] include the CIB effect created by IPE. 

The other NL0 contribution (Table [5]) is due to the static Coulomb interaction which is ~ and, 
therefore, also of order two. The Coulomb potential has considerable impact on the ^5*0 phase shifts at low 
energies and creates CIB effects of similar size as IPE in P and D waves (cf. Tables V and VI and Fig. 5 of 
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Ref. [H]). The calculation of phase shifts in the presence of the Coulomb potential is explained in Appendix 
A.3 of Ref. p]. 

At NNL0 (j/ = 3), CSB contributions to the charge-dependence of the pion-nucleon coupling constant 
em^/A^) occur. However, since empirically there is no clear evidence for such charge-dependence [10411132] . 
this contribution is ignored. Moreover, at NNL0, there is a non-derivative CSB contact term 

cx eml{NT3N){NN) , (4.80) 



which is crucial for the fit of the CSB splitting of the ^5*0 scattering length, Eq. (4.68) 



The leading CIB contact interactions are of electromagnetic origin and have the structure 



=2 



^ iNT3N){NT3N) . (4.81) 

They are obviously of order v = A (N^L0) and needed for the fit of the CIB splitting of the ^S^ scattering 
length, Eq. (4.701. Besides this, there are numerous other N'^L0 contributions (cf. Table [5]). One-loop 



electromagnetic corrections to the pion-nucleon coupling constant create effects of order e^(47r)^ in the 
IPE, which will be ignored for the reasons given above. Corrections due to pion-mass difference in the NLO 
2PE can be estimated to be 

^mlQ^_Q^ (4.82) 



As explained in Refs. |187| I188j . these corrections can be calculated most conveniently by separating the 



NLO 2PE as given in Eqs. (4.9) and (4.10) into an isoscalar and an isovector piece: 

V2. = V^, +T,-T2 Vl . (4.83) 

The isoscalar part ¥2-^ is given by [187] 

2 1 

^2"^ = 3 V^^{m^±,m^±) + - t/2°^(TO^o,m^o) « V^^{m^,m^) , (4.84) 

where the arguments represent the masses of the two exchanged pions, which are given in Table [2] and with 
m,r denoting the average pion masses defined by 

2 1 

fh^ = - m^± + -m^o = 138.0390 MeV. (4.85) 
For the isovector part V^^^ one obtains [1871 1188j 

^1 ^ f y'2'7r{n^7r±,m^±) for pp and nn, 

^'^ 12 V^^ {m^± , TO^o ) - KjI {m^i ,m^±) ~ {m^" , m^o ) for np, / = 1 ; 

and for np,I — 0, we have V27r — V2Tr{'fh-K,i^TT)- In conjunction with the work performed in Ref. [68.' , it was 
found that the CIB effects due to pion-mass splitting in the NLO 2PE potential are negligibly small in P and 
higher partial waves when also the pion-mass dependence of the polynomial terms is included (cf. discussion 
in Appendix A of Ref. |188j ). However, in the ^Sq state, this effect is non-negligible (causing AaciB ~ —0.5 
fm), but can be absorbed by the non-derivative CIB contact. Thus, from a procedural point of view, the 
CIB effect from pion-mass splitting in the NLO 2PE can simply be ignored (which is the procedure applied 
in Ref. [68]). 

Nucleon-mass splitting in intermediate states of 2PE diagrams creates CSB as discussed in detail in 
Refs. [1891 1190] . When derived just for the NLO 2PE diagrams, this effect is N'^LO. Accurate calculations 
conducted within the framework of conventional meson theory [179| have shown that this effect is very small 
in P and higher partial waves for diagrams with only nucleon intermediate states and moderate when single 
A excitation is taken into account (which corresponds to order N'^LO in A-less ChPT). The effect in ^Sq 
is non-negligible, but can be absorbed by the non-derivative CSB contact term. All existing models for the 
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Table 6: Columns three to five display the x'^ /datum for the reproduction of the 1999 np database |173| (subdivided into 
energy intervals) by various np potentials. For the chiral potentials, the x'^ /datum is stated in terms of ranges which result 
from a variation of the cutoff parameters used in the regulator functions. The values of these cutoff parameters in units of MeV 
are given in parentheses. Ti^b denotes the kinetic energy of the incident nucleon in the laboratory system. 



Tlab bin 


^ of np 


Idaho 


Juelich 


Argonne 


(MeV) 


data 


N^LO [SH] 


N^LO |171| 


^18 |174j 






(500-600) 


(600/700-450/500) 




0-100 


1058 


1.0-1.1 


1.0-1.1 


0.95 


100-190 


501 


1.1-1.2 


1.3-1.8 


1.10 


190-290 


843 


1.2-1.4 


2.8-20.0 


1.11 


0-290 


2402 


1.1-1.3 


1.7-7.9 


1.04 



NN interaction at N'^LO ignore this CSB. However, it would be worthwhile to calculate this effect accurately 
for chiral 2PE and compare to the comprehensive results obtained from conventional meson models |179| 
where it was found that the CSB splitting of the singlet scattering length can be explained by nucleon-mass 
difference alone. 

Nucleon masses also enter the outer legs of NN scattering diagrams and the LS equation. Here, practi- 
tioners use, in general, the accurate proton and neutron masses and the proper relativistic formula for the 
relation between the on-shell CM momentum and the kinetic energy of the incident nucleon in the laboratory 
system (cf. Appendix A. 3 of Ref. [13]). 

Irreducible ttj exchange |191) causes CIB of order N^L0 (^ e^/(47r)^). It is a moderate effect in ^5*0 
{AaciB ~ —0.35 fm) and small in P and higher partial waves (cf. Ref. [13]). Some N^LO potentials [68] 
do include this TT'y contribution. Corrections to the ttj graphs (which are N'*L0) and even 2ttj diagrams 
(N^L0) have recently been calculated by Kaiser |192H194] and found to be astonishingly large. No NN 
calculation has yet included them. 

Finally, at N^L0, several corrections to the long-range electromagnetic interaction occur. The leading 
relativistic correction to the static Coulomb potential |176[ I195j is most conveniently included by replacing 
the fine-structure constant a by 

Other electromagnetic contributions to pp scattering are two-photon exchange, the Darwin-Foldy term, 
vacuum polarization, and the magnetic moment interaction |139l I174| . In the case of np scattering, the 
electromagnetic interaction consists only of the magnetic moment contribution. The electromagnetic inter- 
actions are important contributions to the scattering amplitude. Note, however, that in the calculation of 
the strong nuclear phase shifts the electromagnetic interaction is only of relevance when its distortion of the 
wave functions affects the nuclear phase shifts in a non-negligible way. It is well known that this effect is 
large for the Coulomb potential in essentially all partial waves, but it is negligible for all other electromag- 



Table 7: Same as Table [6| but for pp. 
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Figure 14: Neutron-proton phase parameters as described by various chiral potentials at N^LO. The (red) solid and the dashed 
curves are calculated from Idaho N^LO potentials 1681 with A = 500 and 600 MeV, respectively; while the (green) dash-dotted 
and the dotted curves are based upon Juelich N^LO potentials 11711 with cutoff combinations 600/700 and 450/500 MeV, 
respectively. Partial waves with total angular momentum J < 1 are displayed. Empirical phase shifts (solid dots and open 
circles) as in Fig. [s] 



netic interactions in all partial- waves, except ^Sq below 30 MeV, where the effect can be calculated fairly 
model-independently and has been tabulated in Ref. |176j . Thus, in practice, it is sufficient to calculate pp 
phase shifts with only the Coulomb effect and np phase shifts without any electromagnetic effects (which is 



the way the phase shifts published in Refs. [TU [T3] and shown in the tables of Appendix F are calculated) 



4-6.3. A quantitative NN potential at JSPLO 

After previous sections have thoroughly prepared the terrain, we will now present a quantitative NN 
potential at N-^LO, including details of construction and results. We choose the Idaho N^LO pot entia l [55] 
as example. The isospin symmetric part of the chiral NN potential at N^LO is defined in Eq. (4.57) and 
the isospin violating terms were discussed in the previous section. Numerous parameters are involved which 
can be subdivided into three groups: ttN LECs, NN contact parameters, and the cutoff parameter A of the 
regulator Eq. (4.63). For A we choose initially 500 MeV. Within a certain reasonable range, results should 
not depend sensitively on A (cf. discussion in Section 4.5). Therefore, we have also made a second fit for 
A = 600 MeV. 



Data fitting and results for NN scattering. The fitting procedure starts with the peripheral partial waves 
because they involve fewer and more fundamental parameters. Partial waves with L > 3 are exclusively 
determined by IPE and 2PE because the N^LO contacts contribute to L < 2 only IPE and 2PE at N^LO 
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Figure 15: Same as Fig. [14] but J = 2 phase shifts and J < 2 mixing parameters are shown. 



depend on the axial-vector coupling constant, qa (we use gA — 1-29), the pion decay constant, /tt — 92.4 
MeV, and eight low-energy constants (LECs) that appear in the dimension-two and dimension-three TriV 
Lagrangians, Eqs. (2.571 and (2.58). The LECs are listed in Table |2j where column 'iViV Potential' shows 
the values used for the present N'^LO potential. In the fitting process, we varied three of them, namely, C2, 
C3, and C4. We found that the other LECs are not very effective in the NN system and, therefore, we left 
them at their central values as determined in ttTV analysis. The most infiuential constant is C3, which — in 
terms of magnitude — has to be chosen on the low side (slightly more than one standard deviation below 
its TiN determination), otherwise there is too much central attraction. Concerning C4, our choice C4 — 5.4 
GeV^^ lowers the phase shift (and slightly the ^^3) bringing it into closer agreement with the phase 
shift analysis — as compared to using the nN value C4 — 3.4 GeV^^. The other F waves and the higher 
partial waves are essentially unaffected by this variation of C4. Finally, the change of C2 from its TriV value 
of 3.28 GeV~^ to 2.80 GeV~^ (our choice) brings about some subtle improvements of the fit, but it is not 
essential. Overall, the fit of all J > 3 waves is very good. The _F-wave phase shifts are, in fact, described 
better than in the perturbative calculation shown in Fig. [8] because the regulator moderates the attractive 
surplus, thus, simulating correctly higher order contributions beyond the present order. 

We turn now to the lower partial waves. Here, the most important fit parameters are the ones associated 
with the 24 contact terms that contribute to the partial waves with L < 2 (cf. Section 4.3 and Table |4]). 
In addition, we have two charge-dependent contacts which are used to fit the three different ^5*0 scattering 



lengths. Op 



and a„ 



In the optimization procedure, we fit first phase shifts, and then we refine the fit by minimizing the 
obtained from a direct comparison with the data. We start with pp, since the pp phase shifts and data are 
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Table 8: Scattering lengths (a) and effective ranges (r) in units of fm. (a^ and refer to the pp parameters in the presence 
of the Coulomb force, and denote parameters determined from the nuclear force only and with all electromagnetic 
effects omitted.) 



Idaho 


Juelich 








N^LO p] 


N^LO iTT] 


CD-Bonii[T3] 


AV18'174| 


Empirical 


(500) 


(550/600)^ 









pp 



1^ 



nn 
N 



at 

n 



-7.8188 
2.795 

-17.083 
2.876 

-18.900 
2.838 

-23.732 
2.725 

5.417 
1.752 



-7.8003 
2.737 

-16.423 
2.828 

-18.900 
2.770 

-23.613 
2.651 

5.417 
1.742 



'So 

-7.8154 

2.773 

-17.460 
2.845 

-18.968 
2.819 

-23.738 
2.671 

5.420 
1.751 



-7.8138 
2.787 

-17.164 
2.865 

-18.818 
2.834 

-23.732 
2.697 

5.419 
1.753 



-7.8196(26) [176] 
-7.8149(29) [196] 
2.790(14) [nSj 
2.769(14) [inS] 



-18.95(40) [T8n[T82] 
2.75(11) [178] 
-23.740(20) [13] 
[2.77(5)] m 

5.419(7) [ig 
1.753(8) [I3] 



more accurate than the np ones. The pp fit fixes essentially the 7=1 potential. The I = I np potential is 
just the pp one modified by charge-dependence due to nucleon-mass difference, pion-mass splitting in IPE, 
TTj exchange, and omission of Coulomb (as discussed in Section 4.6.2 ). In addition to this, the non-derivative 
contact in the ^5*0 state is changed such as to reproduce the np scattering length. The nn potential is the 
pp one without Coulomb, using neutron masses, and fitting the nn scattering length in the state with 
the non-derivative contact. 

The /datum for the fit of the np data below 290 MeV are shown in Table [6j and the corresponding ones 
for pp are given in Table [Tj These tables reveal that at N'^LO a x^/datum comparable to the high-precision 
Argonne Vig [174] potential can, indeed, be achieved. The Idaho N'^LO potential [S5| with A ~ 500 MeV 
produces a /datum =1.1 for the world np data below 290 MeV which compares well with the x'^ /datum = 
1.04 by the Argonne potential. In 2005, also the Juelich group produced several N^LO NN potentials [171] . 
the best of which fits the np data with a x^/datum — 1.7 and the worse with 7.9 (Table[6]). 



Table 9: Deuteron properties as predicted by various NN potentials are compared to empirical information. (Deuteron binding 
energy B^, asymptotic S state As, asymptotic D/S state rj, deuteron radius r^, quadrupole moment Q, D-state probability 
Po; the calculated r^j and Q are without meson-exchange current contributions and relativistic corrections.) 





Idaho 
N^LO [68] 
(500) 


Juelich 
N^LO [m] 
(550/600) 


CD-Bonn [13] 


AV18[T74] 


Empiricar 


Bd (MeV) 


2.224575 


2.218279 


2.224575 


2.224575 


2.224575(9) 


As (fm-i/2) 


0.8843 


0.8820 


0.8846 


0.8850 


0.8846(9) 


■n 


0.0256 


0.0254 


0.0256 


0.0250 


0.0256(4) 


Td (fm) 


1.975 


1.977 


1.966 


1.967 


1.97535(85) 


Q (fm^) 


0.275 


0.266 


0.270 


0.270 


0.2859(3) 


Pd (%) 


4.51 


3.28 


4.85 


5.76 





"See Table XVIII of Ref. [13| for references; the empirical value for is from Ref. [197) . 
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Figure 16: Deuteron wave functions: the family of larger curves are 5-waves, the smaller ones D-waves. The thick colored 
lines represent the wave functions derived from chiral NN potentials at order N^LO (red solid: Idaho(500) [68^, green dashed: 
Juelich(550/600) |171l ). The thin dashed, dash-dotted, and dotted lines refer to the wave functions of the CD-Bonn 
Nijm-I [12], and AV18 |174| potentials, respectively. 



Turning to pp, the for PP data are typically larger than for np because of the higher precision of pp 
data. Thus, the Argonne Vis produces a x^/datum — 1.4 for the world pp data below 290 MeV and the best 
Idaho N^LO pp potential obtains 1.5. The fit by the best Juelich N^^LO pp potential results in a x^/datum 
= 2.9 and the worst produces 22.3. 

Phase shifts of np scattering from two Idaho (solid and dashed lines) and two Juelich (dash-dotted 



and dotted lines) N'^LO np potentials are shown in Figs. 14 and 15 The phase shifts confirm what the 
corresponding have already revealed. Low-energy scattering parameters are listed in Table [sj A few 
more technical details and tables with numerical values for contact and phase parameters are provided in 
[Appendix F[ 

The deuteron. The reproduction of the deuteron parameters is shown in Table [9j We present results for two 
N^LO potentials, namely, Idaho [BS] with A = 500 MeV and Juelich [171j with cutoff combination 550/600 
MeV. Remarkable are the predictions by the chiral potentials for the deuteron radius which are in good 
agreement with the latest empirical value obtained by the isotope-shift method [ 197] . All TV TV potentials 
of the past (Table [9] includes two representative examples, namely, CD-Bonn [13 and AV18 [174] ') fail to 
reproduce this very precise new value for the deuteron radius. 



In Fig. 16 we display the deuteron wave functions derived from the N^LO potentials and compare them 
with wave functions based upon conventional A^A^ potentials from the recent past. Characteristic differences 
are noticeable; in particular, the chiral wave functions are shifted towards larger r due to softer cutoffs and 
the effect of the contact terms, which explains the larger deuteron radius. 



5. Nuclear many-body forces 

The chiral 2NF discussed in the previous section has been applied in microscopic calculations of nuclear 
structure with, in general, a great deal of success |75H88j . However, from high-precision studies conducted in 
the 1990s, it is well known that certain few-nucleon reactions and nuclear structure issues require 3NFs for 
their microscopic explanation. Outstanding examples are the Ay puzzle of N-d scattering |1981 1199j and the 
ground state of j200j . As noted before, an important advantage of the EFT approach to nuclear forces 
is that it creates two- and many-nucleon forces on an equal footing (cf. the overview given in Fig. [T]). In this 
section, we will explain in some detail chiral three- and four-nucleon forces. We will limit our presentation 
to the isospin-symmetric case; isospin- violating 3NFs are discussed in Ref. [201) . 



50 



(a) (b) (c) 

Figure 17: Thrce-nucleon force diagrams at NLO. Notation as in Fig.[l] 



5.1. Three-nucleon forces 

Nuclear three-body forces in ChPT were initiaUy discussed by Weinberg 
first derived by van Kolck [70] . 

For a 3NF, we have A = 3 and C = 1 and, thus, Eq. (|3.4|) imphes 



and the 3NF at NNLO was 



(5.1) 



1/ = 2 + 2L + ^ Ai . 

i 

We win use this equation to analyze 3NF contributions order by order. 

5.1.1. Next-to-leading order 

The lowest possible power is obviously 1^ — 2 (NLO), which is obtained for no loops {L = 0) and only 
leading vertices {J2i'^i — We display typical graphs in Fig. 17 As discussed by Weinberg [52], the 
contributions from these diagrams vanish at NLO. To see this, let's first look at graph (a), which contains a 
Weinberg- Tomozawa vertex, Eq. (A.43I, that includes a time-derivative of a pion field. Since this diagram 



does not involve reducible topologies, it can be treated as a Feynman diagram in which energy is conserved 
at each vertex, so that the time-derivative yields a difference of nucleon kinetic energies ~ Q^/Mjv instead 
of ~ Q. Thus, the contribution from this graph is suppressed by a factor Q/M^ and demoted to NNLO. 



Graphs (b) and (c) of Fig. 17 are best discussed in terms of time-ordered perturbation theory. Weinberg [52] 



and van Kolck [70] showed that, at NLO, the irreducible topologies of these graphs cancel against the recoil 
corrections from the reducible ones, leaving no net irreducible 3N contribution. What remains is just the 
iteration of the static 2N potentials. In fact, this had been pointed out already by Yang and Glockle |127] 
and Coon and Friar [US] in the 1980's. 

The bottom line is that there is no genuine 3NF contribution at NLO. The first non-vanishing 3NF 
appears at NNLO. 

5.1.2. Next-to-next-to-leading order 

The power v = 3 (NNLO) is obtained when there are no loops [L — 0) and Aj = 1, i.e., A^ = 1 for 
one vertex while A^ = for all other vertices. There are three topologies which fulfill this condition, known 



as the two-pion exchange (2PE), IPE, and contact graphs [70101 (Fig. 18). 

Using the subleading vertices Eqs. (A. 51 1 and (A.52), it is straightforward to derive the 2PE 3N-potential 
to be 

9 A \ ^1 >^ qi){Sj ■ Qj) 

2fJ 2 



T/3NF _ 
^2PE ~ 



E 



{qf + ml){q^+ml) 



ijk 



(5.2) 



with Qi = Pi' — Pi, where pi and p/ are the initial and final momenta of nucleon i, respectively, and 



IciTO^ 2C3 
f2 f2 

J IT J IT 



C4 

J TT 



^abc 



(5.3) 



There are great similarities between this force and earlier derivations of 2PE 3NFs, notably the 50-year 
old Fujita-Miyazawa [202^, the Tucson-Melbourne (TM) [203] . and the Brazil [204] forces. A thorough 
comparison between various 2PE 3NFs is conducted in Ref. [205] resulting in the recommendation to drop 
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Figure 18: The three-nucleon force at NNLO. From left to right: 2PE, IPE, and contact diagrams. Notation as in Fig. [T] 



the so-called "c-term" from the TM force, since it does not have an equivalent in the ChPT derived force, 



Eqs. ([5^, ([5^, giving rise to the construction of the TM' (or TM99) force [206] , 

Notice that Eq. (5.3 1 does not include a C2-term. Due to two time-derivatives, the contribution from the 
C2 vertex is (Q/M^)'^ suppressed and demoted by two orders. Note also that the 2PE 3NF does not contain 
any new parameters, because the LECs ci, C3, and C4 appear already in the 2PE 2NF (Section 4.1.2) and 
are fixed by ttN and/or NN data. 



The other two 3NF contributions shown in Fig. 18 are easily derived by taking the last two terms of the 
A = 1 Lagrangian, Eq. (2.66), into account. The IPE contribution is 



y3NF _ _ n 



9A 



E 



^3 ■ ^3 



(r,; • Tj){a, ■ qj) 



and the 3N contact potential reads 



3NF 



2 ^ 



Tk ■ 



(5.4) 



(5.5) 



These 3NF terms involve the two new parameters D and E, which do not appear in the 2N problem. There 
are many ways to pin these two parameters down. In Ref. |71j , the triton binding energy and the nd doublet 
scattering length ^a„(j were used. One may also choose the binding energies of '^H and ^He [51] or an optimal 
over-all fit of the properties of light nuclei |:99: . Exploiting the consistency of interactions and currents in 
ChPT |207j , the parameter D of the ttNNNN vertex involved in the IPE 3NF can be constrained by p-wave 
pion-production data |208j or electro-weak processes like the tritium /3-decay |209| or proton-proton fusion 
{pp de'^ Vf.) |210j . Once D and E are fixed, the results for other 3N, 4N, etc. observables are predictions. 

The 3NF at NNLO has been applied in calculations of few-nucleon reactions [TU UHl-EZj : structure of 
light- and medium-mass nuclei [98 -^ ICFIl , and nuclear and neutron matter |102[ I103j with a good deal of 
success. Yet, the famous 'Ay puzzle' of nucleon-deuteron scattering is not resolved [7T1[94]. When only 
2NFs are applied, the analyzing power in p-'^He scattering is even more underpredicted than in p-d [7^1211] . 
However, when the NNLO 3NF is added, the p-'^He Ay substantially improves (more than in p-d) ^\ — but 
a discrepancy remains. Furthermore, the spectra of light nuclei leave room for improvement |99| . 

We note that there are further 3NF contributions at NNLO, namely, the \/Mm corrections of the NLO 
3NF diagrams (Fig. 17 1. Some of these terms involve the vertices Eqs. (A. 49 1 and (A.50), and the 1/Mm 
correction of the C4 vertex, Eq. (A.52); others are due to higher order recoil corrections. Several of those 



contributions have been calculated by Coon and Friar in 1986 |128| . These corrections are believed to be 
very small. 

To summarize, the 3NF at NNLO is a remarkable contribution: It represents the leading many-body 
force within the scheme of ChPT; it includes terms that were advocated already some 50 years ago; and it 
produces noticeable improvements in few-nucleon reactions and the structure of light nuclei. But unresolved 
problems remain. Moreover, in the case of the 2NF, we have seen that one has to proceed to N'^LO to 
achieve sufficient accuracy. Therefore, the 3NF at N^LO is needed for at least two reasons: for consistency 
with the 2NF and to hopefully resolve outstanding problems in microscopic structure and reactions. 



5.1.3. Next-to-next-to-next-to-leading order 

According to Eq. (5.1), the value v = A, which corresponds to N^LO, is obtained for the following classes 
of diagrams. 
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(a) 



(b) 



(c) 



(d) 



(e) 



Figure 19: Leading one-loop 3NF diagrams at N^LO. We show one representative example for each of five topologies, which 
are; (a) 2PE, (b) 1PE-2PE, (c) ring, (d) contact-lPE, (e) contact-2PE. Notation as in Fig.[l] 



Figure 20: Sub-leading one-loop 3NF diagrams which appear at N^LO. Notation as in Fig.^ 



3NF loop diagrams at N^LO. For this group of graphs, we have L = 1 and, therefore, all have to be zero 
to ensure v = 4. Thus, these one-loop 3NF diagrams can include only leading vertices, the parameters of 
which are fixed from TriV and NN analysis. We show five representative examples of this very large class of 
One sub-group of these diagrams (2PE graphs) has been calculated by Ishikawa and 
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diagrams in Fig. 

Robilotta [72], and two other topologies (1PE-2PE and ring diagrams) have been evaluated by the Bonn- 
Jiilich group (73]. The remaining topologies, which involve a leading four-nucleon contact term [diagrams 
(d) and (e) of Fig. 19 , are under construction by the Bonn-Jiilich group. The N'^LO 2PE 3NF has been 



applied in the calculation of nucleon-deuteron observables in Ref. [72] producing very small effects. 

The smallness of the 2PE loop 3NF at N^LO is not unexpected. It is consistent with experience with 
corresponding 2NF diagrams: the NLO 2PE contribution to the NN potential, which involves one loop and 
only leading vertices (Fig. [s]), is also relatively small (Fig. [s]). 

By the same token, one may expect that also all the other N^^LO 3NF loop topologies will produce only 
small effects. 



3NF tree diagrams at N^LO. The order v — 4: \s also obtained for the combination L = Q (no loops) and 
Aj = 2. Thus, either two vertices have to carry A; = 1 or one vertex has to be of the A^ = 2 kind, while 
all other vertices are A^ = 0. This is achieved if in the NNLO 3NF graphs of Fig. [18] the power of one vertex 
is raised by one. The latter happens if a relativistic I/AIn correction is applied. A closer inspection reveals 
that all 1/Mm corrections of the NNLO 3NF vanish and the first non- vanishing corrections are proportional 
to l/M^ and appear at N'^LO. However, there are non- vanishing 1/Af^ corrections of the NLO 3NF and 
there are so-called drift corrections |212) which contribute at N'^LO (some drift corrections are claimed to 
contribute even at NLO |212) ). We do not expect these contributions to be sizable. Moreover, there are 
contributions from the A^ = 2 Lagrangian [123J proportional to the low-energy constants di. As it turns 
out, these terms have at least one time-derivative, which causes them to be Q /Mjsi suppressed and demoted 
to N^LO. 

Thus, besides some minor 1/A/^ corrections, there are no tree contributions to the 3NF at N'^LO. 

Summarizing the 3NF at N^LO: For the reasons discussed, we anticipate that this 3NF is weak and will 
not solve any of the outstanding problems. In view of this expectation, we have to look for more sizable 
3NF contributions elsewhere. 

5.14. The 3NF at N^LO 

The obvious step to take is to proceed to the next order, N'^LO or ly — 5. Some of the tree diagrams 
that appear at this order were mentioned already: the 1/A/^ corrections of the NNLO 3NF and the trees 
with one di vertex which are I/Mn suppressed. Because of the suppression factors, we do not expect sizable 
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Figure 21: Disconnected four-nucleon force diagrams at NLO, which cancel against the recoil corrections from corresponding 
iterative diagrams. 



effects from these graphs. Moreover, there are also tree diagrams with one vertex from the = 3 t:N- 
Lagrangian [1221 1213] proportional to the LECs ei . Because of the high dimension of these vertices and 
assuming reasonable convergence, we do not anticipate much from these trees either. 

However, we believe that the loop contributions that occur at this order are truly important. They are 
obtained by replacing in the N'^LO loops (Fig. 19 1 one vertex by a = 1 vertex. We show five examples of 
this large group of diagrams in Fig. [20] This 3NF is presumably large and, thus, what we are looking for. 

The reasons, why these graphs are large, can be argued as follows. Corresponding 2NF diagrams are the 
three-pion exchange (3PE) contributions to the NN interaction. In analogy to Figs. 19 and 20 there are 
3PE 2NF diagrams with only leading vertices (Fig. [7| and the ones with one (sub-leading) Ci vertex (and the 
rest leading). These diagrams have been evaluated by Kaiser in Refs. |57l[58] and (5^, respectively. Kaiser 
finds that the 3PE contributions with one sub-leading vertex are about an order magnitude larger than the 
leading ones. 



5.1.5. 3NF summary 

To make a complicated story short, this is the bottom line concerning 3NF |162j : 

• The leading chiral 3NF (that appears at NNLO) is sizable, improves predictions, but also leaves 
unresolved problems. Therefore, additional sizable 3NF contributions are needed. 

• The chiral 3NF at N'^LO involves only leading vertices and most likely does not produce sizable 
contributions. 

• Sizable contributions are expected from the subleading one-loop 3NF diagrams that occur at N'^LO. 
These 3NF contributions may turn out to be the missing pieces in the 3NF puzzle and have the potential 
to solve the outstanding problems in microscopic nuclear structure. 



5.2. Four-nucleon forces 

Formally, the lowest order four-nucleon force (4NF) occurs for two separately interacting nucleon pairs 
with leading vertices and no loops. Fig. 21 This scenario is characterized by A = 4, C = 2, L = 0, and 



Ai — and therefore, according to Eq. (3.4), has the power v — 2 (NLO). However, similar to the 



NLO 3NF [Fig. [17] (b) and (c)], the 4NF diagrams of Fig. 21 cancel against the recoil corrections from 
corresponding iterative diagrams [70.. The disconnected 4NF diagrams of order four, which are obtained 
for either i = 1 or A^ = 2, also cancel |214j . Thus, we are left with just the connected (C = 1) A = 4 
diagrams for which Eq. ( ]3.4[ ) yields 

z^ = 4 + 2L + ^Ai. (5.6) 



Therefore, a connected 4NF appears for the first time aX v — A (N'^LO), with no loops and only leading 
vertices. Fig. 22 This 4NF includes no new parameters and does not vanish [2141 1215j . Some graphs in 
Fig. 22 appear to be reducible (iterative). Note, however, that these are Feynman diagrams, which are 



best analyzed in terms of time-ordered perturbation theory. The various time-orderings include also some 
irreducible topologies (which are, by definition, 4NFs). Or, in other words, the Feynman diagram minus 
the reducible part of it yields the (irreducible) contribution to the 4NF. The reducible part depends on the 
two-, three-, and four-body scattering equations used. 
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Figure 22: Leading four-nucleon force at N'^LO. 

Assuming a good rate of convergence, a contribution of order (Q/A^)^ is expected to be rather small. 
Thus, ChPT predicts 4NF to be essentially insignificant, consistent with experience. Still, nothing is fully 



proven in physics unless we have performed explicit calculations. Recently, the leading 4NF (Fig. 22 ) has 
been applied in a calculation of the "^He binding energy, where it contributes a few 100 keV '216'. It should 
be noted that this preliminary calculation involves many approximations, but it certainly provides the right 
order of magnitude of the result, which is indeed very small as compared to the full "^He binding energy of 
28.3 MeV. 



6. Introducing A-isobar degrees of freedom 

The lowest excited state of the nucleon is the A(1232) resonance or isobar (a tt-N P-wave resonance 
with both spin and isospin 3/2) with an excitation energy of AM = A/a — Afjv = 293 MeV. Because of its 
strong coupling to the 'k-N system and low excitation energy, it is an important ingredient for models of 
pion-nucleon scattering in the A-region and pion production from the two-nucleon system at intermediate 
energies, where the particle production proceeds prevailingly through the formation of A isobars [TT]. At low 
energies, the more sophisticated conventional models for the 27r-exchange contribution to the NN interaction 
include the virtual excitation of A's, which in these models accounts for about 50% of the intermediate-range 
attraction of the nuclear force — as demonstrated by the Bonn potential [20l I217j . 

Because of its relatively small excitation energy, it is not clear from the outset if, in an EFT, the A should 
be taken into account explicitly or integrated out as a "heavy" degree of freedom. If it is included, then 
AAf ~ rriTT is considered as another small expansion parameter, besides the pion mass and small external 
momenta. This scheme has become known as the small scale expansion (SSE) |218j . Note, however, that 
this extension is of phenomenological character, since AAf does not vanish in the chiral limit. 

In the chiral EFT discussed so far in this report (also known as the "A- less" theory), the effects due to 
A isobars are taken into account implicitly. Note that the dimension- two LECs, the c^, have unnaturally 
large values (cf. Table [2]). The reason for this is that the A-isobar (and some meson resonances) contribute 
considerably to the Ci — a mechanism that has become known as resonance saturation [2X91 . Therefore, the 
explicit inclusion of the A ("A-fuU" theory) will take strength out of these LECs and move this strength 
to a lower order [S31 ISSl I220H222| . As a consequence, the convergence of the expansion improves, which is 
another motivation for introducing explicit A-degrees of freedom. We observed that, in the A-less theory, 
the subleading 2PE and 3PE contributions to the 2NF are larger than the leading ones. The promotion of 
large contributions by one order in the A-fuU theory fixes this problem. 

In the heavy baryon formalism, the leading Lagrangian involving A's reads [Ml ED] (listing only terms 
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Chiral 2N Force 



A- less 



Additional in A-full 



LO 



NLO 



NNLO 



N^LO 



N^LO 




Figure 23: Chiral 2NF without and with A-isobar degrees of freedom. Arrows indicate the shift of strength when explicit A's 
are added to the theory. Note that the A-full theory consists of the diagrams involving A's plus the A-less ones. Double lines 
represent A-isobars; remaining notation as in Fig. [l] 

relevant to our present discussion) 



A,=0 



A{ido - AM) A - ^ (nTSA 



h-cl -^t^-Dt NraN ■ ( NTS A + h.c. 



(6.1) 



where A is a four-component spinor in both spin and isospin space representing the A-isobar and Ha and 
Dt are LECsQ Moreover, are 2x4 spin transition matrices which satisfy S^S^^ = (25*^ — ie^^'^a'^) /i 
and T° are similar isospin matrices with T°'T^^ = {26 



ah 



abc^c 



r'^)/3. Notice that, due to the heavy baryon 
expansion, the mass of the A-isobar, Ma, has disappeared and only the small mass difference AM enters. 

The LECs of the irN Lagrangian are usually extracted in the analysis of tt-N scattering data and clearly 
come out differently in the A-full theory as compared to the A-less one. While in the A-less theory, the 
magnitude of the LECs C3 and C4 is about 3-5 GeV~^ (cf. Table [2|, they turn out to be around 1 GeV~^ in 
the A-full theory 

In the 2NF, the virtual excitation of A-isobars requires at least one loop and, thus, the contribution occurs 
first atu = 2 (NLO), see Fig.[23| The A contributions to the 2PE were first evaluated in Refs. [551IM1E^ 



using time-ordered perturbation theory and later by Kaiser et al. [SS] in covariant perturbation theory. 



5 Our convention for Ha is consistent with Refs. [5il l56l [TOllTOTl and differs by a factor of two from Refs. [2T8l [2211 [223] . 
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Chiral 3N Force 



A-less Additional in A-full 

LO 



NLO 



NNLO 




N^LO 



N^LO 




Figure 24: The 3NF without and with A-isobar degrees of freedom. Arrows indicate the shift of strength when explicit A's 
are added to the theory. Note that the A-fuU theory consists of the diagrams involving A's plus the A-Iess ones. Double lines 
represent A-isobars; remaining notation as in Fig. [T] 



Recently, also the NNLO contributions have been worked out |221j . Krebs et al. |221j verified the consistency 
between the A-full and A-less theories by showing that the contributions due to intermediate A-excitations, 
expanded in powers of 1/AM, can be absorbed into a redefinition of the LECs of the A-less theory. The 
corresponding shift of the LECs C3, C4 is given by 



C3 



-2c4 



9AM 



(6.2) 



Using Ha = Sqa/V^ (large Nc value), almost all of C3 and an appreciable part of C4 is explained by the A 
resonance. 

The studies of Refs. [5511221] confirm that a large amount of the intermediate-range attraction of the 2NF 
is shifted from NNLO to NLO with the explicit introduction of the A-isobar. However, it is also found that 
the NNLO 2PE potential of the A-less theory provides a very good approximation to the NNLO potential 
in the A-full theory. 

The A isobar also changes the 3NF scenario, see Fig. [24j The leading 2PE 3NF is promoted to NLO. 
In the A-full theory, this term has the same mathematical form as the corresponding term in the A-less 
theory, Eqs. (5.2) and (5.3), provided one chooses ci = and C3, C4 according to Eq. (6.2). Note that the 
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other two NLO 3NF terms involving A's vanish [222] as a consequence of the antisymmetrization of the 



3N states. The A contributions to the 3NF at NNLO |222j vanish at this order, because the subleading 
TVAtt vertex contains a time-derivative, which demotes the contributions by one order. However, substantial 
3NF contributions are expected at N'^LO from one-loop diagrams with one, two, or three intermediate A- 
excitations, which correspond to diagrams of order N'^LO, N^LO, and N'^LO, respectively, in the A-less 
theory. 3NF loop-diagrams with one and two A's are included in the Illinois force [224] in a simplified way. 

To summarize, the inclusion of explicit A degrees of freedom does certainly improve the convergence of 
the chiral expansion by shifting sizable contributions from NNLO to NLO. On the other hand, at NNLO 
the results for the A-full and A-less theory are essentially the same. Note that the A-fuU theory consists 
of the diagrams involving A's plus all diagrams of the A-less theory. Thus, the A-full theory is much more 
involved. Moreover, in the A-full theory, l/M^ 2NF corrections appear at NNLO (not shown in Fig. 23), 
which were found to be uncomfortably large by Kaiser et al. [SS]. Thus, it appears that up to NNLO, the 
A-less theory is more manageable. 

The situation could, however, change at N'^LO where potentially large contributions enter the picture. 
It may be more efficient to calculate these terms in the A-full theory, because in the A-less theory they are 
spread out over N^LO, N'^LO and, in part, N^LO. These higher order contributions are a crucial test for 
the convergence of the chiral expansion of nuclear forces and represent a challenging topic for the future. 



7. Conclusions 

The past 15 years have seen great progress in our understanding of nuclear forces in terms of low-energy 
QCD. Key to this development was the realization that low-energy QCD is equivalent to an effective field 
theory (EFT) which allows for a perturbative expansion that has become know as chiral perturbation theory 
(ChPT). In this framework, two- and many-body forces emerge on an equal footing and the empirical fact 
that nuclear many-body forces are substantially weaker than the two-nucleon force is explained naturally. 

In this review, we have shown in detail how the two-nucleon force is derived from ChPT and demonstrated 
that, at N'^LO, the accuracy can be achieved that is necessary and sufficient for reliable microscopic nuclear 
structure predictions. First calculations applying the N'^LO NN potential 68 in the conventional shell 
model [75H77] , the ab initio no-core shell model [78l |79] , the coupled cluster formalism [80H86] , and the 
unitary-model-operator approach [87l 188] have produced promising results. 

We also discussed nuclear many-body forces based upon chiral EFT. The 3NF at NNLO has been known 
for a while [701 EI] and applied in few-nucleon reactions [TT] [5M77] , structure of light- and medium-mass 
nuclei [55^101] . and nuclear and neutron matter |1021 1103"] with some success. However, the famous ^Ay 
puzzle' of nucleon-deuteron scattering is not resolved by the 3NF at NNLO. Thus, one important open issue 
are the few-nucleon forces beyond NNLO ("sub-leading few-nucleon forces") which, besides the Ay puzzle, 
may also resolve some important outstanding nuclear structure problems. As explained, this may require 
going even beyond N^LO. 

Another open question is the convergence of the chiral expansion (of the two- as well as the three-nucleon 
potentials) at orders beyond N'^LO, for which the inclusion of A-isobar degrees of freedom may be useful. 
Furthermore, the non-perturbative renormalizations of the chiral potentials require more work and a better 
understanding. 

Finally, we note that topics of interest we did not discuss include parity-violating nuclear forces and 
consistent electroweak currents. 

Having identified some of the open issues, we hope that this review will be helpful towards future progress. 

If the outstanding problems are resolved within the next few years, then, after 80 years of desperate 
struggle, we may finally claim that the nuclear force problem is essentially under control. The greatest 
beneficiary of such progress will be the field of ab initio nuclear structure physics. 
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Appendix A. Notation, conventions, and Feynman rules 



Appendix A.l. Notation and conventions 

The contravariant space-time four-vector is given by 



and the four-momentum vector reads 



x^" = {t, x) 
jf = {E,p). 



(A.1) 
(A.2) 



We use imits such that h = c=l. 

Greek indices ijl,v, etc. run over the four space-time coordinate labels 0,1,2,3, with a;° = i the time 
coordinate. Latin indices i,j,k, and so on run over the three space coordinate labels 1,2,3. The metric is 
diagonal with 

/I \ 



at' 



-1 



-1 



V 



-1 / 



and the covariant versions of the above-mentioned vectors are 



= g^xux" = {t, -x) , Pfj, = gf.i.p" = {E, -p) , 
where summation over repeated indices is always understood; also 



(A.3) 



(A.4) 



(A.5) 



While for an ordinary three- vector we have, in general, x = {x^,x'^,x^), there is caution in place with the 
(three-dimensional) nabla operator which is defined to be 



V = (Vi,V2,V3) 

The four-momentum operator reads 



dx' 



d_ 
dxi 



Pf^ 



. d 

dx^ 
. d 



id^ = {ido,iV^ = 



.0 : 



d 

'dt^-^ 



The relativistic nucleon field satisfies the free Dirac equation 

{i) - Mn) ^{x) = (7"^^ - Mn) ^{x) = 



(A.6) 

(A.7) 
(A.8) 

(A.9) 



where Mj^ denotes the nucleon mass and 7^ the Dirac matrices which we apply in Dirac-Pauli representation 



/ 
-/ 



r = 



a' 
-cr* 



with I the two-dimensional identity matrix and cr' the Pauli matrices 



a' = 



1 

1 



-i 

1 
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1 
-1 



(A.10) 



(A.11) 



[a\a^] = ctV - aV' = 2^e'•'■*=<T^ {a\a^} = a*a^ + a^a' = 25^^ , (A.12) 

ctV = ie'^'^a'' + 5'^ , (A.13) 

with 

{+1 if {i, j, k) even permutation of (1, 2, 3) 
-1 if odd permutation (A. 14) 

otherwise 

Following convention, we denote the Pauli matrices by t" (a = 1, 2, 3) when operating in isospin space, with 

r V = ie"''^T^ + (A.15) 

Notice that, for isospin components, it does not make sense to distinguish between upper and lower indices 
and, therefore, subscripts and superscripts have the same meaning. 
The Dirac matrices have the properties 

17 ) 7 I = ^9 , 7 =777=7) 7 =777= -7. (A.lbj 



The 75-matrix is defined by 



Commutator of 7 matrices: 



^5=^5 = ^/y^V=( J 0). (A.17) 
{7^7'} = 0, (7')' = 1, 7'^ =7'- (A.18) 

= I [7^ 7l , 7^7" = 9^"" - i^^'^" , (A.19) 



The relativistic Dirac field for positive-energy nucleons is 



*(^) = E/ ■ii^2\l^<P^')^te-'^-^b{p,s,t) (A.21) 
with the Dirac spinor given by 



^P» = V^^ i- X. (A.22) 

where Ep = p° = ^p^ + Mj^. The Pauli spinors Xs and describe, respectively, the spin and isospin of the 
nucleon. Further, h{p, s, t) and {p, s, t) are destruction and creation operators for a nucleon with momentum 
p, and spin and isospin quantum numbers s and t, respectively. They satisfy the anti-commutation relations 

{h{p,s,t),hHp',s',t')}=5,,,,5t^t'5\p-p'), (A.23) 

{h{p,s,t),b{p',s',t')] = {bHp,s,t),h\p',s',t')}=Q. (A.24) 
* = *^7" (A.25) 

In the heavy-baryon formalism, the free field equation for nucleons is, in leading order and using = 
(1,0,0,0), 

idoN{x) = (A.26) 



60 



and the nucleon field is 



Nix) ^Y.j 



dH 



(27r)3/2 



Xs£,te-''-^b{l,s,t) 



where Iq ~Q dX leading order and r /2Mjsi at NLO. 



(A.27) 



(A.28) 



The pion fields are in terms of their cartesian components (i = 1,2,3; no distinction between upper and 
lower index i) 

(fq 1 



■Ki{x) 



(27r)3/2 



-iq-x 



(A.29) 



where u = = y'' (p + and ai and a\ are, respectively, the destruction and creation operators obeying 
the commutation relations 

a,{q),a]{q') \ ^ S.jS^'iq- q') , 



[ai{q),aj{q")] = al{q),a]{q') 
The charged and neutral pion fields are given by 

1 . 



0. 



V2 

1 

7^3 ■ 



[TTl + 1112) 



(tti -1112) , 



Throughout this article, we state amplitudes in terms of the "potential" V which is defined by 

V = iM 



(A.30) 
(A.31) 

(A.32) 

(A.33) 
(A.34) 

(A.35) 



where Ai is the invariant amplitude, calculated according to Feynman rules. The relation of A4 to the 
S-matrix is 



+ {2n)U\p,+p2-p{~p^) 

X M{p{,P2]PuP2) ■ 



(27r)6 \E[E'^EiE2 



(A.36) 



Appendix A. 2. Feynman rules 

The basic formalism underlying the Feynman rules is the Dyson expansion of the S'-matrix [cf. e. g., 
Eq. (6.1.1) of Rcf. [225J]. Since we are applying covariant perturbation theory, we assume Hj = —Cj 
(where Hi denotes the interaction Hamiltonian density and Cj the interaction Lagrangian density) and 
use the usual covariant Feynman propagators. From a procedural point of view, this is acceptable for 
the derivative couplings considered here (see, however, Ref. [226j for exceptions). For a discussion of the 
differences between covariant perturbation theory and time-ordered perturbation theory (also known as "old 
fashioned" perturbation theory), see the evaluation of the football diagram. Appendix B.2.2 below. 

In all one-pion vertices given below, q denotes the four-momentum of an outgoing pion of isospin compo- 
nent a. In all two-pion vertices, qi denotes the four-momentum of an ingoing pion with isospin component a 
and q2 the four-momentum of an outgoing pion with isospin component h. Nucleon momenta are, in general, 
denoted by p and p' . 
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Appendix A. 2.1. Leading order 

The relativistic version of the leading order ttN Lagrangian has been given in Eq. (2.41 1. We consider 
the axial-vector (AV) and the Weinberg- Tomozawa (WT) couplings, which involve one and two pion fields, 
respectively. The relativistic AV interaction Lagrangian is given by 



Cav ^-^^ri5T^ -d^TT. 



(A.37) 



For simple Lagrangians like this one, the vertex is just i times the Lagrangian stripped off the fields which 
generates {q out) 



9a an 



In the heavy baryon formalism, the AV Lagrangian reads [cf. Eq. (2.53)] 



Cav = -^Nt-{S-^)i,N 



with vertex {q out) 



The relativistic WT coupling term is 

CwT 

implying the vertex (gi in, q2 out) 



9 A a- ^ 



-— *7''r • (tt X a^Tr)* 



1 

4il 



The Heavy Baryon version is 



with vertex {qi in, q2 out) 



1 



CwT = - 775" iV T • (tt X d^Tr) N 



4/^ 



e^'^T^l^ql + ql) 



The relativistic nucleon propagator reads [cf. Eq. (A. 9)] 



^-MN + ie p2 _ + ie 



and the leading order heavy baryon version of the nucleon propagator is given by [cf. Eq. (A.26)] 

i 



The pion propagator is [cf. Eq. (2.32)] 



10 



(A.38) 

(A.39) 
(A.40) 

(A.41) 

(A.42) 
(A.43) 
(A.44) 

(A.45) 

(A.46) 
(A.47) 
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Appendix A.2.2. Next-to-leading order 

The "fixed" part of the dimension-two HB Lagrangian is given in Eq. (2.56). It leads to a nucleon kinetic 
energy correction 

^ 2 

(A.48) 



P 



2M 



N 



and produces the one-pion vertex {q out) 



gA 



and the two-pion vertex (qi in, q2 out) 



1 



'{p + p')-{qi + q2)- 



(A.49) 



(A.50) 



The vertices Eqs. (A.49) and (A.50) are the first relativistic corrections of the vertices Eqs. ( |A.40 ) and 



(A.44), respectively. These corrections can also be obtained by starting from the corresponding relativistic 
vertices Eqs. (A. 38) and (A.42), sandwiching them between Dirac spinors, and making a 1/Mn expansion. 

vo-pion 

(A.51) 



The "contact" part of the dimension- two Lagrangian, Eq. (2.57), gives rise to the following two-pion 
vertices (gi in, q2 out): 

10 



iS'''' 



P 



9\ 



4M, 



N 



9192 + 2c3gi^q^ 



and 



P 

J TT 



C4 



1 



(A.52) 



For more Feynman rules, see appendix A of Ref. |119j . 



Appendix B. Two-pion exchange contributions to the 2NF at NLO 

An overview of the 2PE diagrams that contribute at NLO was shown in Fig. [3] We will now evaluate 
the various groups of diagrams one by one. 

Appendix B.l. Triangle diagrams 

Appendix B.1.1. The triangles in covariant perturbation theory 

The triangle diagrams are shown in Fig. B.l The contributions from diagram (a) and (b) will be denoted 
by Va and Vb, respectively. According to relativistic Feynman rules, one obtains 

r,2 r ^4; 



16/4 J (2^)4 ^ " ' ^ ' ''"^ '"^ '''^'{l + q)^-ml + le 



l2^^2+,, 'M-p')(~nr^^.l. f^r:M^ :U l'' + l'')r2W2{-p-), (B.I) 



16/4 y (2^)4^ ' ^ ^ . ^ ' [l + qf - ml + le 

u,{p'){~r - g'^)rh.75 r?"" {I'Klpl^uiip) , (B.2) 



P-ml + ie n,^-...jY 



where the momenta are defined in Fig. B.l and q = (0,p ' — p). The sum over the isospin operators is 
^abc^c^b^a _ • Tj and, for the intermediate nucleon, we have "fulbi^N + MM)^p"fb = lv{lj^N — MN)-fp. 
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{Ep,p') {Ep,-p') 



I 



l + q 



kjsf 



l + q 



P <' 



(Ep,p) {E.,,-p) 



{Ep,p) {Ep^-p) 



Thus, 



kN= {E.p-f ,-p' -tj 



kN = [Ep -f,p- fj 



Km = [Jip - I, p 



(a) 



(b) 



Figure B.l: Two-pion exchange triangle diagrams at NLO. 



-ITI-T2 



g\ f dH 



(2/^ + g^)fIi(p')7M"i(p) 



{l + q)^-ml + ie 



1 

P — m'i + ie 



U2 



'k%-M% + ie 



1 



{I + qY — ml + ie 



P — m?. + ie 



kff-M^ + ie 



(B.3) 



(B.4) 



At NLO, only the lowest order in the 1/Mjv expansion is included. The nucleon mass dependence comes 
from the Dirac spinors and the nucleon propagator and performing such an expansion, we find at lowest 
order 



1 11 



kl - M% + ie 



2Mn -l" + ie 



(B.5) 
(B.6) 

(B.7) 



In this approximation, the amplitudes are 



9\ 



^ _ . ^ . 1 1 

" " ■ 4/4 7 (2^)4 (l + q)2 - ml +ie P -ml+ie"^ 



where we used q^ = 0. 



x(<5,„+ieJ™V2")(F)(r + g'") 

dH 1 



9\ 



4/4 J (2^)4 {I + g)2 - ml +ie P-ml+: 



(B.8) 



(B.9) 
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Alternatively, one can derive the NLO triangle diagram contributions by using the HB formalism from 
the outset. Using the Feynman rules for the HB Lagrangians and propagators yields 



HB 



9\ 



16/^ J {2n 
X r|(a2 • I) 



{2l")e 



\ ^abc„ 



(l + g)^ — H- ze — + ie 



Vy 



HB 



16/4 J (2^ 



-T^{-<j2-{l + q)), 



(/ + g)2 — TO^ + P — + ie 



X (ai . [1 + ^)^^r1{~o, ■ I) , 
+ le 



(B.IO) 



(B.ll) 



HB 



Performing the isospin sums and using the identity {a ■ a){a ■ b) — a ■ h + ia ■ {a y. b) ^ we obtain 

" r dH 1 1 

4/4 J (27r)4 + g)2 - m2 + ie P^^Tt^T^ 

r- + +iCT2 ■ {l + q))] , 



(B.12) 



HB 



ITI-T2 



9l 



dH 1 

4/4 J (2^ + (j)2 - m2 + ie W 



1 



(B.13) 



The two sets of Eqs. (B.8), (B.9) and (B.12|, (B.13I obviously agree — as they should. 

Since {l + q)xl = qxl and since the integral over a term proportional to I will yield a result oc q, the 
spin dependent terms vanish and the sum of both triangles is given by 



V = V^ + Vb = -ixi-X2 



9\ 



dH 1 

2/4 J (27r)4 {I + qY - ml + ie 22-^^ + jg 

At this point, it might be of interest to compare with what other authors obtain for the NLO triangle 
diagrams. In the work of Ref. [S3], time-ordered perturbation theory (TOPT) (also known as old-fashioned 
PT) is used in which loop integrals extend only over the three space-dimensions. 

Thus, we will perform the integration over the time component in Eq. ( |B.14[ ) to make a comparison with 
TOPT possible. For convenience, we first apply the transformation I — ~ q), which yields 



V 



ITI-T2 



9l 



'ill 



1 



dH 

(27r)4 {V + g)2 - 4m2 + ie {V - g)2 - 4m2 + ie 



The poles in the lower half plane are at l"^ = u]± — ie with 



{V ± + 4m2 



and, so, we can write the amplitude as 



V = — ITt ■ T-2 



9\ 



dH' 



1 



8/4 ] (27r)4 (?'0)2 - cj2 + (;'0)2 _ ^^2 _^ 



(B.15) 

(B.16) 
(B.17) 



Now we perform the integration over I using the residue theorem 



V 



Tl ■ T2 



Tl • 1-2 



2 

9a 



9\ 



_3L 
16/^ 



dH' 
(27r)3 

dH' 
(2^ 

dH' 
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1 ' 


1 


\2c^^ 


n o 1 




1 CJ_ — CJ 


{^\ 


— uPl 2u;+cj_ 







(B.18) 



Except for a typo in Ref. |Slj and differences in notation, this agrees with the second term of Eq. (20) of 
Ref. [54] as well as the first term of Eq. (4.30) of Ref. [64], which both represent the triangle diagrams. 



Appendix B.1.2. Dimensional regularization of the triangles 

The integral is divergent and, therefore, requires regularization. In Ref. 1541 , this is done by introducing a 
Gaussian cutoff function exp(— P/A^) which, however, makes the results dependent on the cutoff parameter 
A. To avoid such cutoff dependence, we will use dimensional regularization which was first introduced by 't 
Hooft and Vcltman in 1972 [277] 



For this it is convenient to go back to Eq. (B.14) and apply a trick discovered by Feynman [230], 

dx 



1 

ab 



(B.19) 



;o [b+{a-b)x]^' 

which allows to write the product of propagators in terms of a linear combination. Applying Feynman's 



trick in Eq. (B.14 1 yields 



V 



dx 



dH 



V + l-q 



i(27r)4 [P + {2lq + q^)x - ml + ie] 



(B.20) 



The momentum integral is not yet quite ready for dimensional regularization because the integrand mixes 
three- and four-dimensional versions of the loop- momentum I. The dimension of the integration variable I 
must be consistent throughout the integrand. 

Formal (four-dimensional) covariance can be re-covered by re-introducing the four-vector — (1, 0, 0, 0) 
and re-writing the numerator as P -I- / • q = (/ • w)^ — P — I ■ q. After a shift of the integration variable 
1^1 + xq, well-known formulae for the dimensional regularization of covariant integrals can be applied [see, 
e.g., Eqs. (10.23) and (10.32)-(10.34) of Ref . 



Alternatively, one may go for consistency in three dimensions and reduce the momentum integral of 
Eq. ( B.20| ) to the three space-dimensions — which is what we will do here. This has the advantage that the 
integration variable observes an Euclidean metric from the outset (and there is no need for a Wick rotation). 
Moreover, the planar box diagram (see below) requires treatment in three dimensions. 

The integrand of Eq. (B.20 1 has two poles as evidenced through 



P + {2lq + q^)x -ml+ie = 



l° + ^P+ml + {2l-q + f)s 



where we used 5° = 0. The pole in the lower half-plane is 



f - \IP+ml + {2l-q + q^)x + le 

(B.21) 



and, using the residue theorem, one obtains 

V - 



P + ml + {21 ■ q + (p)x — ie , 

1^ + I- q 



(B.22) 



1 r df dH 

^''^^2t^], J J{2^^ {P-llY{P + l\Y 



dH IP + I 



(27r)3 4 (;0)3 



(B.23) 



To proceed towards dimensional regularization, we now generalize the three-dimensional integral to {D — 1) 
dimensions. Moreover, we introduce the renormalization scale A and multiply the integral by A**^^ such 
that the dimension of the amplitude stays the same for any D |228l 1231] , 



V 



dx 



d^-H P + l-q 



(2vr) 



[llf 



(B.24) 



^Accessible introductions into dimensional regularization can be found in Refs. I228II2291 and a very pedagogical presentation 
is provided in appendix A. 2 of the Scherer review 113 . 
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For convenience, we now divide all dimension-full quantities involved in the integral by A and denote all 
dimensionless quantities (except for the dimensionless momentum k) by a tilde; for example, q = q/X- 
Furthermore, we shift the integration variable by introducing the new variable k with Xk = I + xq, so that 
there are no terms linear in the integration variable in the denominator. Thus, we have now 

= \\J P+ml + q^x{l -x)-ie, (B.25) 

and our integral reads 

gl 1 f d^^^k {k-xqf + {k-xq)-q 

8/4 7o J (27r)^-i [k^ + Tnl+q^x{l- x)]^/^ 

Terms with odd powers of k vanish, 

The volume element in (Z? — 1) dimensions is 

d^-^k = k^-^dkdno-i (B.28) 



with 



^ d^D-i = . (B.29) 

Thus, 



27r(^-i)/2 



q\ 2A2 /■! , /-"^ „ k^ ^x{\~x)q'k^-^ 



Using [cf. for example, Eq. (A. 35) of Ref. 113 ] 



(A:2 + s)"/2 2r(f) 



rffc = ^ ' ^ , (B.31) 



we obtain 



^4 r(f )r(^)(47r)(^-i)/2 7o [ [ml + ipx{l - x)](2-^)/2 

We choose D = A ~ r] and then take the limit -q — > 0. To prepare for the latter, we expand the various 
expressions involved up to first order in 77: 

r(^) _ D-l_i-r, 



r(^) 2 2 



n^l = r(|) = --, + 0(,) 
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(47r) 



(4t 



7r)2 (47r)2 V 2 / 



1 



1- |ln(m2 + ^2,.(i_^)) 
1 - I ln(m2 + _ a,))] (^2 ^ ~2^(^ _ 



(B.33) 



Concerning the properties of the F-function, see, e.g., Ref. [232] : 7 = 0.5772 ... is Euler's constant; and we 
used a*^ — exp(elna) « 1 + elna + . . .. 
Hence 



Y 



Ti-T^^^ (l+'l ln(4^)) dx 

' 8/4 r(|) (4^)i V 2 ^ 



3 — 1] 



+ 7-1 + 0(7;)) l-'^lniml + fxil-x)) {ml + q^x{l-x)) 



-~x{l — x)q^ 
and neglecting terms of order 77 



7 + 0(7;)) l-'^lniml+fxil-x)) 



(B.34) 



^^•^^8/4r(|) (47r)i 



+7^^ ^Iq'^ + j dx[ ^ml + ^g''a::(l - x) ] \n{Thl + g^a;(l - a;)) 



h 7 - 1 - ln(47r) 



+ 12^^ 



(B.35) 



The integral over the logarithmic term is 



dx + -q x{l - x) \nim^ + q x{l - x)) 



4. 2 25 ,2 



- W^Q + 37?7^ + -q \n{m^) 
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+ -i8ml + 5f)L{q) 



6 



where 



L{q) — In „ ^ , Tt; 



g 27n^ 



Ami + 9^ ' E^iid g = jifl . 



(B.36) 
(B.37) 



Introducing the abbreviation 



R 



+ 7 - 1 - ln(47r) , 



the final result for the triangle diagrams at NLO in dimensional regularization reads 



9l 



384^2/4 



13 



(18771^ + 5q^)R - 4777^ - — g2 + 2(18777^ + 5q^) In ( ) + (167772 + lOq^)L{q) 



(B.38) 



(B.39) 

Renormalization in a Modified Minimal Subtraction scheme (AfS'-scheme) [228] amounts to omitting the 
i?-term. 



Appendix B.2. Football diagram 

Appendix B.2.1. The football in covariant perturbation theory 

The football diagram is shown in Fig. |B.2| In the HB formalism, one obtains 

dH 



V 



1 1 

i - 



2 16/4 J (27r)4 



(2Z" + g°)e°^'=r{= 



{I + q)"^ — 7772 +76 /2 — 7772 + It 
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-{2f + q^)€^'"^T^ , (B.40) 



iEp,p') {Ep,-p') 



{Ep ,p) {Ep, -p) 

Figure B.2: Two-pion exchange football diagram at NLO. 



where the factor ^ is a combinatoric factor and q ~ {0,p ' — p). Using ^"-bc^c^bad^d _ t2 and — 0, 



Y 



ITI-T2 



dH 



4/4 J (27r)4 (/ + qf ~ml + ie P - ml + le 



(B.41) 



It is again of interest to compare with what other authors obtain for the diagram under consideration. As 
discussed in conjunction with the triangle diagrams, in the work of Ref. |54| . time-ordered perturbation 
theory (TOPT) (also known as old-fashioned PT) is used in which loop integrals extend only over the three 
space-dimensions. Therefore, we will perform the integration over the time component in Eq. (B.41 1 to 
make a comparison with TOPT possible. Substituting I — ^{l' — q), 



V 



ITi ■T2 



ITI-T2 



'All 



dH 



/0\2 



16/4 J (27r)4 {I' + qY - Ami + i( - - + it 



dH' 



with 



16/4 / (27r)4 (^'0)2 - a;2 + (Z" 
= J{V±q)^ + Aml, 



Lot + it 



(B.42) 



(B.43) 



and applying the residue theorem yields 

V ^ Tl - To 



dH' 



1 



1 



16/4 / (27r)3 [ 2cj+ a;2 - a;2 cj2 _ 



Tl ■ T2 



1 



dH' 



uji — a;_ 



32/4/ (2^)3 0,2-0.2 

dH' 1 



32/4/ (2^)3 oj++oj^ 



(B.44) 



This amplitude for the football diagram does not agree with Ref. [53], where the result for the football is 
given by the first term in Eq. (20). In our notation, the result of Ref. ^54] and also of Ref. [64] [cf. Eq. (4.30) 
therein! is 



1 



ORK — 



Tl • T2 



dH' 



1 (w+-w_)2 



128/4 J (27r)3 Lo+uj- a;+ + uj. 



(B.45) 
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Appendix B.2.2. The proper calculation of the football in time-ordered perturbation theory 



We will explain now the reason for the discrepancy between Eqs. (B.44) and (B.45I. In covariant PT 
theory one assumes for the interaction Hamiltonian density Hi = — £/. This is strictly speaking not correct 
for derivative coupling, because in that case one has Hi = —Ci + additional noncovariant terms. For 
relatively simple cases of derivative coupling, the additional terms are of the contact type. Now, in covariant 
PT, there are also noncovariant contributions to the propagator, and it has been shown that for certain 
interactions with derivative coupling these two groups of additional terms cancel. So, as a procedural matter, 
in covariant PT, one can use Hi = —Ci and the usual covariant Feynman propagators, and everything comes 
out right (cf. e. g., Ref. [225], pp. 318-323 therein). 

In TOPT the story is different. Propagators are just non-covariant energy denominators, which cancel 
nothing. Therefore, the full interaction Hamiltonian including all additional terms has to be use to obtain 
the same result as in covariant PT. However, in Ref. |54j Hi = —Ci is assumed in conjunction with TOPT. 
Thus, the contributions from the additional contact terms in the Hamiltonian are missing in the calculation 
of Ref. '54] and that is the reason for the discrepancy with the covariant calculation, as we will show now. 

Using the notation of Ref. |225| . the time-component of the current associated with the Lagrangian 



density Eq. (A.43) reads 



Jfc" = -^iVr^iVe'^^^, (B.46) 
and the additional term in the interaction Hamiltonian is [cf. Eq. (12) of Ref. 51J] 

Hfit) = J d'x^{J^){J^) = J d'x^^NT-Ne-'''7:aNr-'Ne'''''''7ra'. (B.47) 
This generates the following leading contribution to the S'-matrix 

{p[p',\S\p^p2) = H) J dtHf''it) = -tJ d'x^^e^'^e'^''^\p[p',\NT^N7raNT^'Nna^pw^^ 

d'^x 2 Je''''"e"'^"'TfT2^'e*(p'i+PS^-Pi-P=^)=^(0|7ra7ra'|0) > (B.48) 



(27r)6 J 16/4 2 

where the factor of 2 comes from the permutation of the nucleon fields and the exchange term is left out. 
For the pion fields, we obtain 

ka^a' |0) = TTT^ / d^kd^k' , ^ (0| (a„(fc)e-*'=" + al(fc)e'*^) {aa>{k')e-''''^ + al^k'^^'A |0) 

^ ^ d'kd'k' ^ S^jk - k')5aa' ^ I ^^aa' ■ (B.49) 



Hence, 



(27r)3 J V4wfea;fc 
(p[p'2\S\piP2) = j^Jd^x^e'^''e'^''^'T^T^'e^^P'^+'^'^-P^-P^> J 



d^k 1 



(27r)3 2wfc 



5n.n' 



11 [ d^k 1 



2\ P +4mi 



(B.50) 
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where k — 1/2 and e^'^^e 



abc^abc ^c^c 
~1~2 



the T-matrix [of. Eqs. \kM\ and ( |A.35| ], which is 

1 



V 



add 



2 Ti • T2 . From the last equation, one can read off the contribution to 
dH 



Tl ■ T2 



1 



32/4 y (2^) 



(B.51) 



2\ P+4ml 



Obviously, this additional term doesn't depend on external momenta, only on m,, 
For the amplitude, calculated in TOPT and including all terms, we now get 

radd 



V = VoRK 

= VoRK 

= VoRK 

= VoRK 

= VoRK 

= TI-T2 



Tl ■ T2 

Tl ■ T2 

Tl ■ T2 

Tl • T2 
1 



1 



32/4 

1 

32/1 
1 

32/4 
1 

32/4 
dH' 



1 



2\ P+Aml 



dH 

dH' 1 

(27r)3 2tZ^ 

dH' (J_ 

(27r)3 
dH' {uj++uj^f 



4a;_ 



1 



(27r)-^ 0;+ + u)- Au^u). 
1 



(B.52) 
(B.53) 

(B.54) 
(B.55) 
(B.56) 
(B.57) 



32/4 J (2^Y + UJ- 

where we used the transformations I = I' + q and I' n- —I' . This result is now identical to the one obtained 
using covariant PT, Eq. (B.44|, as it should. We note that the term, Eq. (B.51), left out in the calculation 
of Ref. [54], is a contact term oc mj, which affects the m7r-dependence of the 2PE potential. Such terms are 
relevant for the charge-dependence of the nuclear force |188j and for considerations of the chiral limit [129] . 
If these more subtle aspects are not of interest, then omitting the contact Eq. (B.51 1 is acceptable, since 
contact terms with adjustable parameters are added to the theory anyhow (cf. Section 2.3 1. 

Appendix B.2.3. Dimensional regularization of the football 



We start from Eq. (B.41 1, 



V = iTl ■T2 



1 



dH 



1 



iff 



and apply the Feynman trick Eq. (|B.19|) 
V 



1 



Tl ■ T2 



4/4 J (27r)4 (/ + g)2 -ml + ie P 

iff 



dx 



dH 



i(27r)4 [P + {2lq + q^)x - ml + ie 



(B.58) 



(B.59) 



Using the residue theorem for the f integration, we obtain 



V = -X1-T2 



1 

471 



dx 



df dH 



iff 



Tl ■ T2 



4/4 



dx 



i(2^)4 (/0-?0)2(;0 + /0)2 
d^l 1 1 

(27r)3 il? 



(B.60) 



with 11 = ^JP + ml + (21- q + cp)x - ie. 



We treat this divergent integral by dimensional regularization (cf. Appendix B.1.2). For that purpose, 
we extend the integral to [D — 1) dimensions and shift the integration variable to \k — I + xq with A the 
renormalization scale: 



V 



Tl ■ T2 



1 



16/; 



dx 



jD-li 



1 



(27r)^-i 
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+ ml + q'^x{l - a;)]i/2 



(B.61) 



where quantities with a tilde are divided by A and, thus, dimensionless; and k is dimensionless by definition. 



The angular integration yields [cf. Eqs. (B.28I and (B.29)] 



V 



1 



2A^ 



dx I dk 



and applying Eq. (B.31) we obtain 



V = TI-T2 



1 



2 ' 



dx- 



1 



16/4 r(i)(47r)(^-l)/2 Jo ' [m2 +g22.(l_2.)](2-D)/2 



(B.62) 



(B.63) 



Choosing D = A — rj and using relations displayed in Eq. (B.33) results in 

>2 

V = TI-T2 

^ 2 

X [\x^^ '''"^"-^2 , „~2 

Jo 

and, neglecting terms of order 77, leads to 



^— — =- (1 + - InUn)] f-- + 7 - 1^ X 

16/4 r(i)(47r)t V 2 ^ V V ^ J 



1 - j Hml + fx{l ~ x))\ {ml + fx{l ~ x)) 



(B.64) 



V 



1 



16/4 r(i)(47r)t 



+ 7 - 1 - ln(47r) ] [ ml + -f 



6 



+ / dx {ml + q^x{l — x)')ln{ml + q'^x{l — x)) 
Jo 

The integral over the logarithmic term is 

5 



(B.65) 



dx {ml + fx{l - x)) \n{ml + fx{l - x)) = --ml - —f + 2ml + of H^n) + 



18 



q'^ + irhl 



with L{q) given in Eq. ( |B.37 ). 

Thus, for the football at NLO, we finally get 



y = TI-T2 



1 

3847r2/4 



Ir i6ml + q^) - 4to2 _ 5 2 ^ ^g^2 ^ ^2) 
2 



(B.66) 



(B.67) 



with i? as defined in Eq. (B.38) and w given in Eq. (B.37). 

Appendix B.3. Box and crossed box diagrams 

The planar box (pb) and crossed box (cb) diagrams together with our notation are shown in Fig. B.3 
and the corresponding amplitudes are given by 

d'^l i i 



9\ 



16/4 J (27r)4 qf - ml + ie g| - ml + ie 
xui[p ) (-92)^1 71^75 7:2 A^2 , ■ gi^i7M75-»i(pJ X 

y.U2{-p')qlT^lrl5 tr^^^'- (-gl)^2°7p75M2(-p) , 

— ivi J.J -\- te 



(B.68) 



V 



cb 



rf4/ 



i 



16/4 / (27r)4 qf-ml+ ie ql - m^ + ie 
xui(p ) (-92)^1 7^75 7:2 Af2 I ■ g"^r7M75-m(p) X 

XU2(-P ')(-gr)r2°7.75 .y^^- ,'^,2^^^r^- g2"2Sp75^^2(-p) : 
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(B.69) 



{Ep,p') 



{Ep, -P') 



{Ep, p) 



92 

b 



91 
*- 

a 



^N2 



{Ep, -p) 



H If 



s 92 
b N / 



X 



a / 
' 91 



kM, 



{Ep,p) 



{Ep, -p) 



kN, = {Ep - f , I) 
kN, = {Ep + f , -1) 



kN, = {Ep - f , V) 
kN,={Ep-f ,l-p-p') 



(a) (b) 

Figure B.3: Two-pion exchange planar box (a) and crossed box (b) diagrams at NLO. 



where we chose the loop momentum such that its time-component is the energy transferred by the pions 
and its space-components are the ones of an intermediate nucleon; fcjvi — {Ep — I), fcjva = (Ep + —I), 
qi = {f,p — I), and q2 = {f,p' — I), with Ep — ^J-p ^ -|- M^, for the planar box diagram. For the crossed 
box, only the momentum of nucleon 2 is different, k-^^ — {Ep — P ,1 — p — p'). The poles are located at 



11 = EpTEi±ie 

1° - ~EpTEi±ie 

^ = EpTEi±ie 

Zg — =pwi ± ie 



(B.70) 



with Loi = ^{p-iy + ml, L02 = ^{p' -If +ml, El = ^t^ + M^, and Ej = ^ {I - p - p + M%. The 

first three lines state the nucleon poles and the last two lines describe the pion poles. I2 and I2 are the poles 
for nucleon 2 in the box and crossed box diagrams, respectively. We apply the Feynman trick Eq. (B.19I to 
the pion propagators and obtain 



dH 



{p + Ep)^-{q + Epr'^'^'''"'^ P'' 



(B.71) 



V 



'cb 



9\ 



^Ul{p')q2ly 



(3 + 2xi-T2)^ dx J 



dH 



J (2^)4 [/02 _ ^2 + (^2 _ ^2)^ + ^2 

(;o„^^)2_(;o„^^)2'?i7.-i(p)x 

^"^^"^ )'^i^^ao-i?,)2-(ro_i,^)2'^27p"2(-p), 
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(B.72) 



where we used rfr{'T2rf = 3 — 2ti • T2 and Tfrfrf = 3 + 2xi • T2. The pion poles are now given by 



ll^ = ^^{l-x)u:l + xujl±ie. (B.73) 
Performing the 1/Mm expansion yields 

ui^Mn,-Mn)-i^ui « 5,,;^^,2Mjv((5,j+ie'-'V[), (B.74) 
M27r(^iv, - Mjv)7pU2 « 5„„(5p„2Mjv(<5™„ + ie™"V^), (B.75) 

(B.76) 
(B.77) 



1 1 11 

{P ~ - {10 - Epf {10 - E^Y - ill - E^Y AMj, {P ~ te^ ' 



where = ^prJ^ — h ie, which avoids the pinch singularity (cf. discussion on p. 7 of Ref. f51'). Hence, 



In both diagrams we have a contribution from the pion poles, but in the box diagram there is also a 
contribution from the nucleon poles which is enhanced in the l/M^r expansion and gives the iterated IPE. 
We work this out by closing the P contour integral in the lower half-plane and, collecting the contribution 
from the pole at P = results in 



Using 



1 1 



+ *e'^ V[)g^"g5'(<5™„ + ze^'^a^^) = 5^-{v' - 1)5^ ■ (p - l)a2 • (p" - 1)02 ■ [p- , (B.82) 

1 Mn 



2< -P + ie 

and changing / — >■ — /, we obtain 



,pb _ 9aMn o_ _ ^ /■ t?i-a+p')'?2-(;+p')f?i-G+P)^2-a+p) 



(B.83) 



~ 16/^ ~ '^^'^^> j (2^ (p2 _ 12 + j,)((r+ p)2 + ral){{l+p'Y + m2) 

which is the LO iterated IPE. 

Notice that the contribution from the nucleon pole in the crossed box diagram is zero. 

Next, we calculate the contribution from the pion poles closing the P contour integral in the lower 
half-plane with the pole of interest located at P — l'^^ — \/ {\ — x)uj\ + — ie: 

yph ^ ^9a T2) C dx ^ (ffi • 92 + ■ ((72 X • 92 +i(T2 • (92 X gi)) 



64/r 'Jo J (27r)3 {I0,j 



\5 

(B.85) 

^ = -d^(3 + 2Ti-r2) / dx 



rch 3ff^ /q , ^ \ [ .1 f (gi ■ 92 + i^i ■ {q-i x gi))(gi ■ 92 - ^^2 ■ (92 x gl)) 



34 y 

(B.86) 
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We treat the divergent integrals by dimensional regularization (cf. Appendix B.1.21. For that purpose, we 
extend the integrals to {D — 1) dimensions and shift the integration variable to Xk = p — I + xq with A the 
renormalization scale and q — p'—p: 



dx 



k {qi ■ q2 + iSi ■ (92 X qi)){qi ■ q2 + 1^2 ■ (92 X qi)) 



,D-1 



[fc2 4- rfi^ + (p-x{\ — x)]^/"^ 



(B.87) 



V 



cb 



dx 



d k {qi- q2+ «<?i • (92 X qi)){qi ■ q2 - i(J2 ■ (92 X qi)) 



i2n) 



D-l 



[k^ + rn^ + q'^x{l — x)]^^^ 



(B.88) 



where quantities with a tilde are divided by A and, thus, dimensionless; and k is dimensionless by definition. 
Since 



hence 



qi = k- xq, 

q2 = k-{x-l)q, 

91-92 = k"^ + q'^x{x — 1) - {2x ~ l)k ■ q , 
92 X (71 = q X k . 



Terms in odd powers of k vanish. Furthermore, because of 



we have 



Whence 



d^-'k 


fc j kj 




(27r)^-i [P- 


-f r7i2 + q'^x{l — x)]^/^ 


D - 1 




{k-qf 


1 


(27r)^-i [fc2 - 


r to2 + q22;(l _ x)]^/2 


D - 1 




(fc • 9)iCTi ■ {q'x k) _ 


0, 


(27r)^-i [fc2- 


r r7i2 + ^2^;^ _ a;)]'''/2 


d^-^k iai 


• (q X k)i(72 ■ (q X k) 


1 




h ?7i2 -f q'^x(l — a;)]'''/2 


D - 1 









rf^-lfc fc2 

£1-17 (27r)^-i [fc2+m2 _^ ^23.(1 _^)]5/2 



d^-^k 



d^-^k i^{ai X q) ■ (^2 X g)/c2 
D-1 J (27r)^-i [PT^?4T^Ml^^2?)p72 



r>-i 



(27r) 

[fc2 + 5~2^(^ _ i)]2 ^ [(2^ _ 1)2^^2 _ ^ ~) . (^^ ^ 



[A:2 + ml + q^x{l - x)]^/^ 
(3 + 2x,.X2)A^y^ d^Jj^^ 



64/4 

[fc2 + g2^(2; ^ 1)]2 + [(2a; - 1)2^2 + (^^ ^ g) • (^2 X g)] 
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(B.91) 
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(B.93) 



(B.94) 



(B.95) 



(B.96) 



[k^ + rfi2 + g2x(l — a::)]^/2 
Using ((7i X q) ■ ((72 X g) = A~2[g*2^^ . _ . q){(j2 ■ q)], the sum of the two diagrams is given by 

^pb^-^cb ^ Wc{Ti-T2) + VsicTi-a2) + VT{cTi-q)ia2-q) (B.97) 
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with 



Wc ^ - 



16/.* 



A y2 



dx 



jD-1 



Vt -- 

32/4 

Vs - -q'^Vr 



dx 



(27r)^-i 
d^-^k 



[k^ + rn^ + q'^x{l — x)]^/^ 



Applying the volume element in {D — 1) dimensions, Eqs. (B.28) and (B.29), we find 



[A:2 + + g22;(^i _ x)]^/^ 
2 1 f°° 





32/4 r(i2_i)(4^)(i3-i)/2 



and performing the fc-integration with the help of Eq. (B.31 1, we obtain 
3gl 



Wr 



A2 



dx 



r(^)r(^) 



+ 



16/4 r(|)r(^)(4^)(^-i)/2 Jo 
r(i2±i)r(4^)g2 r2:,(^ _ 1) + 



[m2 +q22.(l_2.)](2-D)/2 

r(^)r(g=^)gV(x-i)- 

^ [m2+ 523.(1 _ 2.)] (6-D)/2 



[m2 +q;22.(i_2;)](4--D)/2 

9gi 1 f\ n^n^) 

32/4 r(|)r(i^^)(47r)(^-i)/2 - 1 io ^ ['^iJ + - 



We choose D = 4 - 77 and take the limit 77 ^ [cf. Eq. (B.33)], 



Wr 



4 rl 



+ ASrhi + 6g^(-l + 22a; - 22a;-') 



12g4a;2(x- 1)2 
hi + q'^x{l — x) 



[4:5ml + f{-6 + 105x - 105x^)] \n{rhl + fx{l - x))} , 



Vt = - 



3gl 

128^2/4 



dx {i? + 1 + ln(m^ + 9^a;(l - a;))} 



(B.98) 

(B.99) 
(B.lOO) 



(B.lOl) 
(B.102) 



(B.103) 
(B.104) 



(B.105) 
(B.106) 

with R as defined in Eq. (B.38). And performing the a;-integration, we finally get for the planar box plus 



crossed box diagrams at NLO 
Wc = - 

4 

Vt = - 



9\ 



3847r2/4 



1 



ri? + ln 



(x) 



{90ml + 23^2) + IQml + -q^ 





48m4 



20mi + 23q^ + ^ L{q) 



V _ 3gl 
^^3- 



647r2/4 



with w and L{q) given in Eq. (B.37) 



w 
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(B.107) 
(B.108) 



Appendix B.4-. Summary of 2PE contributions at NLO 

Adding up all 2PE contributions at NLO, namely, the box and crossed boxes Eqs. (B.107I and (B.108), 
the triangles Eq. (B.39), and the football Eq. (B.67), we obtain 



Wc 



384^2^4 1 1 4"^' (%i - 4gi - 1) + q\2?,g\ - lOgi - 1) + 
[6m2 (1551 " 6.gl - 1) + q^{2ig\ I0g\ 1)] In (^) 



-Aml{Ag\ +g\ + l) + ^-{bg\ + 2(ig\ + 5) 



M"''7r 



Vt 



35 



(B.109) 
(B.llO) 



with w and L{q) given in Eq. (B.37), and where we applied renormalization in a Modified Minimal Sub- 
traction scheme (M5'-scheme) ^228) . i. e., the (infinite) i?-terms are subtracted. The results fully agree 
with Ref. [55]. If one is not interested in subtle aspects, like charge-dependence |188j or the chiral limit 
(tojt 0) |129j . one may omit the polynomial terms in the above expressions [cf. Eqs. (4.9) and (4.10l], 
since contact terms are added anyhow (cf. Section 4.3) which have the same mathematical structure and, 



therefore, can absorb the polynomial terms. Thus, the contacts are getting renormalized by the polynomials 
that result from dimensional regularization. 



Appendix C. Two-pion exchange contributions to the 2NF at NNLO 



Two-pion exchange diagrams that contribute at NNLO were shown in Fig. [4] These diagrams differ from 
the NL O diagrams. Fig. [3| by one insertion from the dimension-two Lagrangian, the vertices of which are 
given in [Appendix A.O There are two kinds of second order vertices: The relativistic corrections of the 
leading order Lagrangian, which are proportional to l/Mjy, Eqs. ( A.48 )-( A.50 ), and new contact interactions 
proportional to the LECs c^, Eqs. (A. 51 ) and (A. 52). Since the shapes of the NNLO diagrams are the same 



as in the NLO case, the evaluation is very similar to what we presented in detail in [Appendix B| for the NLO 
diagrams. The results for the irreducible 2PE at third order (NNLO) as derived by the Munich group [55] 
were given in Sec. 4.1.2 Eqs. (4.13)-(4.18|. Note that the irreducible 2PE depends on what is chosen for 
the (reducible) iterated IFF. For this, the Munich group applies the expression Eq. (4.25). Alternatively, 



one may also use the form Eq. (4.26) implied by the BbS scheme |133j . which differs in third and higher 
orders from the Munich expression. Therefore, when the BbS formalism is used, a (irreducible) correction 
term has to be added to the Munich irreducible 2PE in higher orders. In terms of the notation introduced 
in Eqs. (4.25) and (4.26), this correction term is given by: 



Tr(KBW) _ T^(EM) 



d^p" 

(2^ \Ep Er, 



M' 



p2 _ pll'^ _|_ 



d^p" _ 
(2^)3 Ep„Ep{Ep„ +Ep) 
d^p" 1 



^ vUp\p")ViAp",p} (c.i) 

VAp',p")V,^{p",P) (C.2) 



(27r)3 2Mn 



Vi^{p',p")Vi^{p",p). 



(C.3) 



with the la st li ne showing the correction in third order or NNLO, which is what we want to calculate here. 
Using Eq. (4.5) and the relation Ti • T2 Ti • T2 = 3 — 2 Ti • T2, and defining / = p — p" , yields 



, .(KBW) _ .(EM) 
''^27r.it *^2ir,it 



9a fo o-^ \ f '^^^ ai- (l + q) a2- {l + qj ^i-l ^2-1 ,^ 
— (3-2ri -Ti) / 7;;^ ^ ,a , _2 V<-4) 



32/4Af 



N 



(27r)3 (l + g)2 + ^2 p + ml 



(3 - 2 ■ X2) ['d.[^, + t'^: 1 - 1? - ' (C.5) 



32/4M 



N 



(27r)3 [P + {2r-q + q^)x + ml]^ 
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where, in the last hne, we introduced the Feynma n trick. We now ex tend the integral to {D — 1) dimensions 
and make the change of variables Xk = I + xq (cf. Appendix B.1.2), 



. .(KBW) _ .(EM) 



9\ 



32/4A/ 



N 



-(3- 2ri • T2)X^ / dx 



{2n 



,D-1 



Since 



(7i • [fc + (1 - x)q] (72 • [fc + (1 — x)q] Si - {k — xq) (T2 • (fc — xq) 
[fc2 + fn^ + q^x{\ — x)Y 



i^i • [fc + (1 — x)q\ Ui ■ [k — xq) = fc^ + q^x{x — 1) — {2x — l)k ■ q + i(Ti ■ {q x k) 



(C.6) 



(C.7) 



we obtain, making use of Eq. (B.94), 



..(KBW) _ , .(EM) _ 



^4 /"l f d^-^k 



[fc^ + glT(.T - 1)]^ + ^[(2x - - (ai X g)(a2 x g)] 
[fc2 + toJ + q^x{l — x)Y 



Applying (cti x q){S2 x q) = X "^[q^ai ■ S2 - [ffi ■ q){^2 ■ q)], leads to 



Vc 
Wc 

Vt 

Vs 
Wt 
Ws 



/ dx 



32/4M^ Jo 



^^'^ X I dx 



d^-^k [fc' + ex{x - + ^(2x - l)2g2 



(27r 



[fc2 + ml + q'^x{l - 



32/4M^ Jo 
-q^Vr 

-q^Wr ■ 



D-l 



(27r)^-i [k^ + ml + q^x{l - x)] 



(C.8) 

(C.9) 
(C.IO) 

(C.ll) 
(C.12) 
(C.13) 
(C.14) 



Using the volume element in {D — 1) dimensions, Eqs. (B.28I and (B.29), we find 



Vt 



2A3 

32/4A% r(^)(4^(^-i)/2 Jo Jo 
3g\ 2X 

2 



dx dk- 



^ ^.D~2 l2x{x - 1) + + fc^-2q4^2(,^ _ ^^2 



dx dk-, 



[fc2 + ml + q22;(i — x)]^ 
k^ 



32f^MNT{^)i4TT)(D-i)/2D-lJo Jo [k^ + ml + q^x{l - x)]^ 
Performing the fc-integration yields 



(C.15) 
(C.16) 
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r(^)r(i^) 
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[ml + q^x{l-x)](-'-D)/2 

' dx- r(^)r(^) 



(C.17) 
(C.18) 



We choose D = A — ij and take the Hmit 77 — 0, 
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^ 5mt - m2g2(i „ 20x + 20x'^) + q'^xi-l + 17x - 32x'^ + Wx^) 



Vt 



2567r/4M^ Jo 

A 



\J'm\ + (f-x{\ — x) 



2567r/^MAr Jo 
and, performing the x-integration, we finally obtain 



Vc ^ - 
Vt 



3/ 



^ (m^+w^A(g)). 
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Appendix D. Two-pion exchange contributions to the 2NF at N^LO 

The fourth order 2PE contributions consist of two classes: the one-loop (Fig. [5| and the two-loop 
diagrams (Fig.|6|. 

Appendix D.l. One-loop diagrams 

This large pool of diagrams can be analyzed in a systematic way by introducing the following well-defined 
subdivisions. 

Appendix D.1.1. contributions. 

The only contribution of this kind comes from the football diagram with both vertices proportional to 
Ci (first row of Fig. [5]). One obtains |60j : 



Vc 



Wt 



3L(g) 



''2 2 I ~2 A 2 

—w + C3W — Acinij. 
6 



"-2 4 

45 



with L{q) and w defined in Eq. (B.37) and 



967r2/4 



(D.l) 
(D.2) 

(D.3) 



[(C2 - 6c3)g4 + 4(6ci + C2 - 'ic^)q^ml + 6(c2 - 2c;)m^^ + 24(2ci + c:,)ml 



Appendix D.l. 2. Ci/M^ contributions. 

This class consists of diagrams with one vertex proportional to ci and one l/Mjf correction. A few 
graphs that are representative for this class are shown in the second row of Fig. [5] Symbols with a large 
solid dot and an open circle denote l/M^ corrections of vertices proportional to c^. They are part of £^^'^, 
Eq. (2.671. The result for this group of diagrams is [60] : 

Vc 

Wc 
Wt 
Vls 
Wls 



Ciq^L(.q) 
1927r2MA,/4 



[gl[Sml + bq^) + w^] , 



CiL{q) 



192TTUlNf* 



[g^(16m2 +7^2)_^2] ^ 



C2 gA 



8^2Ma,/4 

CiL{q) 



w^L{q), 
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Appendix D.1.3. 



\/Mjj corrections. 



These are relativistic corrections of the leading order 27r exchange diagrams. Typical examples 

for t his la rge class are s hown in row 3-6 of Fig. [5] This time, there is no correction from the iterated IPE, 



Eq. (4.25) or Eq. (4.26), since the expansion of the factor M^/Ep does not create a term proportional to 



The total result for this class is 



Vc 
Wc 



9\ 



+ 8to>-2 -q'~ 2ml) + 




(D.9) 



+4:g\ (^k^{20ml + 7q^ - Wmlw-^) + I6mlw-^ + \2m%w-'^ - 4m*q^w-^ - 5q 
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Wt = 

Vls = 




q^ 32^2^2^/4 



g2 ^ 15367r2Af2,/4 
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/ 17 
4g\ l7ml + —q^+4miw- 



32g\ (^ml + ^g^^ 



4^2M2^/4 V32 
2567r2M2,/4 
327r2M2,/4 ■ 
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165, 
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— 6rn, 
(D.IO) 
(D.ll) 

(D.12) 
(D.13) 

(D.14) 

(D.15) 



Appendix D.2. Two-loop contributions. 

The two-loop contributions are quite involved. In Fig. |6j we attempt a graphical representation of this 
class. The gray disk stands for all one-loop ttN graphs which are shown in some detail in the lower part 
of the figure. Not all of the numerous graphs are displayed. Some of the missing ones are obtained by 
permutation of the vertices along the nucleon line, others by inverting initi al and final states. Vertices 
denoted by a small dot are from the leading order Lagrangian Eq. ( 2.65 1. The solid square represents 



vertices proportional to the LECs di introduced in C^^l^ which is part of £^ Eq. (|2.67[). The di vertices 



-A=2 



occur actually in one-loop NN diagrams, but we list them among the two-loop NN contributions because 
they are needed to absorb divergences generated by one-loop ttN graphs. Using techniques from dispersion 
theory, Kaiser [SU] calculated the imaginary parts of the TV TV amplitudes, Im Va{i^j) and Im Wa{iy^), which 
result from analytic continuation to time-like momentum transfer q — z/i — 0^ with /i > 2m^. From this, 
the momentum-space amplitudes Va{q) and Wa{q) are obtained via the subtracted dispersion relations: 



VT[q) 



2q^ 



IinVc.5(»M) 
g2) 



TT 



lmVT{i^J,) 



/i3(^2 _|_ q2^ ' 



(D.16) 
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and similarly for Wq s t- 
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Figure D.l: The effect of individual fourtfi-order contributions on the neutron-proton phase shifts in some selected peripheral 
partial waves. The individual contributions are added up successively in the order given i n parentheses ne xt to each curve. 
Curve (1) is NNLO ("N2LO") and curve (6) completes N'^LO. For further explanations, sec [Appendix D.3| Empirical phase 
shifts (solid dots and open circles) as in Fig.ls] 



where 



... 1 q 
Ml) = arctan 



and 



(D.28) 



Note that the analytic solutions hold modulo polynomials. We have checked the importance of those 
contributions where we could not find an analytic solution and where, therefore, the integrations have to be 
performed numerically. It turns out that the combined effect on NN phase shifts from W^'^ 
is smaller than 0.1 deg in F and G waves and smaller than 0.01 deg in H waves, at Tiab = 
less at lower energies). This renders these contributions negligible. Therefore, we omit Wq \ V^^\ and Vg'''' 



y^"), and yj") 
300 MeV (and 



in the const ructio n of c hiral NN potentials at order N'^LO (Section 4.6.31 



In Eqs. (D.20) and (D.24), we use the scale-independent LECs, di, which are obtained by combining the 
scale-dependent ones, d[(A), with the chiral logarithm, ln(TO7r/A), or equivalently di — d^{mTr)- For more 
details about this issue, see Ref. |123| . 

Appendix D.3. Impact of individual N^LO contributions on peripheral phase shifts 

The fourth order is obviously very diverse. Here, we will show how the individual fourth-order con- 



tributions impact NN phase shifts in peripheral partial waves. For this purpose, we display in Fig. D.l 
phase shifts for four important peripheral partial waves, namely, ^F^, ^F^, ■^F4, and ^G^. In each frame, the 
following curves are shown: 



82 



(1) NNLO ("N2L0"^ 



(2) The previous curve plus the cf graph, first row of Fig. [sj Eqs. (D.l) and (D.2|, denoted by 'c2' in the 
figure. 



(3) The previous curve phis the Ci/AlN contributions (denoted by 'c/M'), second row of Fig. [5) Eqs. ( [dII 



(4) The previous curve plus the corrections ('1/M2'), row three to six of Fig. [s] Eqs. (D.9 )-(D.15 ). 



(5) The previous curve plus the two-loop cont ributions with out the terms proportional to di ('2-L'); i.e. 
Fig.je] but without the sohd square; Eqs. ( |D.18[ )-( |d^ , but with all di = 0. 



(6) The previous curve plus the terms proportional to di (denoted by 'd' in the figure) with the parameters 
given in Table [2j column 'Peripheral perturbative NN\ 

In summary, the various curves add up successively the individual N^LO contributions in the order indicated 
in the curve label. The last curve in this series, curve (6), is the full N'^LO result. 

The cf graph generates large attraction in all partial waves [cf. differences between curves (1) and (2) in 
Fig. |D.l] . This attraction is compensated by repulsion from the Ci/Mf^ diagrams, in most partial waves; the 
exception is ^F^ where Ci/Mpf adds more attraction [curve (3)]. The corrections [difference between 

curves (3) and (4)] are typically small. Finally, the two-loop contributions create substantial repulsion in 
^F^ and ^G^ which brings ^F^ into good agreement with the data while causing a discrepancy for ■^6*5. In 
^F-^ and '^-Kt, there are large cancelations between the 'pure' two-loop graphs and the di terms, making the 
net two-loop contribution rather small. 



A pivotal role in the above game is played by Ws, Eq. (D.6I, from the Ci/M^ group. This attractive term 



receives a factor of 9 in ^i^3, a factor (--3) in "^Gs, and a factor of 1 in '^^3 and '^^4. Thus, this contribution 
is very attractive in and repulsive in ^G^. The latter is the reason for the overcompensation of the c| 
graph by the CijMjq contribution in ^G^ which is why the final N'^LO result in this partial wave comes out 
too repulsive. One can expect that l/Af^v corrections that occur at order five or six will resolve this problem. 

Before finishing this appendix, we like to point out that the problem with the ^G^ is not as dramatic as 
it may appear from the phase shift plots — for two reasons. First, the "^Gs phase shifts are about one order of 
magnitude smaller than the F and most of the other G phases. Thus, in absolute terms, the discrepancies 
seen in ^G^ are small. In a certain sense, we are looking at 'higher order noise' under a magnifying glass. 
Second, the ^G^ partial wave contributes 0.06 MeV to the energy per nucleon in nuclear matter, the total 
of which is —16 MeV. Consequently, small discrepancies in the reproduction of ^Gs by a potential will 
have negligible influence on the microscopic nuclear structure predictions obtained with that potential. 



Appendix E. Fourth order NN contact potential: Partial wave decomposition 



The contact potential of order four, Eq. (|4.43|), decomposes into partial- waves as follows 



p 



= D.p^{p'^p + p'p^) 
= D.p^{p'^p + p'p^) 



Dso^p^ 

D3Si-^DiP +D3Si-3DiP P 
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V^'H^Ds) = D3^3p'V (E.l) 



with the coefBcients given by 
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Tabl e F.l : LECs for two N^LO fits by the Idafio group [68] using A = 500 MeV and 600 MeV in the regulator function f (p',p) , 
Eq. | |4.63| . The Cg and Ct of the zero-order counterterms given in Eq. \4.38\ are in units of 10* GcV~^; the d, Eq. | |4.40| , 
in lO^^V""'; and the D,, Eq. ( |4.43| , in 10"' GeV"''. 



LEG 


A = 500 MeV 


A = 600 McV 




-0.009991 


-0.009943 




-0.010011 


-0.009949 




-0.010028 


-0.009955 


\_yrp 


0.000523 


0.000695 




0.000543 


0.000701 




0.000561 


0.000707 


Ci 


0.051949 


0.046433 


C2 


0.163034 


0.174354 


C3 


0.003249 


0.006562 


C4 


-0.048954 


-0.050066 




-0.075081 


-0.086978 




-0.013343 


-0.013639 


C7 


-0.225500 


-0.214188 




-0.016674 


-0.004441 


D2 


2.480231 


2.751960 


n 


n Qi ^OAn 


U.UOuDOU 


Di 


-0.811591 


-0.061000 


D5 


0.138064 


0.126881 


De 


1.249801 


1.266717 


D7 


0.148074 


0.268796 


Ds 


-0.153405 


-0.239967 


D9 


-0.516607 


-0.538043 


Dio 


2.379565 


2.429102 


Dii 


-0.125289 


-0.132228 


D12 


0.081807 


-0.053835 


Di3 


0.008037 


-0.081804 


Du 


-1.393312 


-1.204164 


Di5 


0.164684 


-0.005571 



Appendix F. Parameters of N^LO NN potentials 



This appendix provides detailed information on parameters involved in the Idaho N^LO NN poten- 
tials IMl- The ttN LECs were shown already in Table [2| (colu mn 'NN pot ential') . In T able [rT| the LECs 



are listed, which are associated with the contact terms, Eqs. (4.381, (4.40), and (4.431, and, in Table F.2 



the contact LECs are given in terms of partial-wave parameters. Notice that through Eqs. (4.391, (4.42) 



and (E.2) there is a one-to-one correspondence between these two representations of the contact LECs. The 



partial wave parameters appear one order of magnitude larger, because they are multiplied by (47r) in the 
partial-wave decomposition. The partial-wave contact parameters of the Jiilich N'^LO potentials can be 
found in Ref. |171| . There are similarities between the parameter sets. 

At this point, one may raise the question of the naturalness of the LECs. LECs may be perceived as 
natural if they are roughly of the following sizes: 



Cs,T 
C, 



1 _ 

72 ~ 

J TT 

1 



0.01 



10"^ 
GeV2 ' 
10^ 



0.05 



GeV4 ' 



(F.l) 
(F.2) 
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Table F.2: Partial-wa ve LE Cs for two N^LO fits by the Idaho group [68] using A = 500 Me V and 600 MeV in the regulator 
function f(p' , p), E g. 1 4.63 1. The d of the zero-order part ial-w ave counterterms given in Eq. (4.39i are in units of 10'' GeV~^; 
the d, Eq. (4.4ll, in lO^* GeV—*; and the D,, D^, Eq. |E.iL in 10'' GeV-'^. The last column lists the exponent n of the 



regulator function, which is applied to the corresponding partial-wave counterterm. 



Partial-wave LEG 

D^Po 
CiPi 



Di 



Pi 

3Pi 



Ds 



Pi 



Si 



D3 



D3D2 



A = 500 MeV 
-0.145286 
-0.146285 
-0.147167 
2.380 
-2.545 
-16.00 
1.487 
0.245 
0.656 

5.25 
-0.630 

2.35 
-0.118972496 
0.760 

7.00 

6.55 
-2.80 
0.826 

2.25 

6.61 
-1.770 
-1.46 
-0.538 
2.295 
-0.465 

5.66 



A = 600 MeV 
-0.151165 
-0.151467 
-0.151745 
2.200 
-4.890 

-5.84 

1.548 
-0.215 
0.790 

4.40 
-0.488 

3.24 
-0.116210 
0.775 
4.8004 
10.8654 

-2.35 
0.796 

2.86 

5.58 
-1.764 

-1.27 
-0.548 
2.554 
-0.525 

6.26 



n 
3 
3 
3 
2 
2 
2 
2 
3 
2 
2 
2 
4 
3 
2 
2 
2 
2 
2 
2 
2 
4 
2 
2 
2 
4 

2, S*^ 



/(p',p) = 0.5{exp[-(p'/A)4 - (p/A)4] + exp[-(p7A)6 - {p/Af]} is applied. 



0.2 



10^ 
GeV6 



(F.3) 



where we used A = 0.5 GeV for the numerical estimates. Gomparison with Table F.l reveals that the values 
of the contact LEGs that result from a fit to the NN data are in general natural, indeed. Exceptions are 
D2 and DiQ which are an order of magnitude too large. 



CI 



nn 
Sal 



of Table 



F.2 



The charge-dependent LO contact parameters Cfg^ 
follows. Based upon the isospin-symmetric Lagrangian Eq. (|2.6l| and the isospin violatmg ones, Eqs. (|4.80l 



and C?s^ 



can be analyzed as 



and (4.811, the charge-dependent non-derivative contact parameters in the ^Sq state are given by 



a 



sym 

Cnp ^sym 



CciB + C'CSB : 
CciB — C'cSB : 
CciB , 



(F.4) 
(F.5) 
(F.6) 



where Cgym is the same as C15 of Eq. (4.39). The charge-symmetric and the GSB and GIB contributions 
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can then be obtained from 



C*sym 
CcSB 
CciB 



1 

2 

1 
2 



2 (^^PP ~^ Cfin 



" Cn 



2 ^^PP ^ Cnn^ Cn 



Using the numbers from the A — 500 MeV fit hsted in Table F.2 the following values are produced: 



CSB 



= -0.1465 



0.0005 



10* 
GeV2 ' 
10^ 



CciB = 0.0007 



GeV2 
lO'^ 
GeV2 



(F.7) 
(F.8) 
(F.9) 

(F.IO) 
(F.ll) 
(F.12) 



Note that these are partial-wave parameters, which carry a factor of (47r) as compared to the parameters 
of, e.g., Eq. (F.ll. As discussed in Section 4.6.2 the Idaho N'^LO potential j68] does not include explicitly 



the CIB effect from the pion-mass splitting in the NLO 2PE. Thus, the above value for Cqib is not the pure 
contact contribution since it includes a simulation of this omitted (small) effect (causing lS.aciB ~ —0.5 
fm). Moreover, also the CSB effect from nucleon-mass difference in the NLO 2PE is not calculated explicitly 
and, thus, effectively contained in Ccsb- The exact size of this effect when evaluated in ChPT at NLO is 
not known, but could be a large fraction of the total CSB |179j . 

The regulator function f{p'-,p), Eq. (4.63), has two parameters, namely, the 'cutoff' A and the power 
parameter n. Fits were made for A = 500 MeV and 600 MeV [denoted by N3LO(500) and N3LO(600)]. The 
parameter n has to be chosen such that the higher powers generated by the multiplication with the regulator 



cf. Eq. (4.64|] are beyond the order at which we are working, which is v = A. In short, the requirement is 



2ri + i/i > 4 , 



(F.13) 



where Vi denotes the order of a specific contribution to the potential. Thus, n > 3 for LO contributions, 
n > 2 for NLO contributions, etc.. In the case of the Idaho N^LO potentials, n = 4 is used for the IPE and 
n — 2 for all 2PE contributions (which are NLO or higher). For the contact contributions, n is chosen in a 



partial-wave dependent way and in accordance with Eq. (F.13), see Table F.2 
Phase shifts of NN scattering as produced by the Idaho N^O (500) potential 



are shown in Tables 



F.3 



F.6 up to laboratory energies of 300 MeV. 
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Table F.3: pp phase shifts (in degrees) by the Idaho N3LO(500) potential [68] , 



Tiab (MeV) 


'So 




'Pi 








£2 








1 


32.76 


0.13 


-0.08 


0.00 


0.01 


0.00 


0.00 


0.00 


0.00 


0.00 


5 


54.71 


1.58 


-0.90 


0.04 


0.21 


0.00 


-0.05 


0.00 


0.00 


0.00 


10 


55.01 


3.73 


-2.05 


0.17 


0.64 


0.01 


-0.20 


-0.03 


0.00 


0.00 


25 


48.37 


8.57 


-4.91 


0.70 


2.47 


0.10 


-0.80 


-0.23 


0.04 


0.02 


50 


38.75 


11.47 


-8.31 


1.70 


5.87 


0.33 


-1.65 


-0.69 


0.15 


0.10 


100 


25.74 


9.30 


-13.39 


3.68 


11.11 


0.80 


-2.54 


-1.48 


0.41 


0.42 


150 


16.66 


4.34 


-17.63 


5.28 


14.01 


1.22 


-2.83 


-1.96 


0.66 


0.87 


200 


8.73 


-0.87 


-21.48 


6.59 


15.67 


1.56 


-2.76 


-2.16 


0.88 


1.37 


250 


0.51 


-5.78 


-25.43 


7.95 


16.85 


1.81 


-2.35 


-2.11 


1.05 


1.83 


300 


-8.53 


-10.43 


-29.54 


9.16 


17.73 


1.94 


-1.67 


-1.83 


1.14 


2.17 



Table F.4: nn phase shifts (in degrees) by the Idaho N''LO(500) potential I68| . 



Tiab (MeV) 


'So 


'Pn 


'Pi 




'P2 




£2 








1 


57.51 


0.21 


-0.12 


0.00 


0.02 


0.00 


0.00 


0.00 


0.00 


0.00 


5 


60.86 


1.85 


-1.04 


0.05 


0.26 


0.00 


-0.06 


-0.01 


0.00 


0.00 


10 


57.62 


4.10 


-2.25 


0.18 


0.73 


0.01 


-0.22 


-0.04 


0.00 


0.00 


25 


48.89 


8.92 


-5.15 


0.73 


2.66 


0.11 


-0.82 


-0.24 


0.04 


0.02 


50 


38.67 


11.60 


-8.57 


1.76 


6.16 


0.34 


-1.68 


-0.71 


0.16 


0.11 


100 


25.48 


9.16 


-13.66 


3.75 


11.42 


0.81 


-2.55 


-1.50 


0.42 


0.43 


150 


16.34 


4.08 


-17.89 


5.33 


14.31 


1.23 


-2.82 


-1.98 


0.67 


0.88 


200 


8.31 


-1.18 


-21.69 


6.63 


15.99 


1.58 


-2.74 


-2.17 


0.89 


1.39 


250 


-0.06 


-6.11 


-25.56 


7.97 


17.19 


1.83 


-2.31 


-2.11 


1.06 


1.86 


300 


-9.24 


-10.79 


-29.44 


9.13 


18.07 


1.95 


-1.63 


-1.82 


1.15 


2.19 



Table F.5: I = 1 np phase shifts (in degrees) by the Idaho N3LO(500) potential [68]. 



Tiab (MeV) 


^^0 


'Pn 


'Pi 




'P2 




£2 








1 


61.95 


0.18 


-0.11 


0.00 


0.02 


0.00 


0.00 


0.00 


0.00 


0.00 


5 


63.34 


1.63 


-0.94 


0.04 


0.25 


0.00 


-0.05 


0.00 


0.00 


0.00 


10 


59.55 


3.66 


-2.06 


0.16 


0.70 


0.01 


-0.18 


-0.03 


0.00 


0.00 


25 


50.28 


8.16 


-4.86 


0.68 


2.56 


0.09 


-0.74 


-0.20 


0.03 


0.02 


50 


39.82 


10.75 


-8.25 


1.71 


5.98 


0.30 


-1.57 


-0.62 


0.14 


0.09 


100 


26.51 


8.37 


-13.35 


3.75 


11.17 


0.75 


-2.46 


-1.38 


0.39 


0.39 


150 


17.35 


3.36 


-17.58 


5.37 


14.05 


1.15 


-2.76 


-1.84 


0.64 


0.84 


200 


9.29 


-1.86 


-21.39 


6.68 


15.72 


1.49 


-2.70 


-2.03 


0.87 


1.33 


250 


0.89 


-6.76 


-25.26 


8.04 


16.92 


1.74 


-2.29 


-1.98 


1.04 


1.80 


300 


-8.31 


-11.40 


-29.17 


9.19 


17.81 


1.87 


-1.63 


-1.70 


1.14 


2.13 
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Tabic F.6: I = np phase shifts (in degrees) by the Idaho N3LO(500) potential [68]. 



Tlat (McV) 


'Pi 


^^1 




ei 








^G3 


£3 




1 


-0.19 


147.76 


-0.01 


0.11 


0.01 


0.00 


0.00 


0.00 


0.00 


0.00 


5 


-1.55 


118.18 


-0.19 


0.67 


0.22 


-0.01 


0.00 


0.00 


0.01 


0.00 


10 


-3.19 


102.58 


-0.68 


1.15 


0.85 


-0.07 


0.00 


0.00 


0.08 


0.01 


25 


-6.67 


80.45 


-2.81 


1.76 


3.72 


-0.43 


-0.03 


-0.05 


0.55 


0.17 


50 


-10.03 


62.28 


-6.41 


2.05 


8.99 


-1.16 


0.05 


-0.26 


1.61 


0.72 


100 


-14.05 


42.31 


-11.94 


2.32 


17.56 


-2.35 


0.89 


-0.95 


3.50 


2.14 


150 


-17.49 


29.79 


-15.83 


2.67 


22.68 


-3.22 


2.32 


-1.77 


4.83 


3.51 


200 


-21.43 


20.63 


-19.01 


3.16 


24.71 


-3.86 


3.66 


-2.56 


5.66 


4.62 


250 


-26.08 


13.33 


-21.94 


3.72 


24.29 


-4.27 


4.51 


-3.19 


6.02 


5.35 


300 


-31.28 


7.24 


-24.48 


4.02 


22.00 


-4.41 


5.35 


-3.53 


5.91 


5.57 
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